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Abstract 

This  thesis  investigates  the  feasibility  of  a  proposed  hybrid  linear/FVissy  cmitroller  for  non¬ 
linear  plants.  The  proposed  contrdler  concept  is  based  on  the  nse  of  multiple  linearisations  of  a 
nonlinear  plant,  which  describe  the  dynamics  of  pertnrbations  about  equilibrium  points  throughout 
the  desired  envelope  of  operation.  A  bank  of  linear  compensators  is  devdoped,  each  cmresponding 
to  a  linearized  plant  about  a  different  equilibrium.  The  multiple  contitd  signals  generated  by  the 
bank  of  compensators  are  then  weighted  and  summed  using  FSizzy  Lofpc  to  produce  a  comporite 
control  perturbation  signal,  which  is  used  to  drive  the  nonlinear  plant. 

Experiments  were  conducted  to  test  and  refine  this  control  approach.  Anslysis  shows  that  a 
linear/E\izzy  compensator  based  only  on  a  bank  of  linear  compensators  is  not  fearible,  largdy  due 
to  the  small  regions  for  which  the  linearized  modds  were  valid  and  energy  consideratkms  within 
the  plant/controller  system.  The  analysis  itself,  however,  suggests  an  alternate  fmm  for  a  hybrid 
Linear/Fbzzy  approach,  based  on  linear  waveforms  tuned  by  FHizzy  Log^c.  This  conc^t  is  devdoped 
into  the  Modd-Based  Fuzzy  Logic  Controller  (MBFLC). 


INVESTIGATION  INTO  MODEL-BASED  FUZZY  LOGIC  CONTROL 


/.  Introduction 

A  wealth  of  information  is  available  on  the  application  of  contr<d  theory  to  linev ,  time  in¬ 
variant  (LTI)  plants.  Unfortunately,  most  processes  of  interest  in  nature  are  nonlinear.  Though 
optimal  control  techniques  such  as  dynamic  programming  can  be  used  to  derive  optimal  state  feed¬ 
back  controllers  for  nonlinear  systems,  these  algorithms  are  nnm^cally  intensive  and,  in  general, 
don’t  lead  to  closed  form  solutions  [1].  Approximations  such  as  ’’linearization”  and  the  ”small  per¬ 
turbations  hypothesis”  must  therefore  be  invoked  to  enable  8q>plication  of  conventional  LTI  contnd 
methods  to  these  problems.  The  control  designs  based  on  this  linear  hypothesis  will  then  be  valid 
for  the  nonlinear  process  to  the  extent  that  the  linearity  assumptions  are  not  significantly  viidated. 

The  approximations  required  to  obtain  a  linear  problem  formulation  lead  to  dissimilarities 
between  the  plant  dynamics  and  the  dynamics  of  the  modd.  The  resulting  modeling  error  is 
often  neglected  in  engineering  practice  because  the  mathematical  methods  required  to  address  this 
uncertainty  are  quite  sophisticated  [2].  In  applications  such  as  aircraft  flight  contnd,  these  n^ected 
dynaunics  can  become  significant,  resulting  in  system  behavior  not  predicted  by  the  linear  modd, 
as  alluded  to  in  [3:pages  17-26].  Consequently,  a  controUer  based  on  the  linear  modd  will  become 
less  effective  as  the  flight  condition  of  the  aircraft  moves  away  from  the  equilibrium  condition  about 
which  the  plant  dynamics  were  linearized.  Even  when  operating  dose  to  the  linearization  point, 
the  actual  system  will  sddom  exhibit  the  performance  predicted  by  linear  analysis,  espedally  over 
long  periods  of  time  relate  to  the  natural  frequency  of  the  LTI  modd[4:18]. 

One  possible  means  of  addressing  the  inevitable  disparity  between  the  nonlinear  physical 
plant  and  the  linearized(LTI)  plant  mojfl  is  through  the  use  of  Dizzy  Lo^c.  Dizzy  Logic  is  a 
partial  membership  set  theory  developed  by  Lotfi  Zadeh  in  the  mid-1960’s.  Ori^nally  concdved 
for  machine  learning  applications.  Dizzy  Lo^c  is  now  bdng  used  in  many  Adds  of  academic  and 
commerdal  endeavor,  including  control  systems  [5].  Dizzy  Lo^c  is,  in  essence,  a  means  of  repre¬ 
senting  uncertsdnty  about  a  ^ven  system  or  process  without  directly  applying  statistical  methods 
[6]. 
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FVizzy  Logic  has  been  successfully  applied  to  control  system  design,  and  represents  one  means 
of  developing  controllers  for  nonlinear  plants  [7].  reported  applications  of  Fuzzy  Logic  Contrrd 
(FLC)  include  automotive  control  (breaking,  automatic  parking,  automatic  transmissions,  exhaust 
missions  control)  (8, 9],  industrial  automation  [10],  subway  train  braking  [11],  and  flight  control  [12, 
14].  Current  FLC  work  focuses  primarily  on  the  automatic  control  of  plants  traditionally  operated 
by  humans.  Though  nonlinear,  these  plants  exhibit  simple,  consistent  dynamics  with  a  sense  of 
’’directionality”  associated  with  control  actions.  Therefore,  the  effect  of  a  given  contr<fl  action  is 
known,  at  least  qualitatively,  without  the  aid  of  a  system  model.  Even  in  cases  where  models 
are  available,  they  are  largely  overlooked  in  current  FLC  devdopment,  inhibiting  the  application 
of  FLC  to  complex  plants.  It  is  reasonable  to  assume  that  a  hybrid  conventional-Fuzzy  contnfl 
technique,  dubbed  model-based  Fuzzy  control  (MBFLC)  for  this  thesis,  may  provide  better  overall 
control  than  either  control  approach  alone.  This  thesis  addresses  the  integration  of  FLC  concepts 
and  techniques  into  traditional  model-based  controller  design  methods. 

This  research  is  motivated  by  the  dedre  to  devdop  a  full-envelope  aircraft  flight  controUo^. 
Current  flight  controllers  rely  on  gain  scheduling  to  compensate  for  the  nonlinear  aircraft  operating 
in  different  regions  of  the  (uz,oref)  flight  envelope.  MBFLC  represents  one  alternative  approach  to 
gmn  scheduling,  especially  in  flight  regimes  such  as  high  angle  of  attack  flight  where  sensor  inputs, 
and  thus  full  knowledge  of  u  and  a,  are  often  noise-corrupted.  Though  this  work  is  intended  to  lay 
the  groundwork  for  a  full-envelope  flight  controller,  the  results  are  equally  applicable  to  any  system 
for  which  the  agreement  between  physical  plant  and  linearized  model  is  uncertain. 

In  tlus  thesis  report,  hybrid  nonlinear/Fuzzy  and  linear/Fuzzy  controllers  will  be  devdoped 
based  on  Unearizations  of  a  highly  nonlinear  plant.  The  result -ng  hybrid  controllers  will  then  be 
applied  to  tracking  and  regulation  of  the  nonlinear  plant,  rather  than  linearizations  assumed  to 
represent  the  nonlinear  plant.  F\izzy  Logic,  in  this  context,  will  be  used  to  bridge  the  gap  between 
linear  systems  theory  and  nonlinear  control  application.  Said  another  way,  the  language  of  Fuzzy 
Logic  will  be  used  to  quantify  the  way  in  which  linear  controllers  eventually  fail  when  applied  to  the 
control  of  nonlinear  plants.  This  will  then  be  incorporated  into  an  otherwise  deterministic  controller 
structure.  Since  even  the  simplest  MBFLC  controller  will  be  more  complex  than  a  conventional 
linear  controller,  a  commensurate  increase  in  performance,  robustness,  and/or  envelope  of  operation 
is  to  be  expected. 
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This  thesis  is  divided  into  six  chapters: 

1.  Introduction  •  This  chapter  outlines  the  problem  being  addressed,  the  current  state  of  F^uzy 
Logic  C!oatnd,  any  necessary  assumptions,  and  the  proposed  scdution.  The  introduction 
de^<^  the  justification  and  conceptual  basis  for  Model- Based  F\izzy  Logic  Ck>ntrollers. 

2.  Problem  Development  and  Current  Approaches  -  In  this  section  the  criteria  for  a  suitable 
nonlinear  test  plant  and  its  performance  specifications  is  pven.  The  final  form  of  the  non¬ 
linear  plant  will  be  developed  as  well  as  the  linearized  plant  model.  To  gain  control  insights, 
compensators  will  be  devdoped  for  the  nonlinear  plant  using  dynamic  inversion.  Quantita¬ 
tive  Feedback  Theory(QFT),  and  FLC  concepts.  These  designs  will  serve  as  benchmarks  to 
measure  the  true  effectiveness  of  Model- Based  Frizzy  Logic  Contnd. 

3.  Linear  Design  Considerations  -  In  this  chapter,  a  linear  compensator  is  devdoped  for  the 
linearized  version  of  the  nonlinear  plant.  A  linear  time-var3ring  version  of  the  nonlinear 
plant  is  also  devdoped  and  justified.  The  feasibility  and  desirability  of  using  banks  of  linear 
compensators  to  control  linear  and  linear  time-varying  plants  is  explored.  Error  sources 
inherent  in  banked  compensator  designs  applied  to  linear  and  linear  time-varying  plants  are 
identified. 

4.  Nonlinear  Design  Considerations  -  The  effectiveness  of  linear  compensators  when  applied  to 
the  nonlinear  plant  is  quantified.  Time-varying  linear  compensation  is  devdoped  to  contrd 
the  time-varying  linear  plant.  A  nonlinear  compensator  incorporating  time-varying  linear 
models  is  proposed  and  validated.  The  effect  of  the  time-varying  linear  compensator  is  then 
approximated  using  Fuzzy  Logic,  thereby  creating  a  hybrid  nonlinear/F^izzy  controller. 

5.  Modd- Based  F\izzy  Logic  Controller  Devdopment  -  The  effect  of  the  nonlinear  dements  of 
the  hybrid  nonlinear/Fuzzy  controller  is  analyzed  in  light  of  F\izzy  Lope  approximation. 
A  linear/F\izzy  structure  is  developed  to  eliminate  the  need  for  nonlinear  processing.  The 
hybrid  linear/Fiizzy  design  is  modified  to  encompass  full-envdope  operation.  This  is  the  final 
configuration  of  the  Modd-Based  FVizzy  Logic  Controller.  The  performance  of  the  MBFLC 
will  be  compared  with  other  design  approaches. 

6.  Conclusions  and  Recommendations  -  This  chapter  contains  an  overall  assessment  of  the 
promise  MBFLC  shows  based  on  the  research  conducted  for  this  thesis.  Unresolved  issues 
and  reconunendations  for  future  work  are  provided. 


1.1  Background 

FLC  has  the  potential  to  overcome  two  obstades  inevitably  confronted  when  designing  control 
systems  for  real-world  plants:  1)  Plant  uncertainty,  and  2)  Plant  nonlinearity.  ESach  of  these 
difficulties  is  addressed  in  detail  below. 
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1.1.1  Plant  Uncertainty.  Consider  a  plant  whose  transfer  function  is: 


y(»)  _  1 

U(s)  rs+1 


(1.1) 


with  the  value  of  r  being  dependent  on  the  value  of  the  output/state  y  (r  =  f{Y)).  This  is 
a  common  situation  where,  for  example,  the  dynamic  pressure  on  an  aircraft  is  a  function  o( 
its  altitude.  Note  that  this  plant  is  nonlinear,  so  expressing  this  plant  model  using  a  transfer 
function  is  an  abuse  of  notation.  However,  by  assuming  r  is  approximatdy  constant,  is  time- 
varying  independent  of  Y(t),  or  is  unknown,  all  allow  for  linear  or  linear  stochastic  descriptions 
of  this  nonlinear  problem.  An  additional  ^proach  to  linear  analysis  could  be  to  assume  that  r  is 
Fuzzy. 

If  we  wished  to  develop  a  controller  for  this  plant  based  on  a  Fhizsy  assumption,  we  could 
consider  r  to  be  a  Ebzzy  variable.  The  information  available  on  r  could  be  as  rudimentary  as: 


s  IF  y  is  Negative  THEN  r  is  SnuUl 
s  IF  y  is  Zew  THEN  r  is  Medium 
s  IF  y  is  Positive  THEN  r  is  Large 

This  would  approximate  the  true  knowledge  that  r  varies  with  Y.  These  rules  make  up  a  fkzzy 
Inference  system. 

The  terms  in  italics  are  to  be  considered  linguistic  labels  for  FHizzy  Sets.  Hence,  the  Y  input 
to  the  controller  will  first  be  fuzzified.  That  is,  the  antecedent  Y  will  be  converted  from  a  crisp 
(single-valued)  variable  to  a  Fuzzy  (multi-valued)  variable.  The  fuzzification  procedure  involves 
the  input  value  Y  being  evaluated  for  membership  in  each  of  the  sets  defined  over  the  universe  of 
discourse  (total  allowed  variation)  of  y  -  in  this  case  Negative,Zero,  and  Positive.  The  degree  to 
which  a  given  Y  will  be  considered  a  member  of  any  Fuzzy  Set  is  pven  by  a  membership  value,  p, 
which  ranges  from  0  (complete  exclusion)  to  1  (complete  membership).  The  value  which  p  takes 
on  is  determined  by  an  appropriate  membership  function  for  the  set  in  question,  p(Y),  defined 
over  the  entire  universe  of  discourse.  Sample  membership  functions  for  the  sets  Negative,  Zero  and 
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Figure  1.1  Evaluation  of  the  input  variable  Y(in)  =  0.75  for  membership  in  the  Puzxy  Sets 
YiaNegative,  YisZero,  and  YisPosUive 

Positive  are  shown  in  Figure  1.1.  Note  in  the  Figure  that  Y (tn)  =  0.75  is  shown  being  evaluated  for 
membership  in  each  of  the  three  defined  Fuzzy  Sets.  Y  =  0.75  is  a  member  of  N^ative  to  ns  =  0.0, 
of  Zero  to  Hz  =  0.062,  and  of  Positive  to  /tp  =  0.745.  These  three  partial  set  memberships  represent 
Y  =  0.75  in  the  Fuzzy  domain.  Similarly,  the  output  variable  r  is  also  described  by  three  fHizzy 
Sets  that  span  its  universe  of  discourse-  r  it  SmcUly  r  is  Medium,  and  r  is  Large,  as  shown  in 
figure  1.2.  In  this  case,  membership  functions  relate  the  FVizzy  Sets  for  r  to  a  range  possible 
crisp  r  values.  The  membership  functions  defining  the  F\izzy  Sets  for  r  will  not,  in  general,  resemble 
the  membership  functions  for  the  fuzzification  of  Y.  The  Fuzzy  Implicate  r  is  related  to  the  fHizzy 
Variable  Y  through  the  FVizzy  Inference  system  ^ven  above.  Because  the  fVizzy  Inference  rules 
contain  only  a  single  premise  (for  this  example,  "IF  Y  is  Negative”),  the  FVizzy  Set  representing 
the  resultant  will  be  activated  to  the  same  degree  as  its  premise.  This  is  illustrated  in  the  Figure. 
Therefore,  Hn  —  0*0  Hs  =  0-0,  Hz  —  0.062  — »  hm  =  0.062,  ftp  =  0.745  Hl  —  0.745. 
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Figure  1.2  Activation  of  Implicate  Tnxzy  Sets  based  on  membership  values  of  premise  statements 

A  ’’cri8p’’(single  valued)  estimate  of  r  can  be  produced  by  any  of  a  number  of  defuzzification 
routines  [13].  Generally,  these  algorithms  involve  taking  some  sort  of  centroid  of  the  p<dygon  formed 
by  all  of  the  activated  re^ons  of  the  implicate  sets,  as  shown  in  the  Figure.  In  this  manner,  for 
any  crisp  input  Y,  a  crisp  output  r  will  be  produced.  A  compensator  of  the  form  U  =  F{E,t) 
can  then  be  developed  for  the  plant  in  (1.1)  based  on  estimates  of  r  supplied  by  the  Faxxy  Logic. 
This  represents  the  first  MBFLC  scheme  for  feedback  control,  and  is  shown  in  Figure  1.3.  For  an 
overview  of  INizzy  Sets  and  f\izzy  Membership  FVmctions,  refer  to  [7]  or  [15]. 

Of  course,  the  true  utility  of  this  Fuzzy  approach  is  realized  only  when  the  relation  between 
r  and  Y  is  poorly  understood  or  variable.  In  the  case  where  r  =  f{Y)  is  known  or  statistically 
modeled,  other  approaches  will  probably  yield  better  performance,  though  at  higher  computational 
expense.  The  objective  of  this  type  of  hybrid  FNizzy  design  would  be,  not  to  obtain  optimum 
performance  at  a  single  r  trajectory,  but  to  obtain  ’’acceptable”  performance  over  a  broad  range 
of  T  trajectories. 
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Figure  1.3  Block  diagram  showing  closed-loop  compensation  of  plant  from  Equation  1.1  using 
FVizzy  Lo^c  to  provide  an  ’’estimate”  of  r 

The  Fuzzy  Logic  element  in  Figure  1.3  serves  essentially  as  a  nonlinear  look-up  table,  with 
the  qualities  of  the  table  alterable  in  real-time.  If  the  Fuzzy  Logic  estimate  of  r  was  perfectly 
accurate,  then  a  system  response  very  similar  to  the  established  linear  specifications  is  possible.  It 
should  be  noted,  however,  that  even  for  a  known  r  history,  stability  is  not  guaranteed  by  using  a 
linear  compensator  based  on  r  [16].  As  the  estimates  are  d^aded,  closed-loop  system  performance 
degrades  as  well.  Simulation  and  on-line  retuning  would  be  required  to  maximize  performance. 

This  MBFLC  approach  offers  several  benefits  over  conventional  techniques  (parameter  esti¬ 
mation  or  robust  control).  The  controller  is  conceptually  simpler  and  quicker  to  synthesize  than 
other  control  laws,  as  is  the  case  with  most  Fuzzy  Logic  Controllers  (18).  Lessons  learned  from 
testing  can  be  directly  included  in  the  control  law.  The  flexible  structure  of  the  MBFLC  could 
make  the  system  less  susceptible  to  external  disturbances,  and  some  types  of  system  ndse  will  be 
filtered  from  the  feedback  channel  [19].  Finally,  Model-Based  Fuzzy  Lo^c  Contitd,  as  with  other 
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Fuzzy  Lojpc  control  methodologies,  is  well  suited  to  real-time  adaptive  control.  The  drawbacks  to 
MBFLC  include  lack  of  stability  guarantees,  lack  of  a  priori  performance  analysis  techniques,  tuning 
requirements,  and  the  requirement  for  an  on-line  processor  to  perform  the  Fuzzy  Lo^c  operations. 

It  is  reasonable  to  ask  when  it  would  be  beneficial  to  consider  MBFLC  over  statistically 
based  techniques.  The  answer  is  that  the  benefit  of  using  a  MBFLC  is  directly  related  to  the 
amount  of  non-statistical  or  unmodeled  uncertainty  present  in  the  system.  As  will  be  shown  in 
the  next  section,  the  non-statistical  uncertainty  present  in  the  unknown  parameter  r  is  directly 
related  to  the  nonlinear  model  from  which  the  linearized  system  was  derived.  The  idea  that  both 
parameter  uncertainty  due  to  unmodeled  nonlinearities  and  uncertainty  due  to  linearization  of 
nonlinear  plants  can  be  expressed  as  time-varying  parameters  in  otherwise  linear  systems  will  serve 
as  the  justification  for  MBFLC. 

1.1.2  Plant  Nonlinearities.  Consider  the  general  control  problem,  as  it  applies  to  a 
nonlinear  plant.  Assume  a  model  of  the  plant  dynamics-however  complicated-is  available.  It  can 
be  expressed  in  the  form: 


k  =  f{X,U)  (1.2) 

where  2C.  €  92"  is  a  vector  of  system  states,  U_  €  92"*  is  a  vector  of  system  inputs  (controls),  and 
/  is  a  smooth  nonlinear  function  of  2C  ^d  IJ_.  It  is  desired  to  drive  the  system  from  some  initial 
state,  2Cay  Hoi  to  some  final  state,2C/?  U.jt  within  a  fpven  maximum  settling  time  t,  and  with  no 
state  Xi{t)  overshooting  it  final  desired  value  Xkf{i)  by  more  than  a  ratio  M,. 

This  problem,  as  formulated,  cannot  be  directly  addressed  using  LTI  control  theory.  How¬ 
ever,  a  linearized  model  of  the  plant  dynamics  at  certain  points  in  state  space  can  be  derived  by 
expressing  the  state  and  control  vectors  in  terms  of  perturbations  from  a  ^ven  nominal  condition, 
XrS.  [20:pages  83-97].  The  state  and  control  vectors  become: 

K  =  ^  +  (1.3) 
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H  =  ^  +  M 


(1.4) 


where  ^  and  are  determined  by  setting: 

/(li.fii)  =  0.  (1.5) 

That  is,  the  ”tth  nominal  point  ^  is  an  equilibrium  or  trim  point.  The  subscript  t,  t  = 
l,2,...,i\^,  indicates  that  there  could  be  msmy  such  nominal  points.  The  plot  of  these  equilibria 
points  in  an  nzm-dimensional  space  defined  by  the  ranges  of  ^d  U_  is  referred  to  in  this  thesis 
as  a  nominal  locus. 

Inserting  the  expressions  (1.3)  and  (1.4)  into  (2.36)  and  expanding  the  nonlinear  /  function 
by  a  Taylor  series  about  (J^j,  Oj)  yields: 

i  +  i  =  +  H.O.T.  (1.6) 

The  partial  derivatives  are  Jacobian  matrices  and  represents  higher-order  terms  in  the 

Taylor  series.  Two  assumptions  must  be  made  at  this  point  in  order  to  obtain  a  linear  model  of 
plant  dynamics  about  the  tth  equilibrium  point: 

1.  Only  small  perturbations  about  the  equilibrium  condition  will  occur  (’’H.O.T”  terms  are  then 
negligible). 

2.  The  equilibrium  point  is  static  (^  =  0). 

By  applying  these  two  assumptions  and  noting  that  /  evaluated  at  the  equilibrium  point  is  zero 
(by  definition),  (1.6)  becomes  a  linear  equation  in  the  perturbation  variables: 

i  =  AiS.  +  BiU.  (1.7) 

This  is  the  familiar  derivation  of  an  LTI  plant  by  linearization.  Note  the  explicit  dependence  of 
the  (i4,  B)  plant  model  on  the  tth  equilibrium  point.  This  linearization  is  only  valid  at  or  near 
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this  equilibrium  point.  This  model  says  that,  though  the  behavior  of  the  system  is  nonlinear,  the 
movement  of  the  states  relative  to  am  equilibrium  pmnt  can  be  modeled  as  linear,  provided  the 
actual  states  amd  controls  au’e  ”€1086  enough”  to  equilibrium.  At  what  distamce  from  the  nominal 
the  Uneau’  model  no  longer  holds  depends  on  the  extent  of  the  originad  nonlinearities. 

When  designing  a  control  system  for  a  nonlineatr  plamt  such  as  an  aurcradt,  it  is  customary 
to  consider  the  performance  of  the  plant  at  mamy  equilibrium  points.  In  fact,  gain  scheduling 
in  flight  control  systems  is  required  because  aurcradt  exhibit  such  varied  dynamics  about  different 
equilibria  within  their  flight  envelope  [21].  However,  by  observing  that  the  dynamics  change  in 
a  smooth,  continuous  mamner  throughout  the  flight  envelope,  it  becomes  cleatr  that  the  aurcradt 
transfer  function  may  be  considered  linear  but  containing  pau'auneters  varying  as  a  function  of  the 
current  state,  as  in  the  case  of  pau'aimeter  uncertaiinty  discussed  in  Section  1.1.1.  The  linearization 
process  introduces  signiiicamt  non-statistical  uncertainty  into  the  unknown  pa^aimeters  [22]. 

Other  sources  also  contribute  to  the  non-statisticad  uncertainty  aissociated  with  the  amalysis 
of  a  nonlineau^  plant  such  as  am  aircraft.  First,  all  higher  order  dynamics  au'e  assumed  negligible. 
These  terms  aure  only  negligible  in  actuality  when  the  system  is  exatctly  at  the  nominal  condition. 
Second,  though  the  model  is  now  expressed  in  terms  of  perturbations  from  the  nominad,  all  inconung 
measurements  will  be  of  the  full  state  X  and  the  full  control  U_.  The  nominad  conditions  will  have  to 
be  subtrau:ted  from  eau;h  measurement,  and  atny  error  in  the  nominal  condition  will  directly  manifest 
itself  as  a  disturbamce  driving  the  lineau-ized  system  model  [14].  These  errors  will  be  compounded 
by  amy  aictual  pau'auneter  variation  exhibited  by  the  physical  plamt  itself,  as  well  as  by  noise  present 
in  the  signal  paths  [2]. 

As  developed  in  the  previous  two  sections  amd  elsewhere  [4:688-690],  most  error  sources  cam  be 
accounted  for  by  assuming  parameter  and  signad  uncertainty  in  a  lineair  system  model.  The  current 
practice  is  to  ignore  these  uncertadnties  or  to  include  them  by  aissuming  a  lineau^  system  driven  by 
statisticadly  modeled  process  amd  measurement  noises  [4,  2].  In  both  cases  the  controller  is  then 
derived  in  a  lineatr  fashion.  However,  by  defining  a  "Fuzzy  linear  system,”  where  parameters  are 
assumed  to  have  a  certain  non-statisticad  uncertainty,  the  effects  of  linearization  cam  be  explidtly 
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included  in  the  model.  Though  the  idea  of  a  time-varying  linear  system  represented  by  a  Li^lace 
IVansform  is  an  abuse  of  notation,  it  nonetheless  provides  valuable  insight  into  the  operation  of  the 
proposed  Model-Based  Fuzzy  Logic  Compensator. 
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l.t  Proposed  Approach 

Control  systems  based  on  F\izzy  Logic  have  been  developed  and  implemented  for  many  plants 
in  science  and  industry  [5].  Current  EVizzy  controllers  are  based  on  a  mapping  of  the  input/output 
relations  for  a  given  plant.  These  relations  are  formulated  as  a  series  of  IF...THEN-type  rules  which 
comprise  the  Fuzzy  Inference  system  within  the  Fuzzy  Controller.  Fuzzy  Controller  devdopment 
is  therefore  contingent  on  the  availability  of  a  system  control  knowledge  base,  obtained  either  from 
interviews  with  human  operators  or  &om  tests  on  the  actual  hardware.  The  concept  of  a  Fuzzy 
Controller  based  on  a  model  of  the  system  behavior  rather  than  expert  "rules  of  thumb"  has  not 
been  addressed  in  the  literature. 

To  better  understand  the  potential  benefits  of  model-based  Fbzzy  compensation,  consider  the 
control  design  problem  in  graphical  form,  shown  in  Figure  1.4.  This  figure  shows  the  parameter 
space  for  a  scalar  version  of  the  nonlinear  system  discussed  in  Section  1.1.2.  In  this  diagram, 
the  X  and  u  axes  of  the  graph  represent  the  allowable  rang^  of  variation  of  the  state  and  the 
input  respectively.  The  objective,  as  stated  above,  is  to  drive  the  plant  from  an  initial  equilibrium 
state,  assumed  X  =  0,  U  =  0  in  the  Figure,  to  a  final  value  which  is  shown  at  the  end  of  the 
A'(f),f/(t)-Optimal  trajectory  curve.  If  a  system  were  to  traverse  the  optimal  trajectory,  it  would 
be  assured  of  reaching  the  final  state  within  desired  specifications.  Unfortunately,  this  trajectory 
is  a  mathematical  abstraction  and  is  not  known  prior  to  simulation.  Only  the  nonlinear  equations 
describing  the  plant  are  available,  as  well  as  N  linearized  plants  describing  the  dynamics  about  N 
"evenly"  spaced  operating  points  on  the  locus  of  equilibrium  points.  The  objective  of  the  contrdler 
is  to  remain  approximately  on  the  optimal  trajectory  by  using  the  dynamics  information  available 
in  the  form  of  the  linearized  plants. 

Given  that  the  system  is  at  the  point  labeled  "Current  Position"  in  the  Figure,  the  question  is: 
which  linearized  model  yields  the  best  approximation  of  the  dynamics  at  the  current  point?  Plant 
3  yields  a  better  model  of  the  U_  dynamics  while  Plant  2  yields  a  better  model  of  the  X.  dynamics  at 
the  point  in  question.  Therefore,  without  resorting  to  a  nonlinear  compensator,  the  best  controller 
is  based  on  a  weighted  combination  of  these  two  plant  models.  An  implementation  of  this  approach 
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is  shown  in  Figure  1.5.  Here,  a  linear  contrdler  is  developed  for  each  of  the  N  linearized  plants. 
E^ach  compensator  provides  adequate  closed-loop  performance  for  state  trajectories  originating  and 
terminating  ’’sufficiently  close”  to  the  equilibrium  about  which  the  compensator  is  based.  The 
value  corresponding  to  "sufficiently  close”  depends  on  the  extent  of  the  nonlinearities  exhibited  by 
the  plant  and  will  be  referred  to  in  this  report  as  the  region  of  attraction  of  that  compensator.  The 
N  compensators  are  all  implemented  in  the  forward  path  of  the  closed-loop  system,  and  all  receive 
the  error  signal  E.  Which  compen8ator(s)  wiU  actually  drive  the  plant,  however,  is  determined  by 
the  Fuzzy  Supervisor.  The  Fuzzy  Supervisor  will  appropriately  weight  each  compensator  based  on 
a  Fuzzy  assessment  of  system’s  position  in  the  state  space  of  Figure  1.4.  The  FNizzy  Supervisee  wiU 
base  its  decision  on  a  measurement  of  the  current  system  state  Y,  as  shown  in  the  Figure. 

Each  linearized  plant  is  represented  by  a  PHizzy  Set  within  the  Fuzzy  Supervisor.  When  a 
measurement  of  the  current  state  is  input,  the  Supervisor  will  assign  a  membership  value  pi  for  that 
input  to  each  Fuzzy  Set.  The  values  of  p  are  based  on  membership  functions  Pi{Y)  which  represent 


1-13 


Figure  1.5  Weighting  of  a  bank  of  linear  compensators  using  Fuzzy  Sets 

the  "rules  of  thumb”  for  properly  wdghting  the  compensators.  In  this  research,  the  membership 
functions  are  Gaussian-like  and  are  based  on  the  r^ons  of  attraction  of  the  compensators. 

This  simple  control  approach  leaves  many  questions  unanswered,  particularly  for  multi-state 
problems.  First,  there  are  the  linear  questions  associated  with  the  use  of  multiple  compensators  in 
the  forward  path  of  the  .ilxn^dation.  For  example,  the  fact  that  each  compensator  is  valid  within 
its  region  of  attraction  does  not  necessarily  imply  proper  dynamic  contrd  of  a  plant  traversing 
multiple  r^ons  of  attraction.  Even  the  location  of  the  weighting  blocks  (whether  before  or  after 
the  compensator)  will  determine  the  distribution  of  energy  throughout  the  system.  This  will  have 
a  bearing  on  the  transition  from  one  compensator  to  another. 
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Second,  there  are  also  nonlinear  considerations  to  be  addressed.  The  most  pressing  of  these 
is  that  some  appropriate  control  action  must  be  taken  even  when  none  of  the  compmuators  are 
considered  valid.  This  will  be  a  common  occurrence,  especially  for  multi-state  plants,  so  how  this 
situation  is  handled  will  greatly  impact  the  complexity  and  performance  of  the  final  solution.  The 
designer  would  like  the  controller  to  respond  with  caution.  When  the  the  F\izzy  Supervisor  finds  no 
close  correspondence  between  the  estimated  system  state  and  the  EHizzy  sets  correspcmding  to  the 
universe  of  discourse,  the  compensator  should  decrease  the  input  to  the  plant  (Caution  and  other 
characteristics  exhibited  by  optimal  controllers  are  discussed  in  Maybeck  [1:229]). 

The  above  considerations  and  others  will  be  used  to  determine  the  feasibility  of  a  hybrid 
linear/Fuzzy  Logic  controller  based  only  on  linear  compensators  and  FNizzy  Sets.  Analysis  in  the 
chapters  to  follow  show  that  this  configuration  is  not  adequate  to  control  a  nonlinear  plant.  The 
analysis  itself,  however,  suggests  an  alternate  controller  configuration  based  on  linear-like  waveforms 
generated  using  Fuzzy  Logic.  This  alternative  design  for  the  Model-Based  Fuzzy  Logic  Controller 
is  shown  in  Figure  1.6. 

Note  that  this  design  has  only  a  single  linear  compensator,  which  corresponds  to  the  highest 
gain  required  for  any  linearization  under  consideration.  It  was  found  that  the  linear  compensators 
intermediate  to  the  initial  state  and  final  state  of  the  system  contributed  nothing  to  the  successful 
control  of  the  nonlinear  plant.  Therefore,  only  starting  and  terminating  compensators  are  necessary. 
Between  the  two  regions  of  compensator  validity,  a  Fuzzy  ’’squashing  function”  is  required  to  avoid 
overdriving  the  nonlinear  plant. 

Of  the  two  "compensators,”  only  one  is  implemented  as  a  linear  transfer  function  (shown  at 
the  top  in  the  Figure).  The  second  is  a  F\izzy  Logic  emulation  of  the  u(t)  output  which  would 
have  been  produced  by  a  second  linear  compensator  corresponding  to  the  final  desired  state  of  the 
system.  This  "compensator”  is  implemented  as  shown  in  the  lower  portion  of  Figure  1.6.  This 
Fuzzy  Logic  approximation  is  necessary  because  the  required  form  of  the  second  compensator  does 
not  correspond  to  the  linear  compensator  based  on  the  final  equilibrium  state  of  the  system.  The 
required  controller  output,  then,  must  be  determined  through  simulation. 
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Figure  1.6  Final  structure  of  Model- Based  F\uzy  Logic  Contrcdler 


Tbjs  alternate  Model-Based  Fliszy  Logic  Contrdler  induces  Ilnear-like  closed-loop  perfor¬ 
mance  from  a  nonlinear  plant  over  a  specified  envelope  of  operation,  but  at  the  expense  of  most 
linear  compents  in  the  originally  proposed  controller.  The  internal  operation  of  this  concept  con¬ 
troller  is  developed  further  in  the  chapters  to  follow. 
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t.S  Scope  and  Assumptions 

As  this  study  invcdves  both  uon-Iinear  plants  and  non-linear  compensaUHs,  the  issues  of 
stability,  performance,  and  steady-state  error  will  only  be  addressed  through  experimmitation  via 
simulation.  Rigorous  exploration  of  these  topics  is  beyond  the  scope  of  this  study  and,  in  some 
cases,  beyond  the  state  of  current  nonlinear  mathematics  [21,  27]. 

This  study  is  exploratory  in  nature  so  it  cannot  be  exhaustive  in  its  inclusion  of  all  possible 
nonlinear  plants.  The  objective  is  to  address  plants  pertinent  to  aerospace  science.  The  stipulation 
that  /  be  smooth  is  necessary  to  the  description  of  nonlinear  systems  as  linear  perturbation  modds 
with  varying  parameters.  /  need  not  be  smooth  in  general,  though  the  applicability  of  MBFLC  to 
plants  with  discontinuities  within  the  range  o(  operation  is  beyond  the  scope  of  this  study. 

The  design  approach  developed  here  is  germane  to  both  uncertainty  due  to  nonlinearity 
and  uncertainty  due  to  parameter  variation  or  mismodeling.  The  controller  devdoped  herdn  is 
based  only  on  the  former,  though  the  performance  is  also  measured  against  an  uncertain  modding 
parameter.  Encompassing  both  non-statisticalmor  sources  simultaneoudy  is  left  for  later  research. 

Broadly  stated,  the  objective  of  this  research  is  to  devdop  a  full-envdope  contrdler  for  a 
nonlinear  plant  using  only  linear  control  techniques  and  FHizzy  Lo^c.  The  expected  benefits  of  the 
Model- Based  Rizzy  Logic  Control  approach  are: 

1..  Force  a  nonlinear  plant  to  to  exhibit  linear-like  performance  and  conform  to  linear  design 

spedfications. 

2.  Incorporate  available  modds  into  contrdler  structure. 

3.  Enhance  robustness  of  controller  in  the  face  of  unmodded  uncertainties  (system  damage, 

noise). 

4.  Obtain  full  envdope  operation  without  the  need  for  gain  scheduling. 

5.  Introduce  some  a  priori  synthesis  and  analysis  capability  into  FVizzy  Contrdler  Design. 

1.4  Summary 

Engineering  models  of  real-world  systems  are  inevitably  abstractions  of  reality.  By  modding 
the  extent  of  this  abstraction,  a  more  accurate  representation  of  the  physical  system  can  be  obtained. 
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Fuzzy  Logic  is  one  means  by  which  this  non-statistical  uncertainty  can  be  taken  into  account  in  a 
mathematically  tractable  way. 

Flight  controllers  based  on  Fuzzy  Logic  have  been  developed  and  presented  in  the  open 
literature  [19,  12].  These  designs,  however,  are  based  on  the  assumption  that  no  adequate  plant 
models  exist  for  aircraft.  This  is  not  the  case,  though  the  models  are  nonlinear  and  exhibit  complex 
dynamics.  The  objective  of  this  thesis  is  to  explore  the  application  of  Fdzzy  Set  theory  in  an 
environment  where  the  systems  of  interest  are  described  as  essentially  linear  with  non-statistical 
parameter  variation. 

The  result  of  this  research  is  a  hybrid  linear/Fhzzy  controller  which  forces  a  nonlinear  plant 
to  exhibit  linear-like  closed-loop  behavior. 
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11.  Plant  Development  and  Current  Approaches 


t.l  Introduction 

The  goal  of  this  research  area  is  to  explore  intelligent  flight  control,  where  contrdlers  atilise 
concepts  from  both  the  contnd  and  artificial  intdligence  fields.  It  is  hoped  that  these  hybrid 
contrdlers  will  demonstrate  some  advantages  over  current  controllers,  if  not  in  performance  thra 
in  areas  of  robustness,  noise  rejection,  and  ease  of  design.  The  objective  of  this  thesis  effort  is  to 
demonstrate  the  value  of  applying  both  FHizzy  Set  theory  and  linear  systems  theory  to  the  contnd 
nonlinear  plants.  It  is  conjectured  that  the  resulting  contrdler  will  be  more  successful  than  eitha 
F\izzy  Logic  Control  or  Linear  Contixd  alone. 

Because  the  final  envisioned  application  is  automatic  flight  contnd,  it  would  be  reasrmable  to 
develop  controllers  based  on  the  nonlinear  aircraft  equations  of  motion.  However,  this  theas  effort 
was  conducted  using  a  reduced-order  plant  which  is  subject  to  parameter  variation.  The  reasons 
for  this  are  three-fold: 

1.  AvaUability  of  Nonlinear  Models  -  The  use  of  linearized  aircraft  modds  about  a  given  equi¬ 
librium  condition  is  so  pervasive  that  few  simulations  are  available  to  adequatdy  solve  the 
nonlinear  equations  of  motion  in  real-time.  Developing  the  hardware  and  software  interfaces 
required  for  this  type  of  simulation  is  prohibitive. 

2.  Complexity  of  Linear  Controllers  -  The  linearized  models  resulting  from  the  nonlinear  aircraft 
equations  of  motion  would  have,  as  a  worst  case,  eighth-order  characteristic  equatims.  The 
linear  controller  required  to  adequately  compensate  this  plant  would  be  very  complex,  and 
simulations  would  be  ponderous. 

3.  Difficulty  in  Fuzzy  Tuning  -  The  number  of  fhzzy  parameters  present  in  the  FNizzy  controller 
for  a  full-envelope  flight  control  system  would  be  large.  Without  some  predetermined  rules 
for  tuning,  and  perhaps  even  pre-tuned  components,  the  challenge  of  holistic  tuning  could  be 
insurmountable. 

Because  of  these  three  considerations,  this  research  was  conducted  using  a  two-state  nonUnear 
plant.  The  chosen  plant  is  developed  in  this  chapter.  Conventional  ’’benchmark”  approaches  to 
the  control  of  this  nonlinear  plant  will  also  be  discussed.  All  desi^s  are  validated  using  the 
SIMULINK®simulation  environment  for  MATLAB®. 
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f.S  Redwxd-Order  Plant 


The  ftc^owing  criteria  were  used  in  deciding  upon  the  form  of  the  nonlinear  modd  used  in 
this  effort: 

•  The  plant  should  exhibit  significant  nonlinearities  over  the  expected  ranges  of  the  state  and 
input  vectors. 

•  The  plant  should  exhibit  parameter  uncertmnty. 

•  The  plant  should  have  an  infinite  number  of  equilibrium  conditions. 

•  The  locus  of  points  associated  with  equilibria  should  be  smooth  and  monotonic. 

•  The  system  nonlinearities  should  be  odd-symmetric  to  ensure  tracking  of  system  inputs  for 
both  positive  and  negative  input  values  and  preserve  some  features  of  a  linear  plant. 

•  The  form  of  the  linearized  plant  model  should  resemble  the  plants  used  in  other  studies  of 
linear  systems  theory,  to  the  extent  possible. 

•  The  plant  should  be  as  simple  as  possible  to  clearly  demonstrate  the  the  developed  contndler. 
The  nonlinear  plant  chosen  to  fit  these  requirements  is: 


Xi  =  -Xi  +  aX,  (2.1) 

Xi  =  -Xf  +  U^  (2.2) 

Y  =  Xi  (2.3) 

0.5  <  a  <  1.5.  (2.4) 


When  addressing  only  uncertainty  due  to  linearization,  a  is  set  to  1.  The  SIMULINK  modd 
this  plant  is  shown  in  Figure  2.1.  In  subsequent  simulations,  the  nonlinear  plant  is  represented  as 
a  single  block  called  NonlinJ^lant,  with  o  externally  defined. 

i.S  Linearization  of  Nonlinear  Plant 

In  order  to  apply  linear  system  theory  to  this  plant,  a  linearized  modd  of  this  plant  must  be 
considered,  as  discussed  in  Section  1.1.2.  First,  the  states  and  inputs  are  described  in  terms  of  an 
equilibrium  value,  plus  a  perturbation: 


Xi  =  Xi0  +  Xi 


(2.5) 
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Figure  2.1  SIMULINK  Model  of  Nonlinear  Plant  with  Externally  Defined  a  Parameter. 


Xj  =  Xu  +  *1 

(2.6) 

U  =  U,  +  u. 

(2.7) 

The  relationship  between  the  equilibrium  values  Xio,X}o,  and  U, 

right-hand  sides  of  (2.1)  and  (2.2)  equal  to  0: 

can  be  found  by  setting  the 

0  =  — Xi  -f-  Xj 

(2.8) 

0  =  -xf-l-lf* 

(2.9) 

Xt  ^  X, 

^  u*. 


xf 


(2.10) 

(2.11) 


This  implies  that  Xi,  =  X^  =  U,.  Note  that  this  {daat  has  perfect  tracking  even  without  the 
benefit  of  feedback  compensation. 


To  determine  the  linearized  model,  substitute  (2.5), (2.6),  and  (2.7)  into  (2.1)  and  (2.2).  This 
yields: 

+  *i)  =  — ~  *i  +  +  *a  (2-12) 

+  =  -(X,. +  *»)*  + (*?.  +  «)*.  (2.13) 

Expanding  the  cubic  terms  and  rearran^ng: 

Jfi,  +  Xi  =  —Xi,  +  Xje  —  *1  +  (2.14) 

X^o  +  X7  =  -Xu*-3Yi.**i-3Jlfi^J-x?+l7/  +  3(7,**  +  3f?.tt*  +  tt*.  (2.15) 

Applying  the  small  perturbation  assumption  yields: 

+  *i  =  ~  »i  +  *a  (2.16) 

X^  +  x,  =  -JrV-3jr,.**i  +  f7.*  +  3tF,V  (2.17) 

Finally,  by  assuming  the  plant  is  operating  about  a  static  equilibrium  point,  then  =  0 

and  noting  that  the  terms  which  involve  only  equilibrium  quantities  cancel  out,  the  linearized  model 
for  (2.1)  and  (2.2)  is: 

x’l  =  -Xj  +  Xj  (2.18) 

xj  =  -3jr„’x,  +  3l7,’tt  (2.19) 

y  =  Xi.  (2.20) 
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Because  Xx,  =  X^o  =  0,,  a  new  quantity  r  can  be  defined  as  r  =  zXi,  =  3(J,.  Therefore,  the 
linearized  model  becomes: 


*l  = 

-*i  +  *» 

(2.21) 

*J  = 

—TXi  -1-  ru 

(2.22) 

y  = 

*1- 

(2.23) 

As  mentioned  in  Chapter  1,  the  linear  model  is  based  on  two  assumptions:  1)  Small  pertur¬ 
bations  from  the  equilibrium  condition,  and  2)  Static  equilibrium  condition.  As  long  as  a  single 
equilibrium  point  is  chosen,  as  was  the  case  in  developing  the  Banked  Compensators,  r  is  a  con¬ 
stant  and  the  Laplace  domain  transfer  function  for  this  state  space  description  can  be  derived  via 
Gis)  =  C(sl  -  A)-^B  +  D»s: 


P{s)  = 


T 


(2.24) 


Two  points  should  be  mentioned  about  this  plant  selection  which  will  bear  on  the  design  of  a 
successful  controller: 


1.  r  changes  rapidly  relative  to  the  eigenvalues  of  the  linearized  model.  The  eigenvalues  of  the 
linearized  model  are  given  by: 

.  -\±y/\-  12Y^ 

Ai.a  - - ^ - 

Comparison  with  the  expression  for  r  shows  that  the  A  matrix  will  vary  rapidly  with  time 
for  large  step  inputs.  Much  of  the  linear  controller  analysis  wiU  assume  that  the  plant  can  be 
considered  piecewise  LTI,  but  research  has  shown  that  this  is  precisely  the  type  of  plant  which 
frequently  violates  these  assumptions  [16].  Linear  analysis  must  be  viewed  with  caution,  as 
stability  and  performance  are  not  guaranteed.  This  is,  indeed,  a  difficult  nonlinear  control 
problem. 

2.  Plant  is  unstable  for  T  <  0.  As  T  -►  0-f-  the  system  will  become  less  and  less  stable.  This 
will  adversely  impact  control  around  the  ori^n. 


The  SIMULINK  model  of  this  plant  is  shown  in  Figure  2.2.  Note  that  it  is  represented  in  the 
state  space  form  of  (2.21)  and  (2.22).  In  subsequent  simulations,  this  model  will  be  represented  as 
a  single  block  called  LinJ*lant  with  the  parameter  r  externally  defined. 
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Figure  2.2  SIMULINK  model  of  linearized  plant  with  externally  defined  r  parameter 
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2.^  Control  Performance  Characterutiee 

It  is  de6iT«d  to  develop  a  control  system  for  the  nonlinear  plant  of  interest.  The  following 
performance  characteristics  were  arbitrarily  chosen  as  a  basis  for  controller  design: 

•  Peak  Overshoot  (Afp):  1-12  ((  =  0.5594) 

•  MaLximum  Settling  Time  (t,):  1.62  seconds 

•  Disturbance  Rejection:  |y(t)A4titrtmd«t  ^  0.1 

The  desired  response  is  produced  by  the  following  modd: 


Af(s)  = 


19.5  19.5 

+  4.940.5.«!  +  19.5  ~  (j  +  2.47  +  j3.66)(e  +  2.47  -  ;3.66)‘ 


This  response  is  shown  in  Figure  2.3. 


(2.26) 
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TlmtCMOC) 

Figure  2.3  Reeponsc  of  Model  Plant  to  Unit  Step  Input 
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£.5  Standard  Control  Approaches 


Before  considering  the  xynthesis  of  two  control  a.pproacbes,  it  is  iraponaut  to  assess  the 
potential  of  each  technique  alone.  This  analysis  will  not  only  provide  valuable  benchmarks  for 
performance,  but  will  also  provide  insight  into  how  a  MBFLC  compensator  should  operate. 

Three  techniques  were  considered  in  this  research:  Dynamic  Inversion,  Quantitative  Feedback 
Theory,  and  Fuzzy  Logic  Control.  The  performance  of  these  techniques  was  measured  using  the 
SIMULINK^simulation  environment  for  MATLAB®. 

2.5.1  Dynamic  Inversion.  Dynamic  Inversion(DI)  is  a  straightforward  technique  which 
leads  to  nonlinear  controllers  for  many  types  of  nonlinear  plants.  Though  the  objective  of  this  thesis 
is  not  to  address  nonlinear  control  techniques  explicitly,  Fuzzy  Logic  Controllers  axe  themselves 
nonlinear  and  the  operation  of  the  Dl-based  controller  will  serve  as  the  starting  point  for  Fuzzy 
tuning. 

The  objective  of  DI  is  to  reformulate  the  problem  into  a  linear  system  through  the  use  of 
intermediate  quantities.  Linear  systems  theory  can  then  be  used  to  synthesize  a  controller  fot  the 
"linear”  plant.  In  implementation,  the  intermediate  quantities  required  to  obtain  a  linear  system 
model  can  be  solved  for  in  an  ad  hoc  fashion  and  fed  into  the  nonlinear  plant.  For  more  on  Dynamic 
Inversion,  see  Reference  [27]. 

The  nonlinear  plant  of  (2.1)  and  (2.2)  can  be  viewed  as  linear  by  defining  an  intermediate 
control  variable  V  =  -Xf  +  CT®.  This  leads  to: 


Xj  —  —Xi  +  aX2 

(2.27) 

(2.28) 

Y  =  Xu 

(2.29) 

2-9 


II  *94  13:03 


513  476  4055  PftGE . 006 


Flgtm  2.4  Root  Locus  for  G(s)P{a)^  Dynamic  Invernon  Synthesis 
which  is  a  linear  system.  The  transfer  function  is: 

(««) 


It  is  clear  tcom  root  locns  analyius  that  the  desired  poles  of  s  -2.47  +  >3.66,  -2.47  -  >3.66 
cannot  he  obtained  by  gain  adjustment  alone.  It  is  necessary  to  caned  out  the  ^ant  pole  and  s  s  -l 
and  add  a  compensator  pole  farther  to  the  left  in  the  s>plane.  The  proposed  linear  compmisator  is 
therefore: 


(?(s)  = 


^s+1} 

(s  +  5)  ’ 


(2.31) 


The  root  locus  associated  with  this  G(5)/'(s)  is  shown  in  Tignre  2.4.  The  desired  <  =  .5594  is 
shown  on  the  plot  and  can  be  used  to  determine  the  value  for  k.  IVom  this  analysis,  k  s 
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Figure  2.5  SIMULINK  Model  of  G{s)  Using  Dynamic  Inversion 

This  transfer  fancliuu  cau  be  represented  in  state  space  form  using  the  technique  in  [4]  as 
shown  in  Figure  2.5.  This  model  is  referred  to  with  a  single  block  called  DlComp  in  subsequent 
simulations. 

The  controller  described  by  (2.31)  will  accept  the  error  signal  E  as  input  and  will  output 
V,  the  intermediate  control  quantity.  In  order  to  determine  the  correct  input  into  the  plant,  the 
control  U  rmist  be  solved  for  via; 

U  =  i/V^.  (2.32) 

This  equation  will  be  solved  in  a  simulation  block  called  an  interpreter^  thus  yielding  a  nonlinear 
control  law.  Though  the  interpreter  should  only  consist  of  a  sin^e  function  block,  MATLAB  is 
unable  to  perform  the  cubed  root  operation  when  the  argument  is  negative.  Therefore,  the  more 
complex  strncturc  shown  in  Figure  2.t)  was  required. 

The  simulation  of  this  control  system  is  shown  in  Figure  2.7.  The  response  of  the  system  to 
a  unit  step  input  are  shown  in  Figure  2.8.  Also  plotted  in  the  Figure  is  the  response  of  the  model 
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Figure  2.6  Interpreter  Block  Wkich  Converts  V  ioU 

plant.  The  model  respoDse  and  closed'loop  responae  based  on  D1  are  almost  identical.  The  state* 
space  representation  (similar  to  Figure  1.4)  of  the  plant  transfer  from  JTi  s  0,  =  0,  If  s  0  to 

Xi  s  IfU  -  I  is  shown  in  Figure  2.9.  The  nominal  locus  is  also  plotted  on  this  Itgue  as  a 

dashed  line.  Note  that  the  locu^  takes  a  step  up  along  the  If  axis  immediately  upon  onset  of  the  step 
input.  This  moves  the  system  state  away  from  the  nominal  locus,  where  lineaiiaed  approKhaations 
are  valid.  In  three  dimensions,  the  state  trajectory  never  again  approaches  the  eqnilibnnm  bens, 
where  all  compensator  redone  of  attraction  are  centered.  This  is  a  source  of  problems  in  the  banked 
compensator  approach,  which  hinges  on  proximity  to  the  equilibrium  points  for  the  compensators 
to  be  valid. 

The  response  of  the  system  when  the  nonlinear  plant  takes  on  various  values  of  e  is  shown 
in  Fignre  2.10.  The  compensator  is  designed  in  all  cases  for  a  a  i.  This  Figure  shows  that  this 
simulation  is  affected  significantly  though  not  catastroplucally  by  errors  in  modeling  the  nonrincar 
plant. 

Though  dynamic  inversion  is  a  simpb  technique  in  nonlinear  control,  it  will  not  work  in 
many  situations  (27].  A  more  sophisticated  technique  called  progressive  linearization  provides  a 
more  general  solution  to  the  nonlinear  control  problem.  Both  of  these  synthesis  techniques  produce 
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Figure  2.7  Simulation  of  doaed-loop  traddng  coatioUer  devdoped  using  dynunic  inverdoi 

noiilliiear  contzoUen  and  neither  explidtly  addrett  diaturbance  rejection  or  parameter  robeetaeia. 
Distorbance  rejection  can  be  accounted  for  to  lome  degree,  however,  in  the  the  linear  compeSMtor 
design.  A  dynamic  inversion-like  approach  win  serve  as  the  model  for  initial  MBFLC  dostgas. 
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Figure  2.1 0  Response  of  dynamic  inversion  simulation  when  the  nonlinear  plant  takes  on  different 
values  of  the  parameter  a.  a  =  1  in  the  compensator  for  all  simulations. 


2-16 


JhN  11  *94  13:06 


513  476  4055  PAGE. 013 


t.S.t  QmntitativK  Feedback  Theory.  Th«  previous  two  section  addressed  nonltiMar  com¬ 
pensation  for  nonlinear  plants.  Since  tlie  objective  of  this  thesis  is  to  apply  linear  and  Foezy  control 
theory  to  nonlinear  plants,  it  is  important  to  consider  a  purely  linear  techniQue  as  well.  Thoo^ 
linear  system  theory  is  generally  applicable  only  as  long  as  the  linearization  assumptions  hold, 
some  techniques  lead  to  /tneor  controllers  with  some  applicability  to  nonlinear  plants  ow  a  given 
operating  range.  One  such  technique  is  Quantitative  Feedback  Theory  (QFT). 

QFT  was  developed  by  Dr  Isaac  M.  Horowitz  in  the  19703  to  address  parameter  uncertainty 
in  linear  plants  [29].  As  discossed  previously,  the  variation  could  be  due  to  mlsmodtiing,  parametei' 
uncertainty,  or  underlying  nonlinear  djmamics.  QFT  synthesis  yidda  a  linear  compensator  and 
preiilter  which  ensnres  adequate  performance  for  any  one  of  a  family  of  linear  plants.  A  top-level 
view  of  QFT-derived  compensation  is  illustrated  in  Figure  2.11.  Though  the  resulting  dosed-loop 
performance  is  only  guaranteed  for  linear  time«invariant  plants,  the  QFT  compensator  is  often 
capable  of  controlling  even  a  nonlinear  plant  over  a  limited  operating  range. 

The  QFT  design  process  is  initiated  by  first  defining  both  a  family  of  plants  to  be  controlled 
and  an  "acceptable**  region  of  performance.  Each  of  tbe  plants  should  represent  an  extreme  of 
the  anticipated  parameter  variation.  In  the  case  of  the  nonlinear  plant  under  consideration^  the 
parameter  variation  of  the  linearized  moddb  depends  directly  on  the  desired  operating  range  and 
the  parameter  a.  An  initial  QFT  design  was  developed  for  the  refpon  between  F  =  0  and  Y  m  0.5. 
Five  plants  were  used  to  define  this  operating  range:  V  ss  0,1, 0.2, 0.3, 0.4,  and  0.5,  with  a  assumed 
fixed.  All  of  these  plants  are  of  the  form  r(s)  =  js^^- 


The  region  of  acceptable  performance  is  ddlned  by  upper  and  lower  transfer  function  spec¬ 
ifications,  T/t(/(j)  and  r;u,(^)*  The  magnitude  verses  frequency  plots  of  these  transfer  fnnetions 
exhibit  the  range  of  acceptable  closed-loop  responses.  For  this  case,  Tjtv(s)  is  an  underdamped 
response  with  Afp  =  1.12,  and  f,  =  1.62jecs.  T/tc(a)  is  an  overdamped  function  with  t,  a  1.62secs. 


They  are  given  as: 


Ttur(s)  = 


.8I16(s+H.814) 
(s  +  1.96±  j2.86) 


(243) 


Pigaxc  2.11  QFT  coiopcnsatioa 

and, 

39.99 

^  (3  f  3.45)(«+2.7±i2.074)' 

The  resulting  ’’thumbprint"  of  acceptable  responses  is  determined  by  the  frequency  reeponaes 
of  (2.33)  and  (2.31).  This  is  shown  graphically  in  Fignre  2.12.  As  a  general  rule,  the  larger  the 
thumbprint  of  acceptable  responses,  the  simpler  the  design  will  be.  The  thumbprint  can  also  be 
made  so  small,  or  the  range  of  plant  variation  so  large,  that  successfal  design  is  impossible.  The 
additional  poles  and  zeros  in  (2.33)  and  (2.34)  are  to  ensure  that  the  upper  and  lower  thumbprint 
boundaries  ate  always  diver^ng  for  lower  freqnencies,  as  shown  in  Fignre  2.12.  For  any  given 
frequency,  the  maximum  allowable  change  in  the  frequency  response  due  to  parameter  variation 
is  the  difTereuce  between  the  magnitude  of  the  upper  and  lower  bounds.  This  ’alue  is  dedgnated 
^«(M)- 
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Figtir«2.I2  The  frequency  rexpunsee  of  TjRy(e)  and  r/u.(3),  comprising  the  QFT  response 
thumbprint 

The  development  of  the  compensator  G(s)  is  accomplished  on  a  Nichols  Chart  similar  to  that 
shown  in  Figure  2.13  for  the  current  design  example.  The  open-loop  magnitude  and  phase  data 
for  the  least  stable  (closest  to  the  -180^/0  dB  point)  plant  is  plotted  on  the  Nichols  Chart.  In 
this  case,  the  response  for  the  T  0. 1  plant  is  plotted.  The  other  lines  on  the  chart  represent 
boundaxiea  for  certain  frequencies  of  interest.  These  boundaries  correspond  to  the  points  on  the 
Nichols  Chart,  for  any  value  of  phase  shift,  at  which  the  maximum  uncertainty  in  (caused  by 
parameter  variation)  is  equal  to  the  maximum  variation  allowed  by  If  the  plot  of  the  least 

stable  plant  falls  below  or  to  the  left  of  these  performance  bounds  for  the  frequency  of  interest, 
then  the  thumbprint  of  responses  cannot  be  obtained. 

Boundary  lines  can  be  determined  for  both  tracking  response  (input  at  the  reference  point) 
and  for  disturbance  rejection  (input  directly  into  the  plant).  The  boundary  which  imposes  the 
most  severe  restriction  (highest  in  magnitude  on  the  Nichols  chart)  for  a  particular  frequency  of 
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interest  is  plotted  on  the  Nichols  Chart.  It  can  be  shown  that  QFT  control  of  a  nooliuear  plant  will 
always  involve  rejection  of  am  unwanted  input  at  the  plant  [21].  Therefore,  a  good  signal  rejection 
capability  should  be  designed  into  the  compensator  The  boundaries  plotted  in  Figure  2.13  ate 
almost  all  due  to  signal  rejection  requirements.  The  cylindrical  boundary  around  the  - 180*/0  dB 
point  represents  the  desired  stability  margins  of  the  system.  The  Nichols  plot  of  the  compensated 
system  should  pass  to  the  right  of  or  below  this  encircling  boundary. 


Clearly,  the  frequency  response  in  Figure  2.13  does  not  meet  the  established  boundaries. 
The  plot  of  the  least  stable  plant  falls  considerably  below  the  composite  bounds  and  violates  the 
stability  margin  boundary.  Therefore,  a  cascade  compensator  is  required.  With  the  addition  of  a 
compensator  Gfs)  the  response  shown  in  Figure  2.14  is  obtained.  Figure  2.14  shows  that  all  of  the 
tracking  and  disturbance  rejection  bounds  are  being  met.  The  compensator  used  to  obtain  this 
response  is  given  by: 


477000(s  +  0.1)(s  +  Q.3)(s  4-  7.92) 
s(3  + l.44)(s  + 51.13  ±74.5S1)  ' 


(2.35) 


Successful  compensator  design  ensures  that  the  mriaiion  in  frequency  response  dne  to  parameter 
uncertainly  is  less  than  or  equal  to  Snijui).  The  next  step  is  to  ensure  that  this  variation  ocenn 
within  the  absolute  limits  set  by  TruC#)  auJ  in  Figure  2.12.  Figure  2.15  shows  the  the 

frequency  response  of  a  successful  design.  The  outer  lines  are  the  boundaries  established  by  the 
upper  and  lower  tracking  bounds.  The  inner  hues  represent  the  frequency  response  variatioa  of 
G{s)P(s)  in  closed'loop.  The  jP(s)  required  to  obtain  this  response  is  given  by: 


_  13.053(S+.41) 

(s +  . 425)(s+ 3.54  )(s  + 3.56)’ 

Notice  that  the  6teady-.state  gain  of  /'(s)  is  1. 


(2.36) 


The  performance  of  the  ftnal  closed-loop  system  is  shown  in  Figures  2.16  and  2.17.  Fbr 
aarh  Individual  plant,  the  resulting  tracking  and  disturbance  rejection  responses  are  within  deingn 
specifications.  This  illustrates  the  power  of  the  QFT  synthesis  technique. 
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Unfortunately,  this  G(5)(2-3S)  and  F(s)(2.36)  developed  through  QFT  synthesis  do  not  pro¬ 
duce  a  stable  response  for  the  nonlinear  plant  from  which  the  linearizations  were  derived.  The 
response  quickly  shoots  to  F  =  -l-oo.  This  illuitrates  the  challenge  associated  with  using  linear 
design  techniques  to  control  a  time-varying  linear  plant  with  parameters  changing  faster  than  the 
plant  dynamics  (eigenvalues).  A  redesign  of  the  QFT  controller  is  reqiiired. 

This  redesign  is  accomplished  simply  by  allowing  the  a  parameter  to  vary  between  0.5  and 
1.5  over  the  intended  region  of  operation  from  0  to  0.5.  The  effect  of  this  parameter  variation 
is  to  add  linear  plant  cases  that  are  farther  away  from  the  origin  (farther  from  the  singularity  at 
Y  =  0.0).  Fbr  this  alternate  design,  the  plants  chosen  correspond  to  F  =  0.2, 0.5, 0.7, 1.0,  and  2.5. 
The  compensator  based  on  this  new  extended  family  of  plants  is: 


1228(J  -I-  Q.l)(s  +  0.3)(s  +  7.92) 
s(s-H.44)(s-|-210) 


(2,37) 


This  G(8)  yields  the  freqnency  response  shown  in  Figure  2.18.  The  prefilter  given  by  (2.36) 
is  sufficient  for  this  design  as  well,  so  no  F(s)  synthesis  is  necessary.  Figures  2.19  and  2.20  show 
the  responses  of  each  linearized  plant  to  reference  and  disturbance  unit  step  Inputs.  Note  that  the 
linear  performance  for  both  QFT  designs  is  very  similar.  Indeed,  they  both  even  share  the  same 
precompensator. 

The  new  design  is  also  effective  when  applied  to  control  of  the  nonlinear  plant.  Figure  2.21 
shows  the  response  of  the  nonlinear  plant  with  QFT-based  compensation  to  a  0.5  step  input  at 
the  reference,  for  various  values  of  a.  Notice  that  the  o  variation  has  almost  no  perceptible  effect 
on  the  response  of  the  system.  The  plots  fall  almost  exactly  on  top  of  one  another,  showing  the 
insensitivity  of  the  QFT  design  to  parameter  uncertainty. 

All  plots  used  in  this  section  were  obtained  using  the  QFT  CAD  package  developed  by  Oded 
Yaniv  at  the  University  of  Tel  Aviv,  Israel.  For  additional  information  on  QFT  the  reader  is  referred 
to  reffirences  [4]  and  [29]. 
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Figure  2.14  Nichols  Chart  of  worst-case  plant  in  series  with  G(s) 
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Figure  2.15  Frequency  response  plot  showing  the  variation  of  the  family  of  plants  is  inside  the 
maximum  variation  allowed  by  Tgoi*)  Tiu.(») 
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Figure  2.17  Closed-loop  time  responses  of  all  live  linear  plants  to  unit  stq>  at  tlte  plant  input 
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Figure  2.18  Nichols  chart  showing  the  open-loop  frequency  response  of  G(s)P(s)  based  cm  a  mod¬ 
ified  family  of  plants 
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Figure  2.19  Cloeed-loop  time  responses  of  all  live  alternate  linear  plants  to  unit  step  at  the  refer¬ 
ence  input 


Figure  2.21  Response  of  nonlinear  plant  to  unit  step  reference  input  for  a=n.5,l*0,  and  1.5 
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‘ntble  2.1  Univene  of  Discoone  for  FLC 

2.5.S  FStzzy  Logic  Control.  A  Fiuzy  Logic  Controller  can  be  developed  using  tbe  perfor¬ 
mance  of  the  Dynamic  Inversion  controller  as  a  basis  for  FWzy  tnning.  The  first  step  in  constructing 
this  controller  is  to  determine  the  inputs  required  to  provide  adequate  information  on  the  state  of 
the  plant.  In  this  case,  the  linearized  plant  is  second  wder,  so  it  is  reasonable  to  assume  that  E 
and  E  will  be  adequate  to  control  the  nonlinear  plant.  This  assumption  proves  correct. 

E  and  E  will  be  input  into  the  Fizzy  controUer  and  then  fuzzified  as  discussed  in  Chapter  1. 
AU  Fuzzy  membership  functions  in  this  controll^  are  of  the  form: 

=  (2.38) 

where  2  represents  the  mean  of  the  function  and  is  a  variance-lihe  term  expressing  the  spread 
of  the  function.  The  equation  (2.38)  produces  membership  functions  similar  to  those  shown  in 
Figure  1.1.  The  Gaussian  form  for  the  membership  functions  simplifies  implementation  of  the  FLC 
on  SIMULINK  and  allows  for  simple  tuning  of  the  Fizzy  sets  by  altering  i  and  o’.  Note  that  the 
membership  function  is  completely  specified  by  these  two  quantities. 

The  universe  of  discourse  for  E  and  E  suggested  by  the  DI  simulation  is  shown  in  Table  2.1. 
This  data,  however,  assumes  that  the  plants  remain  on  a  ideal  trajectory  between  0  and  1.  Because 
this  may  not  always  be  the  case,  the  universe  of  discourse  for  the  FLC  will  be  expanded  slightly. 
Note  also  that  by  choosing  this  range  for  E  and  F,  the  controller  will  only  be  effective  for  Jtr»/(t)  < 
1.  The  FLC  will  not  necessarily  break  down  for  Jtr»/(t)  >  1,  but  the  results  will  be  unpredictable. 

Next,  the  granularity  of  each  input  must  be  determined.  Granularity  is  simply  the  number 
of  Fizzy  Sets  which  will  be  used  to  quantify  the  universes  of  discourse  of  the  input  variables. 
The  objective  is  to  use  the  fewest  sets  posable  to  adequately  represent  the  data  contained  in  the 
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l^ble  2.2  Membership  Factions  for  FLC 


input.  In  this  case,  that  number  is  3  for  each  variable:  ’’positive”, "zero”,  and  "negative.”  So  once 
fuzzified,  each  input  will  be  considered  to  be  some  combination  of  positive,  zero,  and  nq;aitive.  As 
mentioned  above,  each  membership  function  for  the  linguistic  variables  will  be  completdy  defined  by 
an  appropriate  x  and  a^.  These  values,  chosen  using  "en^neering  insight,”  are  shown  in  Table  2.2. 


By  modifying  the  DI  simulation,  a  time  history  of  and  V.  can  be  generated  for  "suc¬ 
cessful”  plant  control.  This  data  can  be  used  to  generate  a  first-cut  at  the  IHizzy  inference  system 
which  links  the  compensator  inputs  E  and  E  to  the  output  V.  The  following  rules  can  be  generated 
by  inspection  of  the  DIC  simulation  data: 

1.  IF  i;  is  P  AND  ^  is  N  THEN  U  is  PS. 

2.  IF  £;  is  P  AND  E  is  Z  THEN  U  is  PL. 

3.  IF  £  is  P  AND  E  is  P  THEN  U  is  PL. 

4.  IF  E  is  Z  AND  E  is  N  THEN  U  is  NS. 

5.  IF  E  is  Z  AND  E  is  Z  THEN  U  is  O. 

6.  IF  E  is  Z  AND  E  is  P  THEN  U  is  PS. 

7.  IF  E  is  N  AND  E  is  N  THEN  U  is  NL. 

8.  IF  E  is  N  AND  E  is  Z  THEN  U  is  NS. 

9.  IF  E  is  N  AND  E  is  Z  THEN  U  is  PS. 

In  the  above  rules,  P  is  "Positive,”  Z  is  "Zero,”  N  is  "Negative,”  PS  is  "Positive  Small,”  PL 
is  "Positive  Large,”  PS  is  "Positive  Small,”  NS  is  "Negative  Small,”  NL  is  "N^ative  Large,”  and 
O  is  One.  This  Fuzzy  inference  system  addresses  all  possible  combinations  of  the  two  compensator 
inputs  E  and  E. 
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T^ble  2.3  Values  for  U  suggested  by  Dynamic  Inversion  simulation 


In  many  FLC  applications,  the  next  step  would  be  to  define  the  universe  of  discourse  and 
membership  functions  for  the  output  U.  In  this  case,  however,  a  value  of  U  from  the  01  simulation 
can  be  associated  with  each  combination  of  input  arguments.  The  values  of  U  can  still  be  viewed 
as  Fuzzy  Sets,  but  in  the  limiting  case  where  each  premise  maps  to  a  ddta  function  corresponding 
to  a  single  point  in  the  output  universe  of  discourse.  The  values  for  U  are  shown  in  Table  2.3. 

The  Fuzzy  Inference  system  above  is  the  heart  of  the  FLC,  linking  the  fuzzified  contndln 
inputs  to  the  appropriate  outputs.  What  remains  to  be  resolved  is  the  exact  numerical  procedure 
that  will  be  used  to  implement  these  rules  and  how  multiple  control  conclusions  will  be  resrfived  to 
produce  a  single  ’’crisp’*  output.  This  is  accomplished  in  several  steps.  First,  the  AND  statement 
in  the  argument  of  each  rule  must  be  resdved.  The  standard  solution  to  the  Fuzzy  AND  statement, 
proposed  by  Zadeh  and  others  is  [28]: 

C{x)  =  A{x)AN DB{x)  -*■  fteix)  =  Mir,{ftji{x)fHB{*))-  (2.39) 

Therefore,  whichever  term  in  the  premise  is  activated  the  least  will  determine  the  level  of  activation 
of  the  implicated  control  action.  In  this  way,  every  possible  control  action  will  receive  a  membership 
value  which  indicates  the  degree  of  confidence  the  controller  has  that  each  course  of  action  is  the 
"correct”  one. 

Before  being  output  to  the  nonlinear  plant,  U  must  be  defuzzified.  To  accomplish  this,  the 
defuzzification  routine  of  [13]  will  be  used: 

tf(<)  =  (2.40) 

LA=ilk/i 
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TVible  2.4  Values  for  U  implicants  based  on  tuning  through  simulation 

The  numerator  of  this  expression  is  simply  the  sum  of  all  possible  courses  of  action  times  the  d^ree 
of  confidence  the  controller  has  that  the  ^ven  course  of  action  is  correct.  The  denominator  is  a 
normalization  term  to  ensure  that  the  sum  of  all  the  membership  functions  in  1.  Note  also  that 
because  the  membership  functions  are  Gaussian,  every  resultant  control  action  will  be  activated 
to  some  degree  at  all  times.  Then  the  state  of  the  actual  plant  is  far  away  from  the  mean  of  the 
membership  function,  however,  the  activation  will  be  very  small. 

This  FLC  was  implemented  and  tuned  using  SIMULINK.  Not  surprisingly,  the  values  for 
U  suggested  by  the  DI  simulation  did  not  produce  an  acceptable  response.  The  compensator  is 
unable  to  produce  the  required  dynamic  range  of  U  outputs.  Many  possible  changes  could  have 
been  undertaken  to  correct  the  system  response.  However,  rather  than  change  the  structure  of  the 
Fuzzy  Sets  or  resort  to  a  finer  granularity,  it  was  simpler  to  change  the  values  of  the  implicate  U 
values.  The  new  values  are  given  in  Table  2.4. 

The  simulation  set-up  is  shown  in  Figure  2.22,  with  the  FLC  shown  as  a  single  block.  Fuzzy 
Controller.  The  FLC  as  implemented  in  SIMULINK  is  shown  in  Figure  2.23.  F\izzification  of  inputs, 
ANDing  of  the  premises  and  defuzzification  of  resultants  is  all  carried  out  in  this  simulation.  The 
closed-loop  responses  of  the  system  to  step  inputs  of  Ref  =  lu— 1(f)  and  1.5tt-l(f)  are  shown  in 
Figure  2.24. 

The  purpose  of  showing  the  response  of  the  FLC  simulation  to  the  1.5  unit  step  is  to  show  the 
effect  of  reference  inputs  beyond  the  range  of  the  universe  of  discourse  for  which  the  compensator 
was  designed.  In  this  case  the  response  does  not  exhibit  a  0  steady-state  error.  As  is  to  be  expected, 
the  performance  degrades  further  for  larger  and  larger  step  inputs. 
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Figure  2.22  Closed-loop  control  of  nonlinear  plant  with  Fuzzy  Logic  Control 


Consider,  however,  variations  in  the  plant  parameter  >1,  indicating  mismodeling  or  non- 
statistical  uncertainty.  Figure  2.25  shows  the  system  response  to  a  unit  step  given  three  different 
values  for  A.  The  response  corresponding  to  mismodeling  A(A  =  0.5  and  A  =  1.5)  are  much  closer 
to  the  correct  response  than  are  the  responses  for  the  DI  compensator.  This  demonstrates  that 
FLC  offers  advantages  in  robustness  over  other  approaches. 

As  a  final  example,  Figure  2.26  shows  the  response  of  both  compensators  to  a  step  input 
in  the  face  of  a  time-varying  A  parameter.  In  this  case  a  =  0.5  -f-  .2t.  The  FLC  is  more  able  to 
maintain  the  shape  of  the  desired  response  in  the  face  of  this  type  of  plant  variation. 
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Figure  2.23  Fuzzy  Logic  Compensator  for  Nonlinear  Plant 
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Figure  2.25  Time  response  of  FLC/nonlinear  plant  to  step  input  assuming  modeling  errors 
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Figure  2.26  Time  response  of  FLC  and  DI  controllers  to  time- varying  a 


t.S  Summary 

In  this  duster,  n  nonlinear  plant  was  developed  to  explore  hybrid  linear/FNusy  owtrol  a>n- 
cepts.  Desired  performance  chairacteristics  were  established,  and  a  closed-loop  modd  was  presented. 
The  desired  characteristics  are  such  that  compensation  is  required  to  properly  control  tlm  idant. 

In  order  to  establish  a  control  baseline  and  to  develop  insights  into  the  nonlinear  contnd 
problem,  three  conventional  designs  were  presented.  A  purely  nonlinear  technique  was  donon- 
strated  using  dynamic  inversion.  This  will  serve  as  the  paradigm  for  a  nonlinear/FVizzy  controller 
in  Chapter  4.  Second,  a  conventional  F\izzy  Logic  Controller  was  developed  to  show  the  structure 
and  demonstrate  the  insendtivity  to  modeling  errors.  Finally,  a  linear  design  was  developed  based 
on  Quantitative  Feedback  Theory.  This  serves  as  a  benchmark  for  performance  of  other  linear 
designs. 
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III.  Model-Based  Fuzzy  Logic  Control  Linear  Design  Considerations 
3.1  Introduction 

The  previous  chapter  demonstrated  the  potential  of  FLC  controllers  to  adequatdy  control  a 
nonlinear  plant  for  a  single  prespecified  unit  step  input.  The  resulting  FLC  actuaUy  had  advantages 
in  robustness  over  more  conventional  designs.  The  FLC  was  not  based  on  a  plant  modd,  and 
the  relationship  between  inputs  and  the  desired  outputs  was  determined  in  an  ad  hoc  manner 
through  simulation.  The  remainder  of  this  report  addresses  the  use  of  linear  compensates  to 
control  nonlinear  plants  with  Rizzy  Logic  serving  as  an  "interpreter.”  Because  there  are  quantifiabk 
relationships  between  the  linear  compensators  and  the  nonlinear  plant,  the  extent  of  the  nei- 
statistical  uncertmnty  is  reduced  to  some  extent  compared  to  the  Fuzzy  Lofpc  Contrdler,  and 
enhanced  performance  is  to  be  expected. 

One  possible  means  of  applying  linear  systems  theory  to  nonlinear  control  is  through  the  use 
of  a  bank  of  linear  compensators.  As  discussed  in  Chapter  1,  a  nonlinear  plant  can  be  linearized 
about  a  given  trim  condition  with  system  dynamics  described  in  terms  perturbations  from  the 
equilibrium  point.  Nonlinear  plants  such  as  aircraft  have  many  such  equilibria,  corresponding  to 
straught  and  level  flight  at  a  given  airspeed  and  altitude.  As  the  nonlinear  plant  changes  equilibria, 
the  linear  perturbation  model  can  be  expressed  as  a  time-varying  linear  function.  This  function 
can  serve  as  the  basis  for  multiple  linear  compensators  which  can  then  be  used  to  compensate  the 
original  nonlinear  plant. 

In  the  Banked  Compensation  MBFLC  approach,  the  compensators  themsdves  would  all  be 
present  in  the  forward  control  path  and  their  outputs  would  be  weighted  according  to  an  estimation 
of  which  compensator  most  accurately  reflects  the  current  state  of  the  system.  When  none  of  the 
models  are  considered  valid,  combinations  of  outputs  or  Fuzzy  Lo^c-induced  caution  would  drive 
the  system  to  a  region  where  the  controllers  are  ^ain  effective.  Figure  1.5  ^ves  a  conceptual  view 
of  a  Banked  Model- Based  Fuzzy  Logic  Controller. 

This  chapter  addr^ses  the  linear  considerations  associated  with  the  proposed  MBFLC  ar¬ 
chitecture.  Because  it  is  not  clear  if  the  compensator  bank  will,  in  fact,  achieve  the  desired  goal. 
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several  areas  will  be  explored  in  this  chapter  to  resolve  this  question.  The  following  issues  will  be 
explored: 

1.  Linear  Cmnpensation  of  Linearized  Plant 

2.  Use  ot  Multiple  Compensators  in  Linear  Systems 

3.  Effect  of  Compensator  Dissimilarities  on  Cloeed-Lo<q>  Linear  System  Response 

4.  Use  of  Multiple  Compensators  in  Time- Varying  Linear  Systems 

Based  on  these  deugn  considerations,  a  prototype  Model-Based  Fuzzy  Lope  Controller  will 
be  proposed  in  the  next  chapter. 


S.2  Linear  Compensation  of  Linearized  PUmt 

At  any  given  equilibrium  condition,  the  dynamics  of  the  plant  about  that  equilibrium  can 
be  described  through  a  linear  perturbation  modd.  A  linear  compensator  based  upon  this  modd 
can  be  used  to  contitd  the  plant,  provided  the  plant  remains  sufficiently  dose  to  the  equilibrium. 
A  succession  of  linear  compensators  based  on  this  concept  wiU  be  used  to  produce  the  Banked 
Contndler. 


As  stated  in  the  previous  ch^ter,  the  form  dT  the  linearized  plant  from  the  last  chapter,  given 


a  fixed  equilibrium  point  is: 


This  is  the  transfer  function  form  of  the  state-space  linear  modd  of  (2.21)  and  (2.22)  shown  in 
Figure  2.2.  The  plant  consists  of  two  poles,  both  real  for  small  values  eS  r,  and  as  a  comjdex  pair 
for  large  values  of  r.  Figure  3.1  shows  the  locus  of  poles  as  r  increases  from  0.  The  break-away 
point  on  the  real  axis  occurs  at  x  =  -0.5.  The  compensation  approach  taken  is  to  caned  the  poles 
of  the  linearized  plant  with  zeros  in  the  compensator.  The  demoninator  of  the  compensator  will 
contain  one  pole  at  the  origin  and  one  pole  at  z  =  —5.  As  the  compensator  gain  is  increased,  the 
poles  at  0  and  -5  will  come  together  at  x=-2.5  and  break  away  at  90*.  To  obtain  pdes  similar  to 
the  modd  (2.26),  a  gain  of  K  =  ^  is  required.  Note  that  this  approach  breaks  down  as  r  0,  so 
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Figure  3.1  Locus  of  pdes  of  the  linearized  plant  as  a  function  of  r 

a  perturbation  in  K  will  be  required  aaY  -*0.  This  will  yield  the  desired  pair  of  dominant  pdes, 
as  given  in  the  model(2.26). 

For  a  ^ven  value  of  r,  the  form  of  the  compensator  is: 


GW  = 


s*  +  5s 


(3.2) 


This  compensator  is  shown  as  implemented  in  SIMULINK  in  Figure  3.2.  This  model  is 
represented  in  subsequent  simulations  as  a  single  block,  TVComp.  The  compensator  was  converted 
to  state  space  form  to  accommodate  an  externally  defined  r.  As  expected,  for  any  constant  value 
of  r  this  compensator  will  yield  the  correct  closed-loop  response. 

This  compensator,  however,  does  not  meet  the  specifications  for  disturbance  rejection  (|C«,i«rt«»e((OI  ^ 
0.1),  so  a  second  design  must  be  considered.  In  order  to  add  robustness  in  the  face  of  noise  injected 
at  the  input  to  the  linearized  plant,  three  poles  and  three  zeros  will  be  added  to  the  compensator. 
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Figure  3.2  Model  of  Compeneator  G(tf)  with  r  externally  defined 


The  form  of  the  compensator  with  disturbance  rejection  is; 

^  ^is^  +  s  +  T){s  +  45f 
^  ’  s(s  +  4.5)(s  +  200)* 


(3.3) 


Notice  that  in  order  to  achieve  disturbance  rejection  yet  maintain  acceptable  tracking  performance 
the  other  poles  in  the  denominator  also  had  to  be  moved.  While  the  performance  of  the  compensator 
^ven  by  (3.2)  produces  a  response  precisely  as  given  in  the  model,  the  response  of  (3.3)  will  be 
slightly  faster.  The  implementation  is  also  considerably  more  complex.  The  SIMULINK  model  of 
(3.3)  is  shown  in  Figure  3.3. 

Because  the  benefit  of  linear  robustness  to  the  control  of  nonlinear  plants  is  not  clear,  two 
simulations  were  run  specifically  to  explore  this  issue.  Both  simulations  are  shown  in  Figure  3.4. 
The  first  tests  the  ability  of  each  compensator  (^ven  by  (3.2)  and  (3.3))  to  control  the  linearized 
plant  when  r  was  allowed  to  vary  with  the  output  (r  =  31^’).  Because  r  is  a  function  of  Y,  this 


3-4 


Figure  3.3  Model  of  Compensator  G(«)  with  r  externally  defined  which  meets  disturbance  rejec¬ 
tion  specifications 

plant  will  be  nonlinear  and  exhibit  many  of  the  same  characteristics  as  the  full  nonlinear  plant  [21]. 
The  second  simulation  tests  the  ability  of  each  compensator  to  control  the  nonlinear  plant  itsdf  for 
small  steps  away  from  equilibrium. 

As  mentioned  above,  r  is  undefined  at  =  0  so  a  small  constant  is  added  to  the  "Thu”  block 
in  the  Figure.  Though  this  small  perturbation  should  not  effect  the  system  response,  all  umulations 
in  Figure  3.4  start  from  a  nominal  condition  =  1  to  be  sure. 

Figure  3.5  shows  that  there  is  very  little  difference  in  tracking  performance  between  the  robust 
controller  and  the  non-robust  controller  for  the  time- varying  linear  plant.  This  Figure  shows  the 
responses  of  the  simulation  for  a  step  input  of  Ref  =  0.5u_i(t).  Both  compensators  have  similu^ 
responses,  so  the  test  is  inconclusive.  Figure  3.6  shows  that  the  robust  controller  exhibits  superior 
signal  rejection  capabilities,  so  it  is  tempting  to  conclude  that  the  robust  controller  is  preferable. 
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Figure  3.4  TOP:  Closed-loop  simulation  given  linear  plant  with  time-varying  plant  parameter  r. 

BOTTOM:  Closed-loop  simulation  given  full  nonlineekr  plant.  Both  simulations  are 
run  with  each  compensator. 

Consider,  however,  the  response  of  the  full  nonlinear  plant,  shown  in  Figure  3.7.  This  Figure 
shows  the  response  of  the  closed-loop  system  to  a  step  input  of  Ref  =  0.1tt_x(t),  and  it  is  dear 
that  the  non-robnst  controller  yields  significantly  better  results  than  the  robust  controller.  The 
difference  in  performance  is  rooted  in  the  fact  that  the  robust  compensator  has  a  significantly 
higher  gain  than  the  nonrobust  controller  and  thus  a  smaller  region  of  attraction.  As  the  step 
input  becomes  smaller,  the  dosed-loop  responses  of  the  compensators  converge.  This  is  clearly  a 
nonlinear  phenomenon. 
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Based  oa  this  analysis,  most  simulations  to  follow  will  be  conducted  using  the  nonrobust  com¬ 
pensator  design.  Which  compensator  is  actually  superior  depends  on  the  the  desired  application. 
The  robust  controller  is  superior  for  disturbance  rejection.  The  nonrobust  controller  contributes 
much  less  overshoot  for  tracking.  Use  of  the  robust  compensator  is  revisited  in  Chapter  5. 
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Figure  3.6  Closed-loop  disturbance  rejection  for  nonlinear  plant  under  both  robust  and  non-robust 
compensation. 
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Figure  3.7  Closed-loop  system  responses  for  nonlinear  plant  under  both  robust  and  non-robust 
compensation. 
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S.S  Use  of  Multipk  Compensators  m  Linear  Systems 

The  next  question  to  be  addressed  is  how  exactly  to  apply  multiple  linear  compensators  to  a 
single  linear  plant.  Of  special  interest  is  the  compensation  of  a  plant  which  is  not  at  a  quiescent 
condition  at  the  starting  time  t,.  This  situation  would  arise  repeatedly  as  each  new  compensator  in 
the  bank  were  activated  due  to  a  change  in  the  equilibrium  point  of  the  plant.  If  the  act  of  changing 
compensators  disturbs  the  nonlinear  plant,  then  the  banked  configuration  may  not  be  desirable. 

A  proper  treatment  of  linear  control  systems  operating  firom  non-zero  initial  conditions  was 
not  found  in  the  literature  search  for  this  thesis.  Therefore,  several  simulations  were  necessary  to 
determine  the  appropriate  structure  for  the  MBFLC.  Because  both  the  compensator  bank  and  the 
plant  have  the  ability  to  store  energy,  both  of  these  cases  must  be  considered  separately,  and  in 
terms  of  their  interaction. 

S.S.l  Positive  Internal  Compensator  Energy/No  Internal  Plant  Energy.  Figure  3.8  shows 
a  simple  simulation  used  to  determine  the  effect  of  "chaurging”  the  compensator  before  applying 
it  to  the  quiescent  plant.  In  this  Figure  the  compensator  will  be  receiving  the  error  signal  E  at 
all  times.  A  switch  just  beyond  the  compensator  prevents  the  output  of  the  compensator  from 
reaching  the  linear  plant  for  a  prespecified  length  of  time  t^.  FVom  time  0  to  time  U  the  states  of 
the  compensator  wiU  be  receiving  energy  from  the  error  signal,  while  the  plant  is  forced  to  remain 
at  0.  At  t  =  the  plant  will  be^n  receiving  inputs. 

Figure  3.9  shows  the  effect  on  the  system  response  as  U  increases.  Clearly,  as  U  gets  larger, 
the  response  of  the  plant  system  gets  less  and  less  desirable.  This  shows  that  the  effect  of  energy  in 
the  compensator  is  to  underdamp  the  system  in  dircrt  proportion  to  the  energy  in  the  compensator 
at  fj-l-.  This  demonstrates  that  compensator  charging  is  undesirable,  at  least  when  applied  to  a 
quiescent  plant. 
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Figure  3.8  Simulation  to  determine  the  effect  of  charged  states  in  compensator  applied  to  linear 
system 
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Figure  3.9  Closed-loop  response  of  linear  plant  gjiven  varying  "charging”  times  for  the  linear  plant 


S.S.2  No  Internal  Compensator  Energy/ Positive  Internal  Plant  Energy.  The  previous 
result  assumed  the  plant  was  at  0  initial  conditions  at  the  time  of  compensator  activation.  A 
similar  simulation  can  be  run  but  with  the  plant  itself  starting  from  a  nonzero  initial  condition, 
and  the  compensator  starting  from  0  internal  energy.  There  are  two  possible  situations  which 
could  arise:  1)  A  plant  starting  from  a  nonzero  quiescent  position  and  2)  A  plant  starting  with 
non-qulescent  initial  conditions. 

In  both  cases,  linear  systems  theory  can  be  used  to  predict  the  outcome  [31].  For  the  first 
situation,  the  desired  output  form  will  only  be  obtmned  if  an  additional  constant  input  is  added  to 
the  input  from  the  compensator.  The  value  of  this  input  is  equal  to  the  steady-state  input  which 
would  have  induced  the  existing  quiescent  condition.  The  linearized  plant  plant  in  this  research 
effort,  for  example,  is  in  equilibrium  at  xi  =  =  U.  Therefore,  if  the  plant  states  here  have 

initial  conditions  xi  =  x^  =  0.2,  the  required  control  input  to  the  plant  would  be  0.2  plus  the 
control  determined  by  the  linear  compensator.  The  effect  is  simply  to  move  the  starting  point  of 
the  simulation,  and  this  has  no  consequences  for  LTI  systems. 

The  second  case  is  much  more  difficult  to  predict  in  general.  The  complete  response  of  a 
linear  system  is  always  equal  to  the  response  of  the  system  due  to  initial  conditions  (homogeneous 
solution)  plus  the  response  of  the  system  to  driving  terms  (particular  solution)  [17].  When  the 
system  is  not  at  equilibrium  there  will  be  a  homogeneous  solution  intermingled  with  the  particular 
solution.  All  that  can  be  said  with  certaunty  is  that  the  response  will  deviate  from  the  desired 
trajectory  to  some  degree  proportional  to  the  strength  of  the  homogeneous  response. 

This  is  illustrated  using  the  simulation  shown  in  Figure  3.10.  Here,  two  identical  compen¬ 
sators,  Gi{s)  and  62(5),  are  included  in  the  forward  path,  with  the  input  to  the  plant  bring 
determined  by  a  switch.  The  position  of  the  switch  is  determined  by  the  current  output  of  the 
system.  The  switch  will  toggle  from  the  top  compensator  to  the  bottom  compensator  at  T  =  0.5, 
halfway  to  the  reference  input  of  T  =  1.0.  A  second  switch  prevents  Gj{s)  from  receiving  the 
error  signal  until  Y  =  0.5.  The  effect  is  that  Gjis)  will  "inherit”  the  plant  when  the  plant  is  in  a 
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Figure  3.10  SIMULINK  simulation  to  test  the  effect  of  a  quiescent  compensator  inheriting  a  non- 
quiescent  plant 

non-equilibrium  state.  The  state  of  the  plant,  however,  is  not  arbitrary:  it  is  precisely  the  state 
that  C73(s)  would  have  driven  the  plant  through  had  it  been  utilized  from  to* 

The  result  of  this  simulation  is  shown  in  Figure  3.11.  The  optimal  closed-loop  trajectory  is 
also  shown  in  the  Figure  (optimal  being  the  closed-loop  response  using  a  single  compensator  for  all 
time).  The  compensator  with  zero  initial  conditions  is  unable  to  correctly  drive  the  system  when 
the  plant  has  internal  energy  at  the  time  the  switch  is  thrown.  An  identical  compensator,  then 
cannot  simply  ’’pick  up  where  the  first  compensator  left  off.”  A  mechanism  in  the  contr<d  strat^y 
must  be  devised  to  account  for  the  existing  plant  energy. 
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Figure  3.11  Clo«ed-loop  respouse  of  quiescent  compensntcsr/non-quiescent 
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S.S.S  Positive  Internal  Compensator  Energy/Positive  Internal  Plant  Energy.  To  account 
for  the  energy  internal  to  the  plant,  consider  the  simulation  shown  in  Figure  3.12.  In  this  Figure, 
the  plant  is  linear  and  corresponds  to  the  equilibrium  condition  r  =  3  or  F  =  1.  A  bank  of  two 
compensators-  r  set  to  3  for  both-  feed  into  the  plant,  and  they  are  weighted  at  their  error  signal 
inputs.  In  this  case,  the  only  two  inputs  into  either  compensator  are  possible:  1)  /{  =  0  or  2)  the 
actual  error  signal,  R  =  E.  The  difference  between  this  simulation  and  the  simulation  in  Figure  3.10 
is  that  now  (?i(s)  has  the  ability  to  effect  the  plant  for  all  time. 

The  closed-loop  response  of  this  system  for  various  switching  times  is  shown  in  Figure  3.13. 
Regardless  of  the  time  when  the  compensators  are  switched,  the  closed-loop  response  will  be  unaf¬ 
fected  and  equivalent  to  the  response  of  a  single  compensator  in  series  with  the  plant.  By  driving 
the  error  signal  into  the  <7i(s)  to  0  at  the  time  of  compensator  switch,  but  allowing  the  output  of 
Gi(s)  to  continue  to  drive  the  plant,  the  energy  currently  in  the  plant  is  accounted  for. 

In  this  situation,  the  bank  can  be  thought  of  as  a  single  compensator.  When  Gs(s)  inherits 
the  system,  it  has  no  "knowledge”  of  the  internal  plant  energy.  It  simply  receives  an  error  signal  and 
produces  an  output  "assuming”  the  plant  is  starting  at  a  quiescent  condition.  This  is  far  from  being 
the  case,  and  Gjis)  supplies  a  much  larger  control  input  than  is  actually  required.  This  overzealous 
output  from  Ct3(s)  is  apparent  from  the  response  of  Figure  3.11.  In  this  latest  simulation,  however, 
Gi(s)  actually  corrects  for  this  error  by  introducing  a  negative  control  at  the  summing  junction  to 
the  plant.  This  ensures  that  the  Xi  and  X2  already  present  in  the  system  at  the  time  of  control 
transfer  are  accounted  for  from  the  switch  time  forward. 

It  is  important  to  see  Just  what  Gi{s)  is  doing  to  make  this  simulation  successful.  The  input  to 
Gi(s)  is  forced  to  0  at  t  =  0.25  secs,  but  the  integrators  in  the  compensator  will  continue  to  produce 
an  output,  which  will  eventually  reach  steady-state.  In  this  case,  Gi(tgo)  =  0.5622.  Figure  3.14 
shows  the  input  signal  to  the  plant  (sum  of  control  signals  for  both  compensators)  and  breaks  it 
down  in  terms  of  the  contribution  by  each  compensator  alone.  Either  forcing  the  output  of  Gi(s)  to 
zero  or  replacing  the  compensator  with  a  constant  signal  equal  to  Gi(tco)  have  destabilizing  effects 
on  the  system  output.  Therefore,  the  dynamics  of  Gi(t)  are  important. 
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Figure  3.12  Simulation  of  banked  compensation  technique  for  a  linear  plant 


A  bank  of  linear  compensators  can  be  implemented  in  an  alternate,  but  more  straightforward 
way  by  simply  postweighting  the  compensators,  as  shown  in  Figure  3.15.  In  this  case,  both  com¬ 
pensators  receive  the  error  signal,  though  only  one  influences  the  plant  at  any  time.  One  drives 
the  plant  from  to  while  the  second  from  on.  In  this  case,  the  switching  process  once 
again  has  no  effect  on  the  response,  regardless  of  the  switching  time.  This  result  is  not  surpriung. 
The  compensators  are  identical,  receiving  identical  error  signal  histories,  and  so  they  produce  the 
same  input  signal  u  to  the  plant.  Regardless  of  which  is  actually  driving  the  plant,  the  response 
will  be  the  same.  This  also  serves  to  illustrate  that  only  when  plant  and  compensator  states  are 
adequately  charged  will  the  response  be  unaffected  by  the  switching  operation. 

An  important  conclusion  can  be  drawn  from  this  analysis.  There  must  be  internal  energy  in 
the  compensator  to  account  for  internal  energy  in  the  plant.  Whether  the  internal  compensator 
energy  is  provided  by  the  previous  compensator  (pre-wdghted  case)  or  built  up  in  the  compensator 
about  to  inherit  the  system  (post-weighted  case)  is  a  question  of  implementation.  When  a  linear 
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Figure  3.13  Response  of  Identical  Compensators  in  a  banked  configuration 

compensator  inherits  a  system  that  is  not  at  equilibrium,  undesired  oscillations  will  occur  unless 
an  added  signal  is  included  account  for  the  non-equilibrium  condition.  FVirther,  the  form  of  this 
compensation  signal  is  approximated  by  the  linear  response  of  Gi{s)  in  Figure  3.14  for  t  >  0.5  secs. 
This  general  conclusion  will  be  the  basis  for  Model-Based  Fuzzy  Logic  Controller  design  in  Chapter 
5. 


3.4  Effect  of  Compensator  Dissimilarities  on  System  Response 

In  the  above  simulations,  the  two  compensators  in  the  bank  were  identical.  By  using  dissimilar 
compensators,  the  question  of  compensator  autonomy  can  be  addressed.  The  ideal  scenario  would 
be  to  have  each  hand-over  point  be  equivalent  to  a  quiescent  condition.  In  this  case,  each  new 
compensator  could  start  "fresh,”  without  regard  to  the  past  history  of  the  system.  Unfortunately, 
the  effects  of  unanticipated  internal  energy  are  once  again  introduced  into  the  system. 
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Figure  3.14  Output  of  compensator  bank,  shown  in  terms  of  each  compensator  contribution 


To  test  compensator  autonomy,  one  of  the  compensators  in  Figure  3.12  was  changed.  A 
second  compensator  was  developed,  based  on  the  same  linear  plant,  but  with  a  settling  time  t,  =  2.5 
seconds  rather  than  1.62  seconds.  The  transfer  function  for  this  new  compensator  is: 


Galt{^)  = 


2.44(s»  -f  s  +  3) 
s{s  +  3) 


(3.4) 


Gi(s)  was  replaced  by  this  new  compensator.  The  results  of  this  simulation  for  various  switching 
times  are  shown  in  Figure  3.16.  The  responses  at  the  far  right  and  far  left  correspond  to  =  oo 
secs  (only  G^LTi^)  drives  the  plant)  and  =  0  (only  Ggls))  respectively.  In  the  intermediate 
simulations  the  effect  of  the  two  different  settling  times  is  appairent,  but  the  correct  waveform  is 
preserved.  Notice,  however,  that  none  of  the  responses  influenced  by  both  compensators  have  the 
correct  overshoot,  even  though  both  are  designed  to  the  same  (  specifications.  This  is  the  effect  of 
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Figure  3.15  Simulation  applying  a  bank  of  identical  linear  compensators  to  a  linear  plant 

the  second  compensator  inheriting  a  plant  which  is  in  an  unanticipated  state:  the  states  are  not 
along  the  trajectory  which  the  compensator  would  have  induced  had  it  been  in  control  since  (  =  0. 

This  effect  is  inverted  when  GjiL‘r{s)  is  substituted  for  In  this  case  the  response  starts 

out  much  faster  due  to  (7i(s),  then  slows  down  due  to  This  leads  to  an  overshoot  smaller 

than  is  desired  when  both  compensators  are  involved  in  controlling  the  plant. 

If  the  control  transfer  point  behaved  like  an  equilibrium  point,  then  the  overshoot  could  be 
calculated  by  Af  =  YAtTrantjw  +  1-12  *  (1  —  kOitTran«/(r)>  Because  the  actual  overshoot  does  not 
follow  this  relation,  the  conclusion  is  that  ^he  switch  point  must  not  imply  an  equilibrium  point. 
Therefore,  error  in  overshoot  must  be  due  to  a  mismatch  between  the  internal  energy  in  the  plant 
and  the  internal  energy  of  the  compensator.  Since  (7i(s)  is  the  only  compensator  element  with 
internal  energy  at  the  switching  time,  it  must  be  unable  to  account  for  the  energy  it  induced  in 
the  plant  when  <?3(s)  attempts  to  drive  the  plant  along  a  different  state  trajectory.  The  output  of 
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Figure  3.16  Simulation  of  banked  compensation  technique  for  unlike  compensators 

Gi{s)  for  all  time  after  the  error  signal  is  driven  to  zero  must,  therefore,  be  adequate  only  when  the 
compensator  continues  along  the  same  trajectory  G’i(s)  would  have  induced  had  it  been  in  control 
for  all  time. 

A  similar  analysis  was  carried  out  using  the  post- weighting  approach  of  Figure  3.15.  Inserting 
mismatched  compensators  into  this  configuration,  led  to  large,  lightly-damped  oscillations  in  the 
closed-loop  response.  This  suggests  that  the  final  Model-Based  Fuzzy  Logic  Controller  should 
incorporate  pre-weighted  compensators  to  minimize  the  effect  of  modeling  errors. 
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S.5  Use  of  Multiple  Compensators  in  Time-Varying  Linear  Systems 

As  was  mentioned  in  the  introduction,  the  nonlinear  plant  can  be  approximated  by  a  time- 
varying  linear  model  by  relaxing  the  fixed  equilibrium  assumption  of  the  linearization  process.  The 
resulting  transfer  function  is  time-varying  and,  in  fact,  nonlinear  as  the  value  of  r.  However,  by 
conveniently  neglecting  the  fact  that  r  is  a  function  of  Y,  the  plant  can  (and  will)  be  considered 
linear,  time-varying.  The  plant  will  be  referred  to  in  this  report  as  the  linear  time-varying  linear 
plant  (LTV).  Simulation  shows  that  a  compensator  which  adequately  controls  the  LTV  plant  will 
also  control  the  nonlinear  plant  to  the  extent  that  the  small  perturbation  assumption  is  not  violated. 

Figure  3.17  shows  a  simulation  to  test  the  effectiveness  of  a  single  linearized  compensator 
against  the  LTV  plant.  The  plant  is  operating  from  a  nominal  value  of  C{t)  =  1.  Figure  3.18 
shows  the  response  of  the  system  for  step  inputs  of  various  magnitudes.  Notice  that  the  response 
steadily  decays  as  the  system  moves  away  from  the  nominal  value.  The  largest  step  input  for 
which  the  compensator  will  produce  an  acceptable  closed-loop  output,  will  be  referred  to  as  the 
re^on  of  attraction  for  that  compensator.  That  is,  as  the  long  the  plant  is  operating  within  a 
compensator’s  region  of  attraction,  that  compensator  will  be  capable  of  inducing  an  acceptable 
closed-loop  response. 

Based  solely  on  the  use  of  compensators  within  their  regions  of  attraction,  a  non-fuzzy  Banked 
Compensator  for  the  LTV  plant  can  be  postulated,  as  shown  in  Figure  3.19.  Unlike  the  Banked 
Compensator  for  the  linear  system,  shown  in  Figure  3.12,  Figure  3.19  is  designed  to  switch  at  a 
particular  output  value  rather  than  a  particular  time.  In  this  case,  the  objective  of  the  system  is  to 
drive  the  LTV  plant  from  y(t)  =  1  to  Y{t)  =  1.2.  The  bank  consists  of  two  controllers,  one  about 
y(t)  =  1  and  one  for  Y{t)  =  1.2.  The  switch  will  change  from  one  compensator  to  the  other  at 
y(f)  =  1.1,  halfway  between  compensators.  Even  at  y(f)  =  1.1,  both  compensators  demonstrate  a 
reasonable  capability  to  control  the  plant.  Therefore,  these  compensators  have  overlapping  regions 
of  attraction. 

The  closed-loop  response  of  this  system  is  shown  in  Figure  3.20.  The  banked  compensator 
has  had  the  opposite  effect  than  was  intended.  The  response  is  actually  worse  than  the  response 
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Figure  3.17  Simulation  testing  the  validity  of  linear  compensators  on  time- varying  linear  plant 

of  the  single  compensators  alone,  shown  in  Figure  3.  IP  Clearly,  an  analog  to  the  error  source 
identified  last  section  is  also  at  work  here.  Though  the  compensators  are  all  designed  to  the  same 
specifications,  they  are  both  customized  for  a  single  equilibrium  point.  Once  the  LTV  plant  moves 
away  from  the  equilibrium  a  modeling  error  is  introduced  between  the  plant  and  the  compensator. 
This  results  in  unmodeled  control  effects.  Auxiliary  processing  is  required  to  produce  successful 
non-Fuzzy  Banked  Compensation  for  the  LTV  plant. 
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Figure  3.19  Prototype  MBFLC  for  Time- Varying  Linear  Plant 
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S.6  Controiler  Constdemtions  Based  on  Linear  Analysis 

There  are  at  least  three  possible  causes  for  the  failure  of  the  banked  compensator  above,  for: 

1.  E!ach  compensator  is  not  valid  over  its  entire  rejpon  of  operation.  Even  though  the  compen¬ 
sators  provide  reasonable  performance,  the  time  of  hand-over  the  plant  is  off  the  trajectory 
which  the  inheriting  compensator  would  have  induced  had  it  been  driving  an  equivalent  linear 
plant  the  entire  time.  This  introduces  a  non-statistical  error. 

2.  Each  compensator  is  unable  to  properly  account  for  the  dynamics  induced  while  driving  the 
nonlinear  plant.  Because  the  plant  is  changing  with  time,  none  of  the  compensators  will 
be  valid  for  all  time.  The  compensator  will  be  accounting  for  the  dynaunics  it  had  induced, 
assuming  the  plant  is  precisely  the  linear  plant  for  which  the  compensator  is  designed.  This 
is  another  source  of  non-statistical  error. 

3.  The  region  of  attraction  is  a  limitation  that  applies,  not  Just  to  the  equilibrium  condition  of 
the  plant,  but  to  the  magnitude  of  the  error  signal  into  the  compensator  as  well.  The  error 
signal  received  by  the  first  compensator  will  be  telling  the  compensator  to  drive  the  plant 
to  a  point  outside  of  its  re^on  of  attraction.  Therefore,  the  compensator  output  will  not  be 
valid,  as  it  will  be  based  on  the  assumption  that  the  LTV  plant  does  not  change  between 
K  =  1  and  y  =  1.2. 

These  considerations  account  for  the  poor  performance  shown  in  Figure  3.20. 

These  Umitations  cannot  be  completely  overcome  using  linear  control  theory.  Therefore,  some 
degree  of  uncertainty  will  be  inevitable  in  a  banked  controller.  It  will  be  the  function  of  FWzy  Lo^c 
within  this  controller  to  account  for  these  shortfalls  in  linear  theory.  An  interesting  extension  to 
the  above  simulations  is  to  allow  the  r  upon  which  the  compensator  is  based  vary  as  a  function  of 
the  plant  output.  This  approach  leads  to  time- varying  linear  compensators  and  is  explored  in  more 
detiul  in  the  next  chapter 
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S.7  Summary 

This  chapter  demonstrates  that  banked  compensation  is  possible  in  the  linear  sense  if  each 
compensator  in  the  bank  is  able  to  account  in  some  way  for  the  energy  it  injects  into  the  plant  while 
in  contrtd.  This  can  be  accomplished  either  by  prewmghting  the  compensators  and  allowing  each 
compensator  to  affect  the  plant  for  all  time,  or  by  post-weighting  the  compensators  and  allowing 
energy  to  build  up  in  inheriting  compensators.  For  pre-weighting,  the  compensator  input  must  be 
driven  to  zero  in  order  to  counter  the  tendency  of  later  compensators  to  overdrive  the  system.  For 
post-weighting,  each  compensator  is  receiving  an  identical  error  signal,  so  every  compensator  will 
have  built  up  internal  state  energy  at  the  time  of  transfer. 

Banked  compensation  only  proves  successful  for  the  simple  case  of  a  linear  plant  and  identical 
linear  compensators.  Under  these  conditions  the  resulting  performance  will  be  the  same  as  that 
of  a  single  compensator  in  series  with  the  plant.  Undesired  dynamics  are  induced  by  dissimilar 
compensators,  time-varying  plants,  and  modeling  errors.  A  prototype  non-F\izzy  banked  compen¬ 
sator  using  dissimilar  compensators  to  drive  a  time-varying  plant  was  unsuccessful  due  to  these 
unmodeled  effects. 
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IV.  Model-Based  Fuzzy  Logic  Control  Nonlinear  Design  Considerations 
4‘1  Introduction 

The  last  chapter  developed  a  linear,  non-Foszy  Banked  Compensator  which  performed  wdl 
for  linear,  time-invariant  plants.  However,  the  compensator  was  inadequate  when  applied  to  the 
LTV  plant,  a  simplified  version  of  the  nonlinear  plant.  It  is  clear  that  a  more  sophisticated  struc¬ 
ture  is  required,  involving  the  analysis  of  both  the  nonlinearities  exhibited  by  the  plant,  and  the 
effectiveness  of  linear  control  approaches  at  addressing  these  nonlinearities. 

This  chapter  will  addressed  the  following  areas  of  concern: 

1.  Effectiveness  of  Linearized  Plants  for  Nonlinear  Control 

2.  Time- Varying  Compensation  of  the  Linear  Hme- Varying  Plant 

3.  Time- Varying  Compensation  of  Nonlinear  Plant 

4.  Banked  Compensator  Approximation  of  Time- Varying  Compensation 

5.  Fuzzy  Weighted  Bank  Control  of  Nonlinear  Plant 

6.  Summary 

4.i  Effectiveness  of  Linearized  Plants  in  Nonlinear  Control 

Lineair  perturbation  models  are  routinely  used  by  control  engineers  to  develop  controllers  for 
nonlinear  systems  [17].  However,  it  is  reasonable  to  question  just  how  effective  controllers  based  on 
perturbation  models  actually  are.  Linear  analysis  will  certainly  be  valid  for  operating  conditions 
very  close  to  the  equilibrium  condition  upon  which  the  linearization  is  based,  but  the  actual  r^on 
of  linear  behavior  will  be  heavily  dependent  on  the  the  extent  of  the  underlying  nonlinearity. 

To  address  this  question,  the  concept  of  region  of  attraction  has  been  introduced  for  this 
thesis.  The  region  of  attraction  is  a  sphere  in  state  space  within  which  the  dynamics  demonstrated 
by  the  linearized  model  represent  a  reasonable  approximation  of  the  behavior  of  the  nonlinear  plant. 
The  r^ons  of  attraction  for  the  linearized  model  and  for  the  compensator  based  on  that  modd 
are  considered  one  and  the  same.  For  example,  consider  the  nonlinear  state  equation  y  =  sinx.  In 
this  case,  y  may  be  assumed  linear  for  |z|  <  x/6  (y  tv  x).  Therefore  the  r^on  of  attraction  for  the 
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linearized  model  and  for  the  compensator  based  on  that  model  would  be  |z|,|y|  <  Of  coarse, 
this  nonlinearity  is  much  more  benign  than  that  considered  in  this  thesis.  The  regions  of  attraction 
for  the  nonlinear  plant  under  consideration  will,  in  general,  be  much  smaller. 

The  boundaries  of  the  re^ons  of  attraction  will  be  based  on  the  controller  being  capable  of 
delivering  good  performance,  with  "good”  defined  as: 

•  Within  5%  of  correct  overshoot  (1.064  <  Mp  <  1.176). 

•  Within  4%  of  final  value  within  specified  settling  time  {t,  =  1.62  secs). 

By  analyzing  Equations  (4.11)  and  (4.12)  it  is  reasonable  to  assume  that  linearized  plants 
would  become  more  effective  as  the  equilibrium  values  become  larger.  As  Xio,  Xjo,  Ug  get  larger,  the 
perturbations  contribute  a  smaller  percentage  of  the  overall  state  energy.  As  =  Ug  -*  0 

the  linearized  model  should  become  essentially  worthless  at  estimating  the  response  of  the  nonlinear 
plant.  This  is  borne  out  in  simulation. 

As  the  objective  is  to  drive  the  nonlinear  plant  from  0  to  1,  linearized  plants  in  this  region 
should  be  examined  first.  In  each  case  the  effectiveness  of  the  compensator  at  a  ^ven  equilibrium 
condition  is  determined  as  follows: 

1.  A  linear  compensator  is  developed  using  Equation  (3.2)  for  the  appropriate  linearized  plant. 

2.  The  compensator  is  simulated  in  closed-loop  with  the  nonlinear  plant. 

3.  The  simulation  is  set  up  to  begin  from  rest  conditions  on  or  near  the  equilibrium  point  for 
the  plant  on  which  the  compensator  is  based. 

4.  Successively  larger  steps  away  from  the  equilibrium  are  attempted. 

5.  The  largest  step  input  with  good  performance  is  recorded  as  the  region  of  attraction  for  that 
compensator. 

Table  4.1  shows  the  results  these  tests.  For  each  equilibrium  point,  three  values  are  ^ven: 
originating,  terminating,  and  traversing.  These  labels  refer  to  starting  point  of  the  simulation 
versus  the  equilibrium  point  of  the  compensator.  For  originating  simulations,  the  simulation  starts 
from  the  equilibrium  point  and  steps  up.  For  terminating  simulations,  the  simulation  starts  from 
an  equilibrium  below  the  compensator  and  attempts  to  step  up  the  equilibrium  for  which  the 
compensator  is  designed.  For  traversing,  the  simulation  attempts  to  drive  the  system  from  an 
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Equilibrium 

Originating 

Region  of  Attraction 

Terminating 

IVaversing 

0.1 

0.000025 

0.000027 

0.000025  1 

0.2 

0.00022 

0.00024 

0.00023  1 

lEHIHH 

0.00082  1 

10.4 

0.0022 

0.0024  1 

lESHI 

0.004 

0.0043 

0.0042 

0.0075 

0.0078 

0.0076  1 

0.7 

0.0097 

0.0097 

0.0097  1 

0.8 

0.01 

0.012 

0.01  1 

0.9 

0.013 

0.014 

0.013  1 

1.0 

0.018 

0.02 

0.018  1 

1.1 

0.02 

0.03 

0.022  1 

1.2 

0.03 

0.05 

0.038  1 

Table  4.1  Redons  of  attraction  for  linear  compensators  for  originating,  terminating,  and  traverse 
modes  of  operation 


equilibrium  below  the  compensator  to  a  final  steady-state  value  above  the  equilibrium  point  for 
the  compensator.  In  all  traversing  simulations,  the  originating  and  terminating  equilibria  are 
symmetrically  distributed  about  the  compensator  equilibrium. 

Prior  to  the  0.7  equilibrium  point,  overshoot  is  the  limiting  factor.  After  0.7,  settling  time 
determines  region  of  attraction.  Beyond  approximately  Y  =  1.4,  system  response  will  always  be 
good  by  the  established  criteria.  However,  the  response  wiU  be  highly  nonlinear  for  large  step  inputs. 
Appearance  of  the  closed-loop  response  will  be  the  determining  factor  for  regions  of  attraction  above 
1.4. 

Notice  in  the  table  that  in  all  cases  originating,  terminating  and  traversing  values  are  very 
similar.  This  indicates  that  the  model  is  good  within  the  region  of  attraction,  regardless  of  where 
you  start  and  finish  within  the  region.  Because  the  "originating”  column  is  generally  the  most 
restrictive,  this  value  will  be  used  when  referring  to  re^on  of  attraction  in  the  future. 

The  number  of  compensators  required  for  a  Banked  MBFLC  is  not  clear.  For  initial  designs, 
the  number  will  be  established  by  requiring  that  at  least  one  compensator  be  considered  valid  at 
any  point  between  the  origin  and  reference  value.  Table  4.1  shows  that  the  number  of  compensators 
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T^ble  4.2  Redons  of  attraction  for  linear  compensators  in  Banked  MBFLC 

required  to  ensure  this  can  be  very  large,  especially  near  the  origin  where  the  strong  nonlinearity 
forces  the  regions  of  attraction  to  zero.  Perhaps  more  importantly,  the  it^ons  of  attraction  impose 
severe  limits  on  the  magnitude  of  the  error  signal  for  which  the  compensator  will  respond  properly. 
Linear  compensators  near  the  origin  will  provide  very  little  information  as  to  the  dynamics  of  the 
nonlinear  system. 

Because  of  these  factors,  initial  Banked  MBFLC  development  will  take  place  away  from  the 
nonlinear  distortions  near  the  origin.  The  system  will  be  driven  from  P  =  3  to  T  =  4,  rather  than 
from  Y  =  0  to  y  =  1.  This  provides  for  more  effective  and  homogeneous  linearized  modds.  The 
results  obtained  within  this  range  will  then  be  reproduced  over  a  smaller  range  closer  to  the  origin. 
The  re^ons  of  attraction  for  the  r^on  between  T  =  3  and  y  =  4  are  shown  in  Titble  4.2. 

As  mentioned  above,  establishing  the  region  of  attraction  becomes  more  objective  as  Y  gets 
larger.  As  the  system  moves  further  and  further  from  the  origin,  the  compensators  require  less 
gain  and  so  fail  less  dramatically.  For  large  step  inputs  the  system  response  gets  very  distorted  but 
never  actually  violates  the  boundaries  established  for  a  "good”  system  response  in  the  last  section. 
Figure  4.1  shows  steps  of  varying  magnitudes  away  from  the  equilibrium  pmnt  at  y  =  4.0.  At  no 
time  does  the  response  violate  any  of  the  good  performance  boundaries  (shown  in  the  Figure),  but 
the  quality  of  the  response  drops  off  significantly.  In  order  to  minimize  unwanted  dynamics,  the 
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Figure  4.1  Determination  of  region  of  attraction  for  Y  —  4.0 

region  of  attraction  for  Y  =  4.0  was  set  at  0.55,  as  given  in  the  'Dtble.  Similar  trade-offs  were  made 
for  all  other  equilibrium  points. 

The  compensators  can  only  be  considered  "good”  as  long  as  they  do  not  receive  error  signals 
larger  than  their  regions  of  attraction.  Error  signals  larger  than  the  region  of  attraction  will 
induce  unmodeled  dynamics  and  oscillations  in  the  system  caused  my  mismodeling  of  the  plant, 
as  discussed  in  Sections  3.6  and  3.7.  Fortunately,  the  undesirable  dynamics  in  the  re^on  between 
y  =  3  and  Y  =  4  will  be  much  slower  and  contmn  less  energy  than  those  induced  by  higher-gain 
compensators  near  the  ori^n.  An  effective  banked  MBFLC  compensator  should  be  capable  of 
quelling  these  dynamics  and  inducing  a  smooth  response  within  the  established  ps^ameters.  An 
effective  design  at  higher  outputs  could  then  be  "slid  down”  to  encompass  lower  r^ons  of  the 
system  output. 
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4.3  Time- Varying  Compensation  of  the  Linear  Time-Varying  Plant 

As  was  mentioned  in  the  last  chapter,  the  Linear  Time- Varying  Plant(LTV  plant)  is  a  sim¬ 
plified,  but  still  nonlinear  repre..entation  of  the  full  nonlinear  plant.  The  LTV  plant  is  analogous 
to  relaxing  the  static  equilibrium  point  assumption  of  the  linearization  procedure.  Control  of  the 
LTV  plant  is  equivalent  to  control  of  the  nonlinear  plant  to  the  extent  that  the  small  perturbation 
assumption  holds. 

The  non-Fuzzy  Banked  Compensator  proved  unable  to  adequately  drive  the  LTV  plant  to  a 
given  reference  value.  An  adequate  controller  can  be  formulated,  however,  by  considering  a  single 
time- varying  compensator  in  the  context  of  nonlinear  control  theory.  The  resulting  controller  can 
then  be  approximated  using  only  linear  control  theory  and  Fuzzy  Logic. 

This  controller  can  be  developed  by  considering  a  compensator  containing  an  inverse  of  the 
linearized  plant  dynamics,  as  was  developed  in  Chapter  3  (3.2).  In  this  approach,  the  compensator 
will  attempt  to  null  out  the  dynamics  generated  by  the  LTV  plant  and  replace  these  dynamics  with 
the  response  of  an  ideal  system  model  trident  in  the  compensator.  The  nulling  portion  of  the 
compensator  is  based  on  canceling  the  poles  of  the  linearized  plant.  Because  the  linearized  plant 
model  varies  according  to  the  equilibrium  point,  the  compensator  should  be  made  time- varying  as 
well.  This  approach  depends  heavily  on  an  accurate  model  of  the  plant,  and  is  generally  not  used 
in  linear  designs.  The  TVcomp  block  with  t  provided  in  real  time  based  on  Y  will  be  referred  to 
as  the  Time- Varying  Compensator  (TVC). 

Figure  4.2  shows  a  SIMULINK  simulation  in  which  a  single  Nulling  Compensator  is  used 
in  the  forward  path.  Internally,  this  compensator  is  the  same  as  the  linear  banked  compensators 
considered  last  chapter,  but  with  the  value  of  r  determined  in  real-time  based  on  the  output. 
Therefore,  at  any  instant  ti  (tj  >  0),  the  appropriate  value  of  t  (r  =  3y(ti)’)  will  be  fed  into  the 
compensator.  This  ensures  a  proper  control  input  to  the  LTV  plant  for  all  time  [30], 

Figure  4.3  shows  the  response  of  this  plant  to  a  step  reference  signal  Ref  =  3  +  «-i(0*  I® 
this  case,  the  response  is  equivalent  to  that  of  the  model  closed-loop  response  M(s),  shown  in 
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Figure  4.2  Control  of  Time- Varying  Linear  Plant  Using  Timf»- Varying  Li.  iear  Compensator 

Figure  2.3.  The  time-varying  compensator  is  capable  of  controlling  the  LTV  plant  throughout  its 
envelope  of  operation  (V  >  0). 
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Tim  {Ha} 

Figure  4.3  Closed-loop  response  of  LTV  plant  and  Hme- Varying  Linear  Compensator 
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4.4  Nulling  Compensation  of  the  Nonlinmr  Plant 

It  has  been  claimed  that  the  LTV  plant  is  the  equivalent  of  relaxing  the  static  equilibrium 
point  assumption  of  the  linearization  process.  If  this  is  true,  then  a  controller  applicable  to  the 
LTV  plant  should  also  be  acceptable  for  the  full  nonlinear  plant,  as  long  the  smaU  perturbation 
assumption  holds. 

Figure  4.4  verifies  that  this  is,  indeed,  the  case.  The  Figure  shows  the  response  of  the  nonlinear 
plant  about  F  =  3  for  various  step  reference  inputs.  For  small  inputs,  where  the  small  perturbation 
assumption  holds,  the  response  of  the  nonlinear  plant  resembles  that  of  the  LTV  plant.  As  the 
step  reference  input  becomes  larger,  the  system  response  decays.  The  next  step  is  to  develop  a 
mathematical  approach  to  overcome  the  limitations  imposed  by  the  small  perturbation  assumption. 

In  order  to  develop  an  overall  control  approach,  the  original  linearization  of  the  nonlinear 
plant  must  be  revisited.  The  nonlinear  plant,  once  again,  is  described  by: 


Xt 

Xi 

Y 

and  the  linearized  plant  is  described  by: 

*1 

*a 

y 


-Xt  +  Xi  (4.1) 

-Xf  +  U*  (4.2) 

Xi  (4.3) 

— *1  +  (4-4) 

—TXi  +  TU  (4.6) 

*1  (4.6) 

3y».  (4.7) 


As  mentioned  several  times,  the  linear  model  is  based  on  two  assumptions:  1)  Small  pertur¬ 
bations  from  the  equilibrium  condition,  and  2)  Static  equilibrium  condition.  As  long  as  a  single 
equilibrium  point  is  chosen,  as  was  the  case  in  developing  the  banked  compensation  approach,  r  is 
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Figure  4.4  Closed-loop  response  of  nonlinear  plant  and  Time- Varying  Linear  Compensator  for 
various  step  inputs 


a  constant  and  the  Laplace  domain  transfer  function  is: 


F(s)  = 


r 


(4.8) 


This  plant  model  can  be  used  as  the  basis  for  a  time- varying  compensator  by  simply  aUowing  r  to 
vary  with  the  desired  change  in  equilibrium.  However,  this  violates  one  of  the  assumptions  upon 
which  (4.8)  is  based,  so  the  linearization  process  should  be  reexamined  with  the  static  equilibrium 
assumption  relaxed.  This  is  easily  accomplished  by  expressing  (2.1)  and  (2.2)  in  terms  ai  the 
equilibrium  conditions  plus  a  perturbation,  as  was  done  in  the  initial  derivation. 


Xi,  -f  xi  =  —Xu  —  +  Xu  +  (d-9) 

Xu  +  *3  =  — (J^l#  +  *l)*  +  (1^»  +  **)*•  (d*10) 
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Expanding  the  cubic  terms  and  rearranging  yidds: 


+  +  *1 

Jfa.  +  *a  =  -jr».*-3Jri,**j-3Jr»^;-xf  +  tf,*  +  3(7,’a  +  3l7,«»+a*.  (4.12) 


Applying  the  small  perturbation  assumption  yidds: 

jti0  +  *1  =  —Xi,  +  Xa«  —  Xi  + 1}  (4*13) 

Xu  +  Xi  =  -jru’-3jr,/x,  +  l7.*  +  3^.V  (4.14) 

Note  that  the  ^i^and  terms  are  nonzero  and  therefore  cannot  be  eliminated  from  4.13  and  4.14. 
However,  the  terms  on  the  right-hand  sides  of  these  equations  invdving  only  the  equilibrium  terms 
still  caned  each  other  out.  This  shows  that  rdaxing  the  static  equilibrium  assumption  induces  no 
additional  dynamics  in  this  particular  plant,  and  (4.8)  holds  true  for  all  r. 

Note  also  that  the  left-hand  side  expressions  in  4.13  and  4.14  ate  equal  to  X\  and  le- 
spectivdy.  This  implies  that,  as  long  as  the  small  perturbation  condition  holds,  any  compensator 
based  on  the  time  varying  plant  will  also  ensure  adequate  performance  for  the  nonlinear  plant.  This 
explains  the  success  of  the  time- varying  compensator  at  controlling  the  LTV  plant. 

For  large  perturbations  away  from  the  equilibrium  condition,  the  small  perturbation  assump¬ 
tion  breaks  down  and  dynamics  due  to  the  nonlinear  terms  in  (4.12)  become  significant.  This  is 
the  cause  of  the  ultimate  failure  of  the  time-varying  compensator  to  contrd  the  nonlinear  plant. 
However,  the  fact  *  'lat.  the  optimal  state  space  trajectories  for  the  states  and  the  perturbations  in 
the  states  are  equivalent  assuming  small  departures  from  the  nominal,  that  is  Xx„«m(0  =  A'in«m(0t 
=  A'2nom(0>  snggcsts  an  important  rdationship.  Namdy,  that  the 

state  trajectory  of  the  time-varying  Unear  modd  will  be  the  same  as  the  state  trajectory  for  the 
nonUnear  model  in  general.  The  inputs  required  to  produce  these  trajectories  in  the  nonUnear  case, 
though,  will  vary  according  the  effect  of  the  nonUnear  terms  in  4.14.  For  small  departures  from 
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equilibrium,  U  fa  u.  At  the  tracking  signal  moves  further  from  the  quiescent  point  of  the  system, 
however,  tt(t)  and  U(t)  diverge  in  a  nonlinear  fashion. 

Fortunately,  the  divergence  of  tt(t)  and  U{t)  can  be  determined  in  real-time  and  corrected 
for.  This  is  accomplished  in  a  two-step  process.  First,  the  desired  dynamics  must  be  determined 
based  on  the  output  of  the  TVC  and  the  known  linearized  model  of  the  nonlinear  plant.  Seomd, 
the  correct  input  U  to  the  nonlinear  plant  must  be  determined  based  on  the  nonlinear  plant  modd. 
These  two  steps  equate  to  a  model  fdlowing  procedure  and  can  be  accomplished  on-line  for  the 
class  of  nonlinear  systems  including  the  plant  currently  under  consideration  [2]. 

An  illustration  of  this  nonlinear  contrd  scheme  is  shown  in  Figure  4.5.  In  this  case,  the 
nonlinear  plant  is  described  by  (2.1)  and  (2.2),  and  the  TVC  is  given  by  (3.2).  The  deured 
dynamics  X]  can  be  determined  by  analyzing  the  second  equation  of  the  linearized  plant  (2.22). 
The  first  equation  need  not  be  considered  as  it  is  not  directly  impacted  by  the  control  input,  only 
by  Xj.  The  desired  x%  at  any  time  t  can  then  be  determined  by  measuring  Y  to  determine  r  and 
Xi.  These  quantities,  along  with  the  tt  produced  by  the  TVC  can  be  used  to  sdve  for  the  desired 
X]  using  2.22.  The  SIMULINK  block  performing  this  function  is  shown  in  Figure  4.6. 

It  is  tempting  to  argue  that  the  actual  form  of  the  desired  dynamics  is  simply  Xj  s  ru  because 
the  time-varying  nature  of  the  compensator  ensures  that  Xi  =  0  for  all  time.  However,  simulation 
shows  that  the  system  cannot  be  considered  to  have  actually  shifted  equilibrium  points  until  the 
dynamics  go  to  zero.  Therefore,  though  the  determination  of  the  ideal  state  trajectory  is  based 
on  a  TVC  which  assumes  a  changing  equilibrium  point,  the  — rxi  term  in  2.22  is  still  required  to 
provided  adequate  retarding  of  the  X]  state  model  as  the  desired  final  state  is  reached. 

The  final  step  is  to  solve  for  the  required  Cf,  and  this  is  accomplished  through  (2.2).  Rear¬ 
ranging  this  equation  to  solve  for  V ,  yields: 

V  =  ^(Xa-l-A^,’).  (4.15) 
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Figure  4.5  F\ill-envelope  nonlinear  control  scheme  using  trajectory  following 
By  assuming  that  Xi(t)  =  Xi{t)f  as  was  argued  above,  the  result  is: 

U  =  ^(ar'a  +  A-,*).  (4.16) 

The  SIMULINK  implementation  of  (4.16)  is  shown  in  Figure  4.7.  Due  to  a  difficulty  in 
getting  MATLAB  to  perform  the  cubed-root  operation  when  the  operand  was  negative,  a  modified 
simulation  was  necessary.  The  two-branch  approach  to  taking  the  cubed-root  when  the  operand  is 
positive  or  negative  overcomes  this  limitation  with  MATLAB. 

SIMULINK  runs  demonstrate  that  this  system  performs  precisely  as  required  by  the  system 
specifications.  Figures  4.8  and  4.9  show  the  results  of  simulations  starting  from  various  initial 
conditions  and  with  varying  step  sizes.  In  all  cases,  the  specifications  are  met  to  the  d^ree  of  the 
accuracy  of  the  TVC  design  for  the  time-varying  linear  system. 

This  system  also  has  disturbance  rejection  capabilities.  Figure  4.10  shows  a  i^ulation  prob¬ 
lem  with  a  unit  step  disturbance  introduced  at  the  plant  input.  Figure  4.11  shows  the  system 
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Figure  4.6  Determination  of  the  desired  Xa  based  on  compensator  outputs 


performance  based  on  SIMULINK  simulation.  This  plot  also  shows  the  disturbance  rejection  capar 
bilities  of  the  nonlinear  system  when  a  time- varying  robust  linear  compensator  is  substituted  fw  the 
non-robust  compensator.  The  robust  compensator  performs  as  well  as  the  nonrobust  coinpensator 
for  both  tracking  and  regulation  in  this  contrd  configuration. 
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Figure  4.7  determination  of  U  based  on  desired  Xj 
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Figure  4.8  Closed-Loop  simulations  of  nonlinear  controller  using  time-varying  linear  compensa¬ 
tion,  starting  from  F  =  0 
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Figure  4.9  Closed-Loop  simulations  of  nonlinear  controller  using  time-varying  linear  compensa¬ 
tion,  starting  from  Y  =  I 
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Figure  4.10  SIMULINK  simulation  to  test  the  disturbance  rejection  capabilities  of  nonlinear 
trd  scheme 


Figore  4.11  Disturbance  rejection  capabilities  of  nonlinear  controller  for  unit  step  lifted  at  plant 
input 
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4-5  Banked  Compensator  Approximation  of  Time- Varying  Compensation 

The  above  controller  does  meet  desired  performance  expectations,  but  at  the  cost  of  all 
established  design  limitations.  The  stated  objective  of  this  research  is  to  create  a  controller  nsing 
only  linear  compensators  and  Fuzzy  Logic.  The  challenge,  then,  is  to  approximate  the  above  system 
using  only  these  elements. 

The  first  step  in  accomplishing  this  will  be  to  eliminate  the  need  for  the  time-varying  com¬ 
pensator.  This  wiU  be  accomplished  by  using  a  bank  of  linear  compensators  weighted  by  Fuzzy 
Logic.  Figure  4.12  shows  a  simple  compensator  bank  for  a  Model- Based  Fuzzy  Lope  Controller. 
In  this  case,  five  linear  compensators  are  being  used,  corresponding  to  the  equilibrium  points 
Y  =  3,y  =  3.35,  Y  =  3.5,  Y  =  3.75  and  Y  =  4.  The  optimal  number  of  compensators  required 
for  a  given  application  depends  on  the  nonlinearities  exhibited  by  the  plant  between  the  starting 
output  and  the  objective  output.  This  will  be  explored  in  more  detail  in  the  next  chapter. 

When  the  system  is  at  V  =  3  or  V  =  4,  there  is  no  doubt  as  to  which  compensator  should  be 
applied  to  the  plant.  The  situation  is  not  so  clear  when  the  system  is,  for  example,  at  T  =  3.887. 
In  this  case,  the  best  output  would  be  a  mix  of  both  compensators.  This  is  the  function  of  the 
Fuzzy  Supervisor  shown  in  the  Figure. 

The  internal  structure  of  the  Fuzzy  Supervisor  block  is  shown  in  Figure  4.13.  It  is  referred 
to  as  a  supervisor  rather  than  a  controller  because  it  serves  only  a  very  specific  function:  To 
emphasize  the  compensator  which  most  nearly  reflects  the  true  state  of  the  system.  It  has  the 
additional  feature  that,  when  the  system  is  between  regions  of  attraction,  it  will  output  a  weighted 
average  of  the  control  actions  advocated  by  each  relevant  compensator. 

The  weighting  of  the  compensators  can  occur,  as  discussed  last  chapter,  either  before  or  after 
the  compensators  themselves.  Both  approaches  have  advantages,  but  post- weighting  proves  superior 
for  hybrid  FNizzy/nonlinear  applications.  The  reason  is  simple:  there  is  very  little  uncertainty  in 
this  system.  Preweighting  makes  for  a  more  volatile  bank  of  compensators  and  allows  the  Fuzzy 
Supervisor  more  control  authority.  In  this  application,  however,  the  objective  is  simply  to  reproduce 
the  time-varying  compensator  as  faithfully  as  possible. 
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Figure  4.12  Simple  compeusator  bank  for  MBFLC 

The  precise  operation  of  the  Fuzzy  Supervisor  is  as  follows: 

1.  The  current  value  of  Y  is  fed  into  the  Fuzzy  Supervisor. 

2.  Y  is  evaluated  for  membership  into  the  FWzy  Sets  Y  =  3,y  =  3.35,  Y  =  3.5,  Y  =  3.75  and 
y  =  4.  The  form  of  the  memberships  functions  for  these  sets,  given  in  (2.38)  and  repeated 
here,  is: 

p(*)  =  e-^,  (4.17) 

These  functions  are  evaluated  over  the  universe  of  discourse,  given  the  means  and  variances 
shown  in  Figure  4.13,  are  shown  in  Figure  4.14. 

3.  In  the  block  Normx  the  activation  value  of  each  FHizzy  Set  is  divided  by  the  total  activation 
of  all  IHizzy  Sets  in  the  the  Fuzzy  Supervisor.  This  scales  all  wmghting  outputs  so  that  thdr 
sum  is  one. 

4.  The  normalized  weights  are  output  to  the  compensator  bank.  The  activation  of  the  f\izzy 
Set  y  =  3  is  used  to  weight  the  compensator  for  r  =  27.  The  activation  of  y  =  4  is  used  to 
weight  the  input  to  the  compensator  for  r  =  48,  etc. 
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Figure  4.13  Internal  Structure  of  F\izzy  Supervisor 

This  procedure  is  the  same,  regardless  of  the  number  of  compensators  the  FHizzy  Supervisor 
is  wdghting.  The  Fuzzy  supervisor  could  be  replaced  by  a  look-up  table  in  final  implementation, 
provided  no  on-line  retnning  is  desired. 


Figure  4.14  Activation  of  f\izzy  Membersliip  functions  for  F\izzy  Supervisor 
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A  comparison  of  the  Fuzzy  Weighted  Bank  and  the  time-varying  compensator  is  shown  in 
Figure  4.15.  This  plot  shows  the  closed-loop  compensator  output  histories  (u)  for  the  time- varying 
compensator,  the  Fuzzy  Weighted  Compensator  Bank  and  a  single  linear  compensator  for  r  =  27. 
The  time- varying  compensator  and  Fuzzy  Weighted  Bank  generate  almost  identical  u  time  histories, 
and  therefore  are  virtually  indistinguishable  in  closed-loop  simulation.  The  response  of  the  nonlinear 
plant  in  closed-loop  simulation  with  the  Fuzzy  Weighted  Compensator  Bank  is  shown  in  Figure  4.16. 
The  response  of  the  ori^nal  time- varying,  nonlinear  system  over  the  same  interval  is  also  shown  in 
the  Figure,  though  they  are  almost  identical. 

Figure  4.17  shows  the  state  space  trajectory  of  the  hybrid  nonlinear/Fuzzy  Bank  controller. 
This  plot  is  very  similar  to  the  trajectory  generated  by  the  dynamic  inversion-based  compensator 
in  Chapter  2.  Notice  that,  after  an  initial  jump  in  the  value  of  ff,  the  system  departs  from  the  locus 
of  equilibrium  points  (shown  as  a  dotted  line),  and  does  not  return  until  the  system  is  approaching 
the  desired  final  equilibrium.  This  form  of  the  response  is  reasonable,  as  the  dynamics  of  the 
plant  go  to  zero  as  the  state  approaches  the  equilibrium  locns.  However,  Figure  4.17  predicts 
that  linear  compensators  intermediate  to  the  initial  and  final  conditions  of  the  simulation  will  not 
provide  information  relevant  to  control  of  the  plant.  The  optimal  state  trajectory  will  be  outside 
the  regions  of  attraction  for  these  compensators.  This  is,  indeed,  the  case,  as  the  final  form  of  the 
MBFLC  (given  in  Chapter  1)  utilizes  the  equivalent  of  only  two  compensators. 
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Figure  4.15  Time  history  of  «  outputs  for  time-varying  compensator  (TVC),  frizzy  Compensated 
Bank,  and  single  linear  compensator.  The  TVC  and  single  linear  compensator  plots 
and  nearly  identical. 
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^.6  Summary 

In  this  chapter,  the  effectiveness  of  linear  controllers  applied  to  nonlinear  {dants  was  quanti¬ 
fied.  Regions  of  attraction  were  established  to  estimate  the  number  of  compensators  that  mif^t  be 
required  in  a  compensator  bank. 

A  time- varying  linear  control  strategy  was  develop  to  adequatdy  compensate  the  LTV  |daat. 
The  successful  design  was  then  applied  to  the  nonlinear  plant.  It  was  found  that  the  time- varying 
linear  approach  was  effective  for  the  nonlinear  plant  for  step  inputs  with  small  amplitudes.  Ad¬ 
ditional  nonlinear  processing  was  developed  to  create  a  full-envelope  nonlinear  contndler  for  the 
nonlinear  plant. 

The  time- varying  linear  compensator  in  the  nonlinear  controller  was  then  replaced  with  a  bank 
of  linear  compensators  pre- weighted  with  F\izzy  Logic.  The  performance  of  this  hybrid  compensator 
was  found  to  be  indistinguishable  from  the  oripnal  time-varying  linear  case. 

The  state  trajectory  generated  by  the  nonlinear/Fuzzy  Bank  contrtdler  is  similar  to  the  tra¬ 
jectory  based  on  other  nonlinear  techiuques.  Tlus  suggests  that  the  optimal  trajectory  path  will 
lead  outside  the  redone  of  effectiveness  of  intermediate  linear  compensators,  and  favors  a  two- 
compensator  bank. 
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V.  Development  of  Model-Baaed  Pozzy  Logic  Controller 


5.1  Introduction 

A  controller  was  developed  last  ch^ter  which  ntilues  both  linear  system  modds  and  FHizsy 
Logic  to  control  a  nonlinear  plant.  Unfortunately,  the  controller  also  nltimatdy  relied  on  a  nonlinear 
model  of  the  plant  to  convert  the  xj  trajectory  produced  by  the  Frizzy  Weighted  Bank/zHtfoLtn^ 
blocks  to  a  U  plant  input.  The  linear/F\izzy  dements  of  the  controller  serve  only  to  provide  a 
’’modd”  for  the  plant  to  follow.  The  actual  contrd  signal  was  determined  using  nonlinear  mathe¬ 
matics  and  an  assumed  equivalence  between  the  dynamics  of  the  linearized  plant  and  the  nonlinear 
plant  starting  from  equilibrium.  It  should  further  be  noted  that  the  Xj  history  can  be  produced 
using  other  approaches  as  well.  In  fact,  any  mechanism  which  produces  a  valid  state  trajectory 
could  be  inserted  in  the  controller  in  place  of  the  linear/Fnzzy  dements. 

The  objective  of  this  chapter,  then,  is  to  remove  all  nonlinear  processing  dements  from  the 
controller  other  than  Fuzzy  Sets  or  equivalent  look-up  tables  based  on  F^zy  representations.  This 
objective  can  be  accomplished  in  two  ways: 

•  Approximate  the  functionality  of  the  Interpreter  block  using  FHizzy  Lope.  This  could  be 
accomplished  using  the  technique  d  Wang  [23]  and  would  result  in  a  contrdler  very  similar 
to  Figure  4.5,  with  the  nonlinear  processing  block  replaced  by  a  FHizzy  Logic  approximation. 

•  Approximate  the  entire  controller  as  Fuzzy  with  implicates  of  the  FVizzy  IF...THEN  rules 
supplied  by  the  Fuzzy  Weighted  bank.  In  this  approach,  the  structure  of  the  contrdler  is 
essentially  that  of  a  FNizzy  Logic  Ck>ntroller.  The  recommended  control  actions,  however,  are 
supplied  by  linear  compensators,  rather  than  a  look-up  table. 

Examples  of  Frizzy  Lo^c  Approximations  of  nonlinear  functions  can  be  found  in  the  literature 
[23].  Therefore,  this  thesis  will  follow  the  second  route  to  obtain  an  MBFLC. 

This  chapter  contains  the  following  sections: 

1.  Structure  of  Banked  Model- Based  FHizzy  Logic  Controller 

2.  FHill  Envelope  Banked  Modd-Based  FHizzy  Logic  Control 

3.  Contribution  of  Linear  Compensators  to  MBFLC  Design 

4.  Modd-FoUowing  Hybrid  Compensators 

5.  Summary 
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S.t  Structure  of  Banked  Model- Based  Fussy  Logic  Controller 

The  objective  of  the  MBFLC  is  to  apply  linear  compensators  to  the  control  of  a  nonlinear 
plant,  using  Fkizzy  Logic  to  account  for  the  errors  committed  by  ^>plying  linear  control  the(»y  to 
a  nonlinear  problem.  By  identifying  specific  linear/nonlinear  disjunctions  in  the  last  turn  ch^ters, 
the  development  of  the  necessary  f\izzy  Logic  elements  is  straightforward.  The  operation  ol  the 
controller  is  easily  summarized:  When  agreement  between  linear  and  nonlinear  theory  is  good,  then 
apply  the  control  recommended  by  the  linear  compensators.  When  agreement  is  poor,  modify  the 
control  signal  to  ensure  that  the  system  responds  with  caution. 

There  are  three  sources  of  error  that  need  to  be  addressed  in  the  MBFLC  design: 

1.  The  state  trajectory  which  produces  the  desired  closed-loop  response  exists  outside  of  the 
regions  of  attraction  for  compensators  intermediate  to  the  initial  and  final  conditions. 

2.  The  error  signal  is  larger  than  the  compensator  re^ons  of  attraction. 

3.  The  system  will  not  be  at  equilibrium  when  the  output  reaches  the  desired  final  output  due 
to  linear/nonlinear  mismatches. 

One  approach  to  overcome  these  errors  would  be  to  use  a  dedicated  F\izzy  Logic  element  for 
each  compensator.  In  this  way,  the  bank  would  consist  of  a  number  of  Linear  Compensator/FHizzy 
Compensator  pairs,  where  each  pair  would  operate  independently,  in  the  spirit  of  Figure  3.12. 
However,  the  analysis  of  banked  compensation  in  Chapter  3  does  not  support  this  approach.  Sim¬ 
ulations  show  that  reasonable  performance  of  individual  compensators  over  small  regions  does  not 
necessarily  imply  good  performance  over  larger  r^ons  using  Banked  Compensation.  For  the  linear 
design  presented  in  Chapter  3,  the  Banked  Compensator  actually  did  a  poorer  job  than  would  have 
any  individual  compensator  acting  alone. 

This  motivates  the  designer  to  consider  the  MBFLC  as  a  single  compensator,  rather  than  as 
a  confederation  of  individual  compensators.  A  top-level  perspective  is  shown  in  Figure  5.1.  The 
Banked  MBFLC  requires  the  inputs  of  both  a  FVizzy  Logic  Controller  {E  and  E)  and  a  linear  banked 
controller(£!  and  Y).  The  controller  will  be  of  the  same  basic  structure  as  the  FLC  in  Chapter  2, 
but  with  one  important  difference:  The  required  control  actions  are  obtained  from  a  bank  of  linear 
compensators,  rather  than  from  a  look-up  table. 
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Figure  5.1  B'nkod  Model-Based  Faz^  Logic  Contnd  paradigm 


By  conttdering  the  Banked  MBFLC  as  a  dngle  unit,  each  error  source  can  be  addressed 
independently,  rather  than  on  a  compensator-by-compmisator  basis.  This  is  the  approadi  adopted 
in  this  research.  EUuJi  of  the  error  sources  above  is  discussed  below  to  determine  an  appropriate 
control  strategy  by  application  of  Fuzzy  Logic  prindples.  The  universe  of  discourse  f<v  MBFLC 
development  was  chosen  to  be  from  F  =  0.6  to  y  =  0.7,  so  that  nonlinear  effects  would  be  evident. 

5.2.1  Limited  Effectiveness  of  Intermediate  Linear  Compensators.  It  was  hoped  that  the 
system  could  be  driven  in  a  linear-like  manner  using  only  information  available  from  linearized 
models  of  the  nonlinear  plant.  Nonlinear  analysis  last  ch^ter,  however,  demonstrated  that  the 
deured  state  trajectory  for  the  nonlinear  plant  under  conmderation  exists  outside  the  r^ons  of 
attraction  for  intermediate  compensators.  FHirther,  the  fact  that  the  system  is  away  from  equilib¬ 
rium  throughout  the  transition  creates  an  energy  mismatch  between  compensators  and  the  plant, 
as  discussed  in  Chapter  3. 

The  success  of  the  MBFLC,  then,  for  any  ^ven  transition,  will  hinge  on  three  elements:  An 
initiating  compensator,  a  terminating  compensator  and  a  FHizzy  Logic  Controlled  transition  between 


the  two.  At  this  point  it  is  assumed  that  the  required  initiating  and  terminating  contrdllers  are 
the  linear  compensators  based  about  the  appropriate  equilibrium  values.  This  assumption  will  be 
tested  in  this  chapter,  and  will  ultimatdy  prove  false. 

Because  the  MBFLC  is  envisioned  as  a  fulUenvelope  contrdler,  many  linear  compensators 
will  be  required  to  adequately  represent  all  pcMsible  combinations  of  initial  and  final  conditions. 
In  theory,  any  number  of  compensators  could  be  developed  for  a  desired  envelope.  For  this  thesis, 
the  nimber  of  linear  compensators  used  is  determined  by  the  desire  to  have  at  least  one  valid 
compensator  at  any  point  between  the  starting  condition  and  the  final  condition.  FYom  Ihble  4.1, 
the  re^ons  of  effectiveness  between  Y  =  0.6  and  Y  =  0.7  are  about  0.009.  Therefore,  to  cover  a 
universe  of  discourse  of  0.1,  approximately  11  compensators  are  required.  This  number  was  rounded 
to  ten  for  ease  of  implementation. 

To  determine  the  appropriate  control  action  for  regions  of  operation  between  compensator 
equilibrium  values,  a  Fuzzy  Supervisor  will  be  utilized.  The  Fuzzy  Supervisor  was  developed  last 
chapter  to  approximate  the  time-varying  linear  compensator.  Though  linear  analysis  demonstrated 
that  pre-weighting  the  compensators  provided  more  tolerance  for  modeling  errors,  nonlinear  anal¬ 
ysis  showed  that  post-weighting  better  emulated  the  time-varying  compensator.  Therefore,  post¬ 
weighting  will  be  used  for  this  iteration  of  the  MBFLC. 

Figure  5.2  shows  a  compensator  bank  for  a  Model-Based  Ihzzy  Lo^c  Controller  between 
Y  =  0.6  and  Y  =  0.7.  The  internal  structure  of  the  FVizzy  Supervisor  for  this  application  is  the 
same  as  that  shown  in  Figure  4.13,  but  with  10  FVizzy  Sets  being  evaluated  simultaneously.  The 
tuning  values  for  the  Fuzzy  Sets  are  given  in  Table  5.1. 

If  the  Fuzzy  Weighted  Compensator  Bank  is  simulated  in  closed-loop  with  the  nonlinear  plant, 
the  result  is  shown  in  Figure  5.3.  Notice  that  overshoot  is  too  large  and  the  compensator  does  not 
cause  the  plant  to  settle  out  within  the  required  settling  time  of  1.62  seconds.  These  errors  are  due 
to  all  three  of  the  sources  mentioned  above.  To  determine  the  effect  of  intermediate  compensators  on 
the  response,  this  simulation  was  repeated  with  all  but  the  initiating  and  terminating  compensators 
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Table  5.1  Tuning  Parameters  for  Fuzzy  Supervisor  Membership  Fhinction 


disconnected.  This  had  only  a  slight  settling  effect  on  the  system,  so  the  other  errors  must  dominate. 
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with  universe  of  discourse  from  V  =  0.6  to  K  =  0.7 


'niiie(aocs) 

Figure  5.3  closed-loop  response  of  nonlinear  plant  when  driven  by  Frizzy  Weighted  Compensator 
Bank  from  Y  =  0.6  to  V  =  0.7 
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Fuzzy  Set 

C2S1 

Variance 

Implicate(  Penalty) 

E  is  0.0 

0.000021 

1.0 

E  is  0.014 

0.014 

0.000021 

0.6 

E  is  0.028 

0.000021 

0.56 

E  is  0.042 

0.042 

0.000021 

0.48 

E  is  0.056 

0.056 

0.000021 

0.44 

E  is  0.07 

0.07 

0.000021 

0.41 

E  is  0.084 

0.084 

0.000021 

0.38 

E  is  0.098 

0.098 

0.000021 

0.37 

E  is  0.112 

0.112 

0.000021 

0.35 

E  is  0.126 

0.126 

0.000021 

0.34 

Table  5.2  Tuning  Parameters  for  Fuzzy  Limiter  Membership  Functions 


S.S.S  Pull  Error  Signal  Input  to  Linearized  Plants.  The  error  signal  E  is  based  on 
error  between  current  and  desired  position.  By  feeding  E  directly  Into  a  compensator  the  tacit 
assumption  is  made  that  the  compensator  being  driven  is  valid  between  the  initial  state  and  the 
final  state.  In  the  case  of  Banked  Compensation,  this  assumption  is  not  true,  and  feeding  the  entire 
E  signal  into  any  single  compensator  contributes  significantly  to  the  errors  evident  in  Figure  5.3. 
To  overcome  this,  E  must  be  scaled  to  account  for  the  limited  validity  of  the  compensators  in  the 
bank.  This  is  the  purpose  of  the  fStzzy  Limiter. 

The  location  of  the  Fuzzy  Limiter  is  shown  in  Figure  5.4.  It  would  be  natural  to  place 
a  Limiter  immediately  before  each  linear  compensator  in  the  bamk.  In  this  case,  though,  given  a 
fixed  starting  point  {Y  =  0.6),  each  error  magnitude  will  be  associated  with  a  specific  compensator. 
Therefore,  a  single  limiter  can  be  constructed,  eliminating  some  complexity.  This  consideration  sets 
a  non-rigorous  lower  limit  on  the  number  of  Fuzzy  Sets  required  vnthin  the  Limiter:  There  should 
be  at  least  one  Fuzzy  Set  for  each  linear  compensator  in  the  Bank. 

The  internal  Structure  of  the  Fuzzy  Limiter  block  is  shown  in  Figure  5.5.  The  tuning  param¬ 
eters  for  the  Fuzzy  Sets  ue  given  in  Table  5.2  Notice  that  this  block  is  not  strictly  Fizzy;  the  final 
output  is  a  direct  function  of  the  input.  In  this  case,  the  Fizzy  Lo^c  is  only  used  to  determine  a 
’’penalty  factor”  for  large  error  signals. 
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The  rules  for  determimug  the  i4>propriate  weighting  factor  have  single  pranises.  The  rules 
represented  in  the  Figure  are: 

•  IF  — Error —  is  0  THEN  Weighting  Fhctor  is  1.0 

e  IF  — Error —  is  0.014  THEN  Wughting  Fhctor  is  0.6 

•  IF  — Error —  is  0.028  THEN  Weighting  Fhctor  is  0.56 

•  IF  — Error —  is  0.042  THEN  Wdghting  Fhctor  is  0.48 

•  IF  — Error —  is  0.056  THEN  Weighting  Fhctor  is  0.44 

•  IF  — Error —  is  0.07  THEN  Weighting  Fhctor  is  0.41 

•  IF  — Error —  is  0.084  THEN  Wdghting  Fhctor  is  0.38 

•  IF  — Error —  is  0.098  THEN  Weighting  Fhctor  is  0.37 

•  IF  — EIrror —  is  0.112  THEN  Weighting  Fhctor  is  0.35 

•  IF  — Error —  is  0.126  THEN  Weighting  Fhctor  is  0.34 

The  operation  of  the  EHizzy  Limiter  is  summarized  as  follows: 

•  The  error  signal  E  is  input  to  the  controller. 

•  The  absolute  value  of  the  error  signal  is  used  to  define  the  universe  of  discourse.  The  Fhzzy 
Limiter  is  only  concerned  with  the  magnitude  of  the  error  signal  and  how  it  compares  with 
the  re^ons  of  attraction  of  the  linearized  compensators. 

•  The  magnitude  of  the  error  is  evaluated  for  membership  in  ten  Fhzzy  Sets.  The  membership 
functions  for  these  sets  are  shown  in  Figure  5.6 

•  The  crisp  weighting  factor  is  determined  by  the  method  of  Mamandi,  which  is  fpven  in  (2.40). 
Simply  stated,  the  membership  value  for  each  set  is  multiplied  by  the  wdghting  factor  implied 
by  that  Set.  These  terms  are  added  together  and  divided  by  the  sum  of  the  membership  values 
for  all  defined  sets. 

•  The  original  error  function,  sign  and  magnitude,  is  multiplied  by  the  wdghting  function. 

It  should  be  mentioned  that  the  objective  of  the  F\izzy  Limiter  is  not  to  place  hard  limits 
on  the  magnitude  of  the  error  function.  This  restriction  would  result  in  a  dosed-loop  system  with 
a  very  slow  response  time.  The  Fuzzy  Limiter  simply  approximates  a  "squashing  function"  which 
scales  the  error  to  acknowledge  the  fact  that  the  linearized  modd  predicts  a  linear  response  to  the 
input  whereas  the  actual  eflfect  of  the  input  will  be  cubic.  The  F\izzy  limiter  could  well  be  a  "Fhzzy 
Booster"  in  situations  where  the  linear  modd  leads  to  under-control  of  the  actual  nonlinear  plant. 

The  values  for  the  penalties  implied  by  each  Fizzy  Set  are  shown  in  Table  5.2.  These  values 
are  based  on  an  estimate  of  the  error  being  introduced  by  a  ^ven  magnitude  Error  dgnal  and  ate 
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refined  through  closed-loop  simulations.  Notice  that  the  penalty  associated  with  the  EitZero  FHissy 
Set  is  1,  implying  no  penalty  for  small  error  signals.  This  will  always  be  the  case  for  a  correctly 
modeled  system. 

The  toning  simulations  start  with  the  smallest  reference  step  away  from  the  initial  equilibrium 
corresponding  to  a  FHizzy  Set  within  the  Limiter.  For  this  step,  the  penalty  factor  is  tuned  so 
that  the  peak  overshoot  exactly  touches  the  Mp  line.  The  process  is  repeated  for  each  step  size 
corresponding  to  a  FHizzy  Set,  until  all  penalties  are  adjusted.  Step  inputs  are  then  simulated  which 
correspond  to  regions  between  the  fuzzy  sets  (such  as  from  an  equilibrium  of  .6043  to  a  final  value 
of  .67).  If  these  steps  do  not  exhibit  the  correct  overshoot,  then  two  actions  can  be  taken  to  correct 
the  problem.  More  FHizzy  Sets  can  be  defined  within  the  compensator,  or  the  variance  term  of 
the  existing  Sets  can  be  altered  to  obtain  adequate  performance. 

The  character  of  the  penalties  depends  on  both  the  location  and  extent  of  the  universe  of 
discourse.  A  universe  of  discourse  far  from  the  ori^n  will  require  smaller  penalties  because  the 
regions  of  attraction  for  these  compensators  are  large.  If  the  universe  of  discourse  itsdf  is  very 
large,  however,  commensurately  large  E  signals  are  possible.  This  could  require  large  penalties  on 
the  E  signal,  even  far  from  the  ori^. 

Figure  5.7  shows  the  closed-loop  response  of  the  nonlinear  plant  to  a  step  input  Rr»f{t)  = 
0.6  -|-  with  the  fHizzy  Limiter  in  place.  It  is  clear  from  the  Figure  that  the  tendency  of 

the  response  to  overshoot  Mp  has  been  eliminated.  All  that  remains  now  are  the  lightly  damped 
oscillations  arising  from  modeling  errors  and  compensator  transfers  occurring  away  from  equilibrium 
points.  This  will  be  addressed  in  the  next  section. 

Notice  that  the  effect  of  the  Fuzzy  Limiter  is  a  time-invariant  function  of  the  Error  and 
the  universe  of  discourse.  The  Fizzy  Limiter  could  be  replaced  with  a  look-up  table  in  final 
implementation,  unless  Fuzzy  Logic-based  optimization  was  desired.  A  concept  for  a  robust  Fizzy 
Limiter  is  presented  at  the  end  of  this  chapter. 

The  Banked  Compensator,  modified  by  the  Fizzy  Supervisor  and  Fuzzy  Limiter,  will  be 
referred  to  henceforth  as  the  Fuzzy  Bank. 
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Figure  5.5  Stmctoie  of  SIMULINK  block  FNizzy  Limiter 
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Figure  5.6  Membership  functions  for  FNizzy  Sets  in  FNizzy  Limiter 


Figure  5.7  Closed-loop  response  of  nonlinear  system  compensated  by  Fuzzy  Bank 
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S.t.S  System  is  Not  at  Equilibrium  when  Output  Reaches  Final  State.  The  Foiiy  Buk 
is  capable  of  producing  the  correct  initial  transient  response  in  the  closed-loop  system,  ipven  an 
adequate  number  of  Sets  within  the  Fuzzy  Limiter.  It  is  unable,  however,  to  settle  out  the  system 
upon  reaching  the  desired  output,  as  shown  in  Figure  5.7.  Dynamics  have  been  generated  in  the 
system  which  are  not  accounted  for  by  the  controller,  the  culmination  of  the  mismodding  between 
the  linear  contndler  and  the  nonlinear  plant.  Farther,  the  error  signal  is  very  small  close  to  the 
input  reference  value.  The  Fuzzy  Limiter,  then,  which  only  weights  the  error  signal,  has  very 
little  contrd  authority.  Therefore,  an  auxiUary  control  input  is  required  to  drive  the  system  to 
steady-state. 

The  undesired  oscillations  can  be  eliminated  by  introducing  an  overall  FSizzy  Logic  architec¬ 
ture  to  the  controller.  Two  rules  are  required  to  ensure  proper  system  response. 

1.  IF  Error  is  Not  Zero  AND  Error  is  not  Positive  THEN  U  =  Uun 

2.  IF  Error  is  Zero  AND  Error  is  Positive  THEN  U  =  Uun  +  C{t) 

As  stated  above,  the  structure  is  essentially  a  Fuzzy  Logic  Controller  with  two  possible  control 
actions:  1)  Output  the  exact  value  of  U  produced  by  the  FHizzy  Bank,  or  2)  Output  the  exact  value 
of  U  produced  by  the  Fuzzy  Bank  p/us  an  auxiliary  C{t).  The  auxiliary  provides  an  additional 
settling  term  when  the  system  is  dropping  towards  the  reference  value  after  the  initial  overshoot. 

The  form  of  C(<)  is  based  on  the  response  of  a  bank  of  pre-weighted  linear  compensators, 
shown  in  Figure  3.14.  Recall  that,  for  pre-wdghted  compensators,  once  a  compensator  in  the  bank 
is  activated  it  wiU  produce  an  output  for  all  time  forward.  The  output  of  (7t(<)  in  the  Figure  3.14 
proved  fully  capable  of  accounting  for  the  fact  that  compensator  G^it)  inherited  a  plant  out  of 
equilibrium.  No  oscillations  were  induced  as  long  as  there  were  no  modeling  inconsistencies.  With 
inconsistencies  present,  the  G3(/)  predicted  by  linear  analysis  was  incorrect.  In  our  case,  modeling 
errors  are  also  inherent  to  the  problem,  so  it  is  reasonable  that  the  plant  contnds  advocated  by  the 
linear-based  FHizzy  Bank  would  ultimately  induce  undesirable  oscillations.  The  question,  then,  is 
will  any  equivalent  Gi{t)  be  able  quell  the  oscillations  induced  the  Fuzzy  Bank. 
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Figure  5.8  Desired  form  for  C(t)  based  on  linear  analysis  of  pre-weighted  compensator  bank 

Tb  answer  this  question  a  fbxsy  Logic  element  was  created  which  would  emulate  the  output  oi 
Giit)  in  Figure  3.14,  but  with  arbitrary  tranuent  response  and  steady-state  output.  The  response 
of  G2{t)  in  Figure  3.14  is  repeated  in  Figure  5.8.  The  initial  spike  in  the  function  will  not  be  as 
important  in  this  application,  because  the  switches  used  in  Chapter  3  have  all  been  replaced  with 
fbzzy  Sets  for  the  MBFLC.  Ikansitions  between  compensators,  therefore,  occur  more  gradually. 
However,  an  initial  surge  will  be  required  to  halt  the  state  derivatives  reducing  to  a  steady-state 
offset  value  based  on  the  requirement  that  the  control  input  to  the  plant  too  he: 

t^(<«o)  =  i?r./(*oo)-n*0).  (5.1) 

For  an  input  Rr,j{t)  =  0.6  +  .l«_i(t),  starting  from  Y{t,)  =  0.6  then  V{t„)  =  0.1. 

The  function  shown  in  Figure  5.8  can  be  approximated  using  the  Fbzzy  Sets  EisO  and 
EisPositive.  The  membership  functions  for  these  two  sets  are  shown  in  Figure  5.9  based  on  the 
objective  underdamped  response.  Notice  that  £ts0  is  activated  initially  when  the  cloeed-lo(q>  re- 
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Figure  5.9  C(t)  function  formed  by  ANDing  of  FNizzy  Sets  EisO  and  EdotisPoaitive 

sponse  passes  through  zero.  EiaO  goes  to  zero  as  the  system  teaches  its  peak  overshoot  then,  is 
driven  to  1  again  as  the  system  settles  out  at  the  desired  reference  value. 

The  Rtzzy  Set  EisPoaitive  indicates  when  the  system  is  approaching  the  reference  value  from 
a  negative  value  of  Error.  This  can  be  seen  from  the  fact  that  E  =  (re/  -  T)  =  —Y.  When  Y 
initially  begins  climbing  Y  is  positive  so  E  is  negative.  E  will  only  become  positive  when  the  Y 
has  reached  its  peak  overshoot  and  started  back  down  towards  the  reference  value.  As  the  system 
settles  out  at  V  =  Ref,  the  Set  EiaPoaitive  reaches  a  steady-state  value  of  membership,  in  this 

C*®®  t*BitPo$itiva  —  0-35. 

The  desired  function  shape  is  obtained  using  the  FHizzy  AND,  discussed  in  Chapter  2.  The 
result  of  the  ANDing  procedure  is  also  shown  in  Figure  5.9.  The  resulting  function  reasonably 
approximates  the  function  given  in  Figure  5.8  for  an  arbitrary  starting  time. 

The  exact  form  of  the  Fuzzy-derived  C(t)  can  be  altered  by  the  tuning  parameters  of  the  Fuzzy 
Sets  and  a  constant  weighting  term  x?  called  the  Fuzzy  Offset.  The  procedure  for  determining  initial 
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estimates  of  these  parameters  is  ^ven  here.  Of  coarse,  simalation  is  required  to  maximise  the  C(t) 
damping  effect. 

The  first  step  is  to  determine  values  for  and 

These  quantities  completely  define  the  FHissy  Sets.  The  procedure  for  EisZero  is  straightfwward. 
C(t)  must  activate  as  the  system  output  nears  the  desired  state,  or  as  E  goes  to  0.  Therefme, 
•^m«s«r*=0.  determines  how  far  away  from  zero  C(t)  will  begin  to  activate.  Fbr  this 

design,  it  was  iteratively  decided  that  EisZero  should  be  activated  to  a  value  (tf  /i  =  0.2  when  the 
response  was  within  5 

The  purpose  of  the  EisPoaitive  Fuzzy  Set  is  two-fold:  To  prevent  C(l)  from  activating  when 
the  system  response  initially  crosses  £  =  0,  and  to  provide  steady-state  constant  into  the  plant.  It 
is  decided,  then,  that  fiBiiPosUivti^)  =  0.1  for  the  value  of  E  induced  by  the  IfKn(t)  when  E  first 
goes  to  zero  and  =  0.35.  The  values  for  a^BUPtiu*,  *>« 

solved  for  by  solving  the  system  of  two  equations  a  two  unknowns  using: 


H(x)  =  e-^^.  (5.2) 

The  value  of  E  induced  by  the  Uiin(t)  when  E  first  goes  to  zero  is  determined  by  simulation. 

X  can,  in  theory,  be  determined  by  the  constraint  given  in  (5.1).  For  for  the  case  shown  in 
Figure  5.9,  for  an  input  i2r«/(0  =  0.6  +  .lt»_i(<): 


y(U-o.i 

^  0.35 


(5.3) 


Unfortunately,  this  equation  most  be  solved  through  simulation,  because  U(t„)  cannot  be  deter¬ 
mined  analytically. 

The  form  of  the  Banked  Model  Based  FVizzy  Lo^c  Controller  incorporating  this  emulated 
(?3(t)  is  shown  in  Figure  5.10.  In  this  Figure,  the  Fbzzy  Bank  is  shown  as  a  single  block.  The 
structure  internal  to  the  block  FuzzyBank  is  identical  to  the  compensator  shown  in  Figure  5.5.  It 
was  determined  that,  through  proper  tuning,  the  oscillations  in  Y(t)  could,  indeed,  be  eliminated 
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Figure  5.10  SIMULINK  implementation  of  Banked  MBFLC 

for  a  given  initial  condition  and  step  magnitude.  This,  therefore,  is  a  successful  form  of  the  Modd- 
Based  I\izzy  Logic  Controller. 

The  membership  functions  for  the  Fuzzy  Sets  EisO  and  EisPositive  are  determined  through 
simulation  and  are  plotted  versus  their  universes  of  discourse  in  Figure  5.11.  The  tuning  data  for 
these  two  Fuzzy  Sets  is  given  in  Thble  5.3.  The  Fuzzy  tuning  parameters  were  chosen  so  as  to 
produce  the  correct  shape  of  the  C(t)  function.  The  magnitude  of  the  C{t)  is  determined  by  the 
Offset  Value  x-  X  is  slso  determined  through  simulation,  and  in  this  case,  x  =  0*13. 

The  response  of  the  final  Banked  MBFLC  is  shown  in  Figure  5.12.  Both  the  initial  transient 
and  settling  characteristics  of  this  system  are  similar  to  a  linear  response.  This  demonstrates  the 
applicability  of  a  hybrid  linear/FHizzy  approach  to  nonlinear  contrd. 
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Figure  5.11  LEFT:  Membership  function  for  E  is  Zero.  RIGHT:  Membership  function  for 
EisPositive 


Figure  5.12  Response  of  nonlinear  plant  in  closed-loop  simulation  with  Banked  MBFLC 
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5.3  Fidl- Envelope  Banked  MBFLC 


The  sources  of  error  identified  in  the  previous  Section  have  now  been  accounted  for  in  the 
MBFLC  design.  Errors  caused  by  intermediate  compensators  were  found  to  be  minimal.  Errors 
caused  by  input  error  signals  larger  than  the  regions  of  attraction  for  linear  compensators  are 
minimized  by  use  of  the  Fuzzy  Limiter.  Other  plant/compensator  mismatch  errors,  which  result 
result  in  lightly  damped  oscillations  near  the  desired  state,  are  addressed  with  an  additive  signal 
derived  using  Gaussian  Fuzzy  Sets. 

Consideration  of  these  error  sources  results  in  a  controller  capable  of  driving  a  nonlinear  plant 
from  a  given  initial  condition  {Y  =  0.6)  to  a  given  reference  value  {Y  =  0.7)  in  a  linear-like  manner. 
The  next  step  is  to  modify  the  compensator  so  that  the  linear-like  behavior  is  exhibited  regardless  of 
the  starting  point  and  ending  point.  A  compensator  which  yields  linear-like  performance  throughout 
a  specified  operating  range  will  be  referred  to  as  a  full-envelope  controller.  This  is  a  term  from  flight 
control,  which  implies  a  flight  controller  applicable  throughout  the  flight  envelope  of  an  aircraft. 

Consider  the  case  where  to  objective  is  to  develop  a  full-envelope  controller  for  the  envelope 
between  Y  =  0.6  and  Y  =  0.7,  using  equilibrium  models  at  Y  =0.6, 0.614, 0.628, 0.642, 0.656, 0.67, 
0.684,  0.698,  0.712  and  0.726.  The  form  of  the  Fuzzy  Bank  (Compensator  Bank,  Fuzzy  Supervisor 
and  Fuzzy  Limiter)  will  be  the  same  as  the  previous  section  (Figures  5.2,  5.4,  and  5.5).  The  Fuzzy 
Supervisor  will  select  as  the  initiating  compensator  that  which  most  closely  corresponds  to  the 
true  starting  equilibrium  of  the  nonlinear  plant.  It  has  been  shown,  however,  that  the  compensator 
based  on  the  final  desired  equilibrium  point  is  incapable  of  driving  the  system  to  steady-state  within 
the  desired  settling  time.  An  auxiliary  signal  C{t)  hsul  to  be  defined  to  quickly  drive  the  closed 
loop  system  to  steady-state. 

The  system  is  very  sensitive  to  the  form  of  the  C(t)  signal.  Small  errors  in  the  shape  of 
the  Fuzzy  Sets  EisZero  and  EisPositive  or  in  the  offset  term  x  degraded  system  response. 

Therefore,  like  the  Fuzzy  Limiter,  the  required  C(t)  function  for  various  step  magnitudes  within  the 
envelope  of  operation  must  be  developed  in  order  to  achieve  full-envelope  performance.  The  proper 
C{t)  for  intermediate  step  inputs  could  then  by  determined  through  interpolation.  The  form  of  the 
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Figure  5.13  SIMULINK  simulation  of  F\ill  Envelope  Banked  MBFLC 

controller,  then,  could  be  expressed  as  a  bank  of  linear  compensators  and  a  bank  of  C(t)-ptodudng 
Fuzzy  AND  elements,  with  a  second  Fuzzy  Supervisor  providing  the  interpolation. 

An  alternate  way  to  simulate  the  full-envelope  compensator  is  through  the  use  of  look-up  ta¬ 
bles.  Figure  5.13  shows  the  form  of  the  iHizzy  Controller  with  a  FHizzy  Offset  look-up  table  included. 
The  correct  value  for  x  is  obtained  based  on  the  magnitude  of  the  reference  step  {Ref  {in)  —  (/^{in)). 
Because  the  Fuzzy  Sets  X  is  0  and  Xdot  is  Positive  must  also  change  based  on  the  input,  step  mag¬ 
nitude  is  also  fed  into  these  blocks.  Look-up  tables  supply  the  values  of  i  and  a*  corresponding 
to  the  appropriate  C(f).  Two  other  changes  are  made  to  minimize  the  number  of  inputs  to  the 
compensator: 

•  The  derivative  block  is  brought  inside  the  compensator 

•  y  is  now  determined  by  V  =  Ref  —  E 

A  new  top-level  view  is  shown  in  Figure  5.14. 

Look-up  tables  are  needed  for  x  (offset),  •„•,*,  and  <r^Ei$P0$iUv  Again, 

this  information  also  could  have  been  expressed  by  Fuzzy  Sets,  using  a  bank  of  C{t)  generators 
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Figure  5.14  Top-level  SIMULINK  simulation  of  nonlinear  plant  in  Full-Envelope  Banked  MBFLC 

and  a  Fuzzy  Supervisor.  The  look-up  tables  serve  to  simplify  the  simulation.  Look-up  tables  are 
limited  in  that  they  can  only  perform  linear  interpolation,  while  iHizzy  Supervisor  can  be  tuned 
to  enhance  performance.  In  this  case,  linear  interpolation  is  not  a  significant  drawback.  The  data 
contained  in  the  look-up  tables  is  given  in  Table  5.4. 

The  response  of  this  system  to  step  inputs  of  various  magnitudes  is  shown  in  Figure  5.15.  The 
principle  difference  between  the  response  of  this  closed-loop  system  and  that  of  the  model  M(s) 
2.26  is  that  the  system  peaks  slightly  later  in  time.  M(8)  is  fully  settled  out  by  t  =  1.62sec(mds 
while  the  Full  Envelope  MBFLC  is  just  reaching  the  reference  value  at  this  time.  Some  of  this 
error  is  attributed  to  the  integration  routine  used  by  SIMULINK  to  solve  nonlinear  equations.  It 
is  possible,  however,  that  certain  plant  characteristics  cause  the  Fuzzy  Bank  to  exhibit  trajectory 
errors. 

Figure  5.16  shows  the  closed-loop  response  of  the  nonlinear  plant  for  various  values  of  the 
parameter  a.  The  Figure  shows  that  the  MBFLC  is  more  sensitive  to  modeling  errors  than  either 
QFT  or  FLC.  This  is  due  to  the  tuning  of  the  Fuzzy  Limiter  which  assumes  a  given  plant  sensitivity. 
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Table  5.4  Look-up  table  data  for  full  envelope  MBFLC 


The  effect  of  modeling  errors,  however,  is  not  catastrophic.  The  effect  of  retuning  the  Limiter  on 
modeling  errors  is  explored  in  the  next  section. 
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Figure  5.15  Cloeed-loop  response  of  FNdl-Envelope  MBFLC  to  step  inputs  of  various  magnitudes 


Figure  5.16  Closed-loop  response  of  I^ill-Envelope  MBFLC  for  a=  0.5, 1.0,  and  1.5 


5.4  Contrilnttion  of  Linear  Compensator*  to  Banked  MBFLC  Design 

The  previous  section  developed  a  Model-Based  Fuzzy  Lofpc  Coatroller(MBFLC)  that  meets 
all  established  design  specifications,  with  the  exception  of  disturbance  rejection.  In  essence,  the 
MBFLC  is  a  confederation  of  n  controllers  consisting  of  an  initial  linear  compensator  and  a  C{t) 
generator.  The  F\izzy  Limiter  is  common  to  all  the  compensator/generator  pairs.  Fuzzy  Logic, 
however,  effects  both  the  input  to  and  the  output  from  the  bank  of  linear  compensators.  Therefore, 
it  is  reasonable  to  question  just  how  critical  a  bank  of  correctly  designed  compensators  is  to  the 
overall  success  of  the  design. 

Several  tests  can  be  run  on  the  MBFLC  to  determine  the  contribution  made  by  the  linear 
elements  of  this  hybrid  controUer.  One  is  to  determine  the  performance  of  the  Banked  MBFLC  as  a 
function  of  the  number  of  compensators  in  the  bank.  The  bank  currently  uses  enough  compensators 
so  that  at  least  one  is  considered  valid  at  any  point  in  the  state  trajectory.  If  full  coverage  is,  indeed, 
a  requirement  then  reducing  the  number  of  linear  compensators  should  be  reflected  by  a  degradation 
of  system  performance.  A  second  test  would  be  to  substitute  the  nonrobust  compensators  in  the 
current  bank  with  the  robust  compensators  devdoped  in  Chapter  3.  If  robustness  in  the  linear 
sense  translates  directly  into  robustness  in  the  nonlinear  sense,  this  would  suggest  a  significant 
contribution  being  made  based  on  linear  theory.  A  third  test  would  be  to  alter  the  a  parameter  of 
the  nonlinear  plant  to  determine  if  mismodeling  can  be  overcome  through  tuning.  These  questions 
are  explored  here. 

5.4.1  Effect  of  Linearized  Plants  on  MBFLC.  Knowledge  of  the  nonlinear  system  is 
resident  in  two  components  of  the  controller:  the  linearized  plants  and  the  FHizzy  Sets.  These 
data  sources,  however,  are  interrelated.  The  Fuzzy  Sets  are  tuned  based  on  simulations,  and  the 
simulations  aire  performed  with  a  given  number  of  linear  compensators  in  place.  It  is  reasonable, 
then,  to  assume  that  there  is  some  redundancy  in  the  compensator  architecture,  and  that  the 
structure  can  be  ultimately  simplified  with  little  penalty  in  performance.  Based  on  the  analysis 
of  the  previous  sections.  It  is  assumed  that  only  the  initiating  and  terminating  compensators 
contribute  significantly  to  the  success  or  failure  of  the  compensator.  In  addition,  it  has  been  shown 
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Figiire  5.17  SIMULINK  sinmlation  of  MBFLC  luing  only  a  single  linear  contndltt 

that  the  linear  compensator  based  on  the  desired  ontpnt  of  the  system  is  unable  to  drive  the 
system  to  the  reference  value  with  zero  steady-state  error.  This  indicates  that  only  the  initiating 
compensator  could  be  considered  a  linear  element  of  the  MBFLC. 

Consider  the  simulation  shown  in  Figure  5.17.  Here  the  bank  has  been  replaced  by  a  an^ 
compensator,  with  r  corresponding  to  F  =  0.6.  This  results  greatly  simplifies  the  simnlati<m, 
allowing  it  to  be  displayed  on  a  single  screen.  The  F\izzy  Supervisor  can  be  removed,  as  there 
are  no  alternate  compensators  to  weight.  The  closed-loop  response  oi  this  qrstem  to  step  inputs 
y  =  0.6  -H  0.025tt-i(f),  Y  =  0.6  0.05«-t(0>  <^<1  F  =  0.6  +  0.1v_i(t)  is  shown  in  Figure  5.18.  As 
anticipated,  the  response  is  not  adversely  affected  by  the  deletion  of  the  other  compensators  in  tlm 
bank.  In  fact,  the  response  has  improved  in  some  respects,  such  as  settling  time.  What  is  surpriang, 
however,  is  that  tlus  increase  in  performance  is  independent  of  the  starting  equilibrium  condition. 
Figure  5.19  shows  closed-loop  responses  for  Y  =  0.65  -f-  0.025«_i(<),  Y  =  0.65  -f  0.05«_i((),  and 
y  =  0.62  +  0.03«_i(t).  All  of  these  exhibit  acceptable  closed-loop  responses. 
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Figure  5.18  Closed-loop  response  of  single-compensator  MBFLC  to  step  inputs  d  varions 
magnitudes 

This  independence  of  the  initial  state  is  not  indq)endent  of  the  choice  of  compensator.  In 
fact,  the  larger  the  value  of  r  (that  is,  the  smaller  the  gain  in  the  single  compensator)  the  more 
sluggishly  the  system  responds.  This  h<dds  for  steps  both  above  and  below  the  equilibrium  point 
of  the  compensator.  This  is  illustrated  in  Figure  5.20.  The  single  compensator  in  Figure  5.17 
is  altered  to  correspond  to  F  =  0.65  (r  =  1.2675).  Figure  5.20  shows  the  closed-loop  response 
of  the  system,  starting  from  equilibrium  both  above  and  below  the  equilibrium  about  which  the 
compensator  is  designed.  In  both  cases  performance  is  slng^sh. 

The  performance  of  the  Single-Compensator  MBFLC,  then  is  a  function  of  compensator 
selection.  The  Fuzzy-emnlated  G3(t)  will  correctly  stabilize  the  system  response  provided  the 
linear  compensator  drives  the  system  to  the  correct  overshoot.  The  Fuzzy  Limiter,  however,  will 
overdamp  all  compensators  but  the  one  with  the  highest  gain  in  the  ori^nal  Limiter  design.  The 
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Figure  5.19  Qaeed-loop  respoase  of  single-compeiuator  MBFLC  to  step  inputs  stnrting  from 
various  initial  equilibria 

conclusion,  therefore,  is  that  the  single  compensator  should  be  the  one  with  the  highest  gain  in  the 
envelope  of  operation.  For  this  design,  the  ’’correct”  single  compensator  corresponds  to  F  =  0.6. 

It  has  now  been  established  that  a  single-compensator  MBFLC  can  be  developed  ^ven  the 
tuned  elements  of  the  full  10-compensator  MBFLC.  It  is  reasonable  to  ask  if  the  same  tuning 
data  could  have  been  obtained  using  the  single  compensator  design.  The  answer  is  yes.  The 
tuning  procedure  can  be  accomplished  even  with  only  a  single  compensator  providing  trajectory 
information.  An  arbitrary  number  of  Fbzzy  Sets  can  be  defined  within  the  Fuzzy  Limiter  and 
the  m<^nitnde  of  the  required  E  penalty  can  be  determined  using  the  progressive  step  procedure 
discussed  in  the  last  section.  The  last  Fuzzy  Set  should  correspond  to  the  largest  expected  input  to 
the  system.  The  look-up  tables  for  x  (offset),  <r^6up„uiw,fXBup,tiu„t  *“<1  must  then 

be  tuned  through  simulation  to  minimize  the  oscillatory  tendencies  of  the  system. 
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Figure  5.20  Closed-loop  response  of  single-compensator  MBFLC  to  step  inputs  starting  from 
equilibria  both  above  and  below  the  equilibrium  of  the  linear  compensator 

The  final  issue  to  resolve  is  the  minimum  number  of  Fuzzy  Sets  in  the  Limiter  required  to 
provide  full-envelope  control  of  the  nonlinear  plant.  Figure  5.21  shows  the  MBFLC  simulation  with 
the  Fuzzy  Limiter  reduced  to  two  Sets.  Because  of  the  greatly  reduced  size  of  the  Fuzzy  Limiter,  its 
special  block  was  removed  so  the  entire  simulation  could  be  viewed  in  a  single  diagram.  The  F\izzy 
Limiter  and  Linear  compensator  could  be  adjusted  to  ensure  the  correct  overshoot.  However,  to 
accomplish  this,  the  requirement  that  the  £rror«s0  Fuzzy  Set  in  the  Limiter  carry  a  penalty  of  1 
by  definition  had  to  be  eased.  By  setting  the  penalty  associated  with  this  Fhzzy  Set  to  0.8  and 
tuning  the  look-up  tables  to  damp  oscillations,  the  response  in  Figure  5.22  is  obtained. 

Figure  5.22  also  shows  the  response  of  the  compensator  to  steps  of  other  magnitudes  starting 
from  the  same  equilibrium  condition.  This  shows  that  the  two-set  Limiter  is  not  adequate  to  contitd 
the  full  envelope  of  this  system. 

The  conclusions  of  this  analysis  are  as  follows: 
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•  Only  one  Uneax  compensator  is  required  in  a  MBFLC.  This  compensatcw  should  correspond 
to  that  portion  of  the  plant  envelope  requiring  the  highest  gain. 

•  In  general,  one  Fuzzy  Set  will  be  required  in  the  Fuzzy  Limiter  for  each  desired  final  state. 
However,  a  break-even  point  occurs  where  the  entire  operating  envelope  is  ”adeqnatdy  cov¬ 
ered”  and  a  step  input  to  any  point  the  envelope  can  be  made.  The  number  of  sets  required 
to  reach  this  point  seems  to  vary  according  to  the  extent  of  the  nonlinearity. 

•  Increasing  the  number  of  linear  compensators  does  not  overcome  a  defidency  in  the  granularity 
of  the  Fuzzy  Limiter. 

•  Removing  all  linear  compensators  or  repladng  the  single  compensator  with  a  simple  gain 
results  in  poor  performance. 

•  The  Fuzzy  elements  play  a  significant  role  fordng  the  nonlinear  system  to  exhilut  linear 
behavior. 
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Figure  5.21  Single-compeneator  MBFLC  with  a  simplified  2-8et  Frizzy  Limiter 
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Figure  5.22  SIMULINK  simulation  of  MBFLC  with  a  two-set  Fuzzy  Limiter  and  a  sin^e  linear 
compensator 
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l^ble  5.5  Tuning  parameters  for  FHizzy  Limiter  driving  robust  linear  compensator 
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T^ble  5.6  Look-up  table  data  for  MBFLC  incorporating  the  robust  linear  compensator 

5.4.  i  Effect  of  Robust  Linear  Compensators  on  MBFLC.  A  second  test  of  the  contribution 
being  made  by  linear  contrcd  theory  is  to  replace  the  linear  compensator  used  in  the  last  section 
with  the  robust  compensator  developed  in  Chapter  3.  If  robustness  in  the  Unear  sense  translates 
via  the  MBFLC  architecture  into  nonUnear  robustness,  this  indicates  a  direct  appUcation  of  linear 
synthesis  to  nonUnear  control  through  f\izzy  Logic. 

Figure  5.23  shows  the  single  compensator  MBFLC  with  the  robust  controUer  inserted.  Five 
equiUbrium  conditions  were  chosen  for  tuning:  =  0.6,0.625,0.65.0.675,  and  0.7.  hNizzy  Sets 

were  defined  within  the  f\izzy  Limiter  corresponding  to  each  of  these  equiUbrium  conditions  and 
the  MBFLC  was  tuned  as  discussed  in  the  last  section.  The  fNizzy  Limiter  tuning  data  is  ^ven 
in  Table  5.5.  The  tuning  data  for  the  look-up  tables  x>  variance,  EisPositve  variance,  and 
EisPoaitve  mean  are  given  in  Table  5.6. 

Notice  the  data  in  the  look-up  tables  are  much  difierent  from  the  case  of  the  nonrobust 
controller.  TI  ough  the  values  for  the  E  and  E  Fuzzy  Sets  are  stiU  nearly  the  same,  the  offset  value 
X  is  significantly  larger.  The  penalties  associated  with  the  Error  signal  in  the  Fuzzy  Limiter  are 


5-33 


J 


Figure  5.23  SIMULINK  simulation  of  single  compensator  MBFLC  using  robust  compensate 

also  an  order  of  magnitude  higher  in  the  nonrobust  case.  These  differences  are  necessary  to  account 
for  the  much  higher  gain  of  the  robust  compensate. 

The  response  of  this  system  to  step  inputs  of  various  amplitudes  is  shown  in  Figure  5.24. 
Performance  is  equivalent  to  the  nonrobust  case  of  the  previous  section.  Figure  5.25,  however,  shows 
the  responses  of  both  the  robust  and  nonrobust  single  compensator  MBFLCs  to  a  step  disturbance 
of  magnitude  0.1  at  the  input  of  the  plant.  This  Figure  shows  that  the  nonrobust  compensator 
actually  exhibits  better  signal  rejection  characteristics  than  does  the  robust  compensator. 

The  disparity  in  performance  is  due  to  the  F\izzy  Limiter.  With  the  FHizzy  Limiter  removed, 
the  single  robust  compensator  MBFLC  meets  the  \YpisTVRBANCB{t)\  ^  |0.1I/i>/STnitBiiwcn(Oi 
specification,  but  cannot  track  a  signal.  With  the  FHizzy  Limiter  inserted,  the  system  can  track 
a  step  input,  by  the  penalties  associated  with  a  large  E  signal  inhibit  signal  rejection.  The  other 
brsmch  of  the  MBFLC  activated  by  EiaO  and  EisPositive  does  not  contribute  to  the  regulation 
problem. 
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Figure  5.24  Closed-loop  tracking  response  of  robust  Single  Compensator  MBFLC  to  step  inputs 
of  various  magnitudes 

For  the  case  of  the  nonrobust  compensator,  a  balance  is  struck.  The  error  gnal  is  not 
attenuated  nearly  as  fast  as  in  the  robust  case.  Therefore,  as  the  E  gets  larger,  the  compensator 
is  still  able  to  drive  the  plant  back  to  the  reference  value  in  a  timely  manner.  Even  in  this  case, 
the  Fuzzy  Limiter  serves  to  inhibit  adequate  tracking.  The  effect  of  the  Limiter,  though,  is  not 
nearly  as  severe.  Notice  that,  for  regulation  with  the  Limiter  removed,  disturbance  response  is 
approximately  that  which  was  determined  through  linear  analysis  in  Chapter  3.  This  equivalence 
of  linear  and  nonlinear  controllers  for  regulation  could  be  an  area  for  future  exploration. 

An  additional  simulation  was  done  to  determine  if  reintegration  of  the  full  compensator  bank 
would  have  a  significant  effect  on  signal  rejection.  The  conclusion  was  that  it  did  not.  Therefore,  this 
analysis  suggests  that  linear  robustr'ess  does  not  contribute  directly  to  robustness  in  the  MBFLC. 
The  robustness  provided  by  the  controller  in  the  linear  case  will  be  observed  in  the  nonlinear  case 
in  inverse  proportion  to  the  weights  which  determine  the  squashing  effect  of  the  Fuzzy  Limiter.  The 
smaller  the  slope  of  this  line,  the  greater  the  re^on  of  attraction  of  the  linearized  compensator, 
and  the  more  robustness  will  be  be  observed  in  the  nonlinear  system. 
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Figure  5.25  Disturbance  rejection  response  of  both  robust  and  nonrobust  MBFLC  to  a  step  dis¬ 
turbance  of  magnitude  0.1  injected  at  the  plant  input.  Also  shown  is  the  response  oi 
the  nonrobust  compensator  with  the  Fuzzy  Limiter  removed. 
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5.4.S  Effect  of  Mismodeling  on  Compensator  Tuning.  For  this  experiment,  the  value  of 
the  a  parameter  within  the  nonlinear  plant  was  altered  to  a  =  1.5.  With  the  MBFLC  tuned  as 
in  the  other  two  experiments,  the  response  is  extremely  underdamped.  The  tuning  procedure  was 
repeated,  both  for  the  Fuzzy  Limiter  and  the  C{t)  functions.  Upon  retuning,  the  closed-loop  system 
could  be  driven  within  the  normal  design  specifications. 

Changing  the  value  of  a  alters  the  equilibrium  condition  of  the  nonlinear  plant.  The  starting 
conditions  on  Xi  and  Xj  were  altered  to  account  for  this.  For  a  =  1.5,  A'i(t0)  =  1.5X3(19).  The 
tuning  parameters  for  the  case  of  a  mismodeled  plant /compensator  are  given  in  Tables  5.7  and  5.8. 
Note  that  the  effect  of  the  mismodeling  was  to  force  the  reduction  of  the  penalty  in  the  Limiter 
associated  with  EisO.O  from  1  to  0.8.  If  this  penalty  had  not  been  reduced  the  values  of  the  penalties 
would  not  have  gotten  steadily  smaller  as  the  error  grew. 

This  experiment  suggests  that  the  design  of  the  linear  compensator  is  not  critical  to  proper 
performance.  There  is  some  design  latitude  as  long  as  the  compensator/Fuzzy  Limiter  pair  induce 
a  smooth  transition  from  the  initial  state. 

This  experiment  also  shows  that  the  MBFLC  architecture  is  potentially  very  robust,  given  a 
mechanism  to  adjust  the  weights  in  the  Fuzzy  Limiter,  based  on  perceived  performance.  Thble  5.8 
shows  that  the  values  determining  C{t)  changes  very  little,  as  long  as  the  Fuzzy  Limiter  and  linear 
compensator  induce  the  correct  overshoot. 
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Table  5.7  Inning  parameters  for  F\izzy  Limiter  driving  the  mismodeled  nonlinear  plant  (a=1.5) 
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Table  5.8  Look-up  table  data  for  MBFLC  incorporating  the  robust  linear  compensator 


5-38 


5.5  ModeUFoUomng  Hybrid  CompensaU>r$ 

The  previous  section  showed  that  there  was  very  little  contribution  to  the  success  of  the 
MBFLC  made  by  its  linear  elements.  Though  the  MBFLC  architecture  is  based  largdy  on  linear 
design  considerations,  the  controller  nltimatdy  relies  on  Fiizzy  Sets,  tuned  through  simulation, 
to  achieve  successful  full-envelope  performance.  Linear  compensators,  it  seems,  have  significant 
limitations,  even  when  employed  in  conjunction  with  FNizzy  Logic. 

Linearized  plant  models  are  not  the  only  way  in  which  linear  theory  can  contribute  to  the 
control  of  a  nonlinear  plant.  An  alternate  approach  is  to  consider  an  actual  model-following.de8ign, 
similar  to  the  nonlinear/Fuzzy  hybrid  controller  from  Chapter  4.  In  this  design,  the  linear  system 
model  developed  in  Chapter  2  is  integrated  dir«:tly  into  the  controller  structure.  FHizzy  Lo^c  is 
then  used  to  drive  the  error  between  the  plant  states  and  model  states  to  zero  at  all  times. 

Linear/Fuzzy  model  following  is  illustrated  in  Figure  5.26.  For  this  controller,  two  signals  are 
fed  back:  Y{Xi)  and  Xt.  Y  is  used  to  drive  the  linear  system  model  contained  in  the  SIMULINK 
block  Model  and  described  by  (2.26).  The  internal  structure  for  the  Model  block  is  shown  in 
Figure  5.27.  This  block  is  essentially  a  closed-loop  simulation  of  the  model  poles  with  the  reference 

signal  and  the  feedback  term  supplied  externally. 

•  • 

Based  on  the  model  response,  a  desired  Xi  trajectory  is  generated.  The  actual  Xi  is  sub¬ 
tracted  from  the  desired  value  and  the  difference  (trajectory  error)  is  input  to  the  block  called 
FStzzyDriver.  The  purpose  of  the  Fuzzy  Driver  is  to  produce  an  error  signal  which  will  drive  the 
error  between  Xi-desired  and  JITi-actual  to  zero.  The  error  signal  is  fed  into  a  linear  compensator, 
and  the  compensator  ultimately  drives  the  plant. 

The  internal  structure  of  the  Fuzzy  Driver  is  shown  in  Figure  5.28.  Within  the  Fazzy  Driver, 
the  Xi  Error  signal  is  evaluated  for  membership  in  three  different  F\iz^  Sets:  Negative,  Zero, 
and  Positive.  The  membership  functions  for  these  sets  are  shown  in  Figure  5.29.  The  membership 
functions  for  Negative  and  Positive  are  open-ended.  That  is,  they  approach  1  as  the  Xx  Error  signal 
approaches  either  positive  or  negative  oo. 
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Figure  5.26  Hybrid  linear/Fuzzy  Controller  using  linear  model-following  and  Linear/FNizzy  contnd 
input  generation 

Each  Ftazzy  Set  has  associated  with  it  an  error  signal,  very  lar^  rdative  to  the  contnd  inputs 
typically  required  to  drive  the  plant.  Small  deviations  in  Xi,  then,  result  in  a  large  correcting 
error  signal  to  the  linear  compensator,  shown  in  Figure  5.26.  The  linear  compensator  drives  the 
nonlinear  plant.  Error  signals  of  ±10  were  chosen  arbitrarily  to  provide  excellent  model-following 
and  signal  rejection  regardless  of  the  value  of  the  parameter  a.  Large  error  signals  associated  with 
small  error  signals  is  typical  of  a  bang/bang  type  controller  [21]. 

The  linear  compensator  is  required  to  hold  the  system  at  steady-state  once  all  dynamics 
have  subsided.  Because  there  is  no  error  signal  associated  with  an  Xi  Error  of  Zero  within  the 
Fuzzy  Driver,  no  control  will  be  generated  when  the  actual  and  desired  trajectories  coincide.  The 
integrators  within  the  linear  compensator  are  needed  to  obtain  a  constant  input  to  the  plant  after 
all  dynamics  have  subsided. 

The  response  of  this  system  to  step  inputs  of  various  magnitudes  is  shown  in  Figure  5.30.  All 
restrictions  on  the  re^on  of  operation  can  be  removed,  resulting  in  a  true  full-envelope  controller. 


Yfm) 

Figure  5.27  Internal  structure  of  SIMULINK  block  Model 


This  strategy  only  breaks  down  for  extremely  large  step  inpnts  which  require  larg^  contnd  inputs 
than  are  provided  for  by  the  Implicates  of  the  Fbzzy  Sets.  Variations  in  the  a  parameter  produces 
no  change  in  the  system  response.  Similarly,  this  controller  demonstrates  almost  complete  rejection 
of  unwanted  signals  at  the  input  to  the  nonlinear  plant. 

The  price  paid  for  this  performance  is  in  extreme  swings  at  the  plant  input.  Figure  5.31 
shows  the  input  to  the  nonlinear  plant  alter  the  fbzzy  Driver  and  the  linear  compensator.  The 
Figure  shows  a  great  number  of  input  spikes  which  would  cause  havoc  in  any  physical  system.  Of 
course,  the  input  spikes  are  exaggerated  by  the  choice  ±10  for  the  implicate  Error  signals  within 
the  Fuzzy  Driver.  For  the  ori^nal  r^on  of  operation,  between  Y=0.6  and  Y=0.7,  implicate  error 
signals  of  ±0.8  are  adequate  to  drive  the  plant  to  within  a  reasonable  approximation  of  the  model. 
Input  spikes  are  still  present  in  the  plant  input,  but  are  not  nearly  a  severe  as  those  shown  in  the 
Figure. 

The  control  spike  problem  can  be  mitigated  by  reinserting  the  FWzy  Limiter  and  inputing  the 
sum  of  the  Driver  output  and  the  Limiter  output  to  the  linear  compensator.  In  this  configuration. 
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the  Fuzzy  Driver  will  only  need  to  generate  a  driving  signal  when  the  plant  model  deviates  firom  the 
model  trajectory.  The  errors  seen  by  the  Driver  will  generally  be  much  smaller,  and  the  generated 
driving  signal  commensurately  smoother. 

This  controller  design,  unfortunately,  invokes  the  state  directionality  argument  that  the 
Model-Based  Fuzzy  Lope  Control  approach  was  seeking  to  avoid.  It  is  assumed  that  the  cor¬ 
rect  sign  (positive  or  negative)  of  the  input  can  be  inferred  by  the  sign  of  the  error  signal.  This 
may  not  always  be  the  case  for  nonlinear  plants  [16]. 

5.6  Summary 

A  successful  MBFLC  can  be  developed  by  accounting  for  the  three  error  sources  identified  in 
Chapters  3  and  4  separately.  The  resulting  design  is  capable  of  achieving  linear-like  behavior  from 
the  nonlinear  plant,  ^ven  fixed  initial  conditions  and  final  states.  With  the  insertion  of  look-up 
tables  or  an  additional  Fuzzy  Bank  of  C{t)  generators,  lull-envelope  controller  performance  can  be 
obtained. 
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Figure  5.29  Membership  functions  f<»^  Fuszy  Sets  within  the  Fuzzy  Driver 

Initial  analysis  indicates  that  Fuzzy  Lo^c  plays  the  significant  nde  in  the  adequate  compen¬ 
sation  of  the  nonlinear  plant.  Only  one  linearized  compensator  was  found  to  be  necessary,  and  its 
correct  design  was  not  critical.  The  addition  of  more  compensators  was  actually  detrimental  to  the 
closed-loop  response  of  the  nonlinear  plant.  Further,  mismodeling  between  the  single  compensator 
and  the  plant  can  be  easily  compensated  for  in  the  Fuzzy  Limiter. 

Because  the  full  bank  of  compensators  is  not  necessary,  this  design  clearly  does  not  emulate  the 
fuzzy /nonlinear  controller  developed  at  the  end  of  the  last  chapter.  In  the  FVizzy/Nonlinear  design 
of  Chapter  4,  the  contribution  of  the  linear  compensators  was  understood  in  the  context  of  the  Time- 
Varying  Linear  Plant.  The  time-varying  linear  compensator  and  it  time  invariance  sq>proximation, 
the  Fuzzy  Bank,  were  able  to  successfully  contrd  the  LTV  plant,  which  exhibits  a  subset  of  the 
nonlinear  dynamics  of  the  fuU  nonlinear  plant  under  consideration.  The  IWzy  Wdghted  Bank 
was  shown  to  closely  approximate  the  time-varying  linear  compensator.  The  current  F\izzy/Linear 
design  was  shown  to  approximate  two  pre-weighted  linear  compensators,  with  the  output  of  the 
second  linear  compensator  being  emulated  via  Fuzzy  Logic. 
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Figure  5.30  Closed-loop  responses  of  Model-Following  hybrid  controller  for  step  inpnts  of  vairions 
magnitudes  and  initial  conditions. 


The  single-compensator  MBFLC  was  found  to  possess  some  disturbance  rejection  capabilities. 
The  F\izzy  Linuter,  however,  actually  inhibits  robustness  in  inverse  proportion  to  the  magnitude 
of  the  penalties  on  E.  This  effect  becomes  significant  when  tuning  the  compensator  for  the  robust 
linear  compensator.  Any  advantages  to  using  the  robust  compensator  are  completely  lost  due  to 
the  weights  required  in  the  Fuzzy  Limiter  to  ensure  tracking  performance. 

An  actual  model-following  hybrid  linear/Fuzzy  controller  was  shown  to  meet  all  performance 
and  signal  rejection  requirements,  though  at  the  expense  of  signal  spikes  at  the  plant  input.  A 
balance  between  desired  performance  and  maximum  allowable  input  variation  could  result  in  a 
successful  hybrid  controller  design  for  real-world  applications. 
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Figure  5.31  Input  signal  to  nonlinear  plant  generated  by  F^i2zy  Driver  and  linear  compensatOT 
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VI.  Conclusions  and  RecommendatUms 


6.1  Introduction 

Moet  systems  of  interest  in  engineering  are  nonlinear,  at  least  at  the  extremes  ol  their  at- 
velopes  of  operation.  Controllers  for  these  systems  are  often  based  on  a  linearized  modd  oi  the 
plant,  or  on  multiple  linearizations  when  the  plant  exhibits  significant  nonlinearities.  This  practice 
leads  to  significant  modeling  error  between  the  physical  plant  and  the  plant  modd. 

The  objective  of  this  research  effort  was  to  devdop  a  hybrid  contrd  ^proach  to  account 
for  the  plant /linear  modd  ambiguity  in  a  meaningful  but  non-statistical  way.  The  non-statistical 
nature  of  the  approach  was  important  due  to  the  difficultly  in  accuratdy  ganging  the  extent  of  the 
modeling  error.  Other  error  sources  such  as  time-varying  parameters  can  also  manifest  themsdves 
in  a  manner  similar  to  modeling  error  and  further  motivates  a  non-statistical  approach.  FVizzy 
Lojpc  concepts  were  chosen  as  a  natural  framework  within  which  to  address  these  errors. 

The  devdoped  control  approach,  dubbed  Modd-Based  Fuzzy  Logic  Control  (MBFLC),  is  a 
logical  extension  of  both  linear  and  Fuzzy  Lc^c  Control  ideas.  In  essence,  the  amtrdler  is  FHizzy 
but  with  a  resident  system  ’’expert,”  represented  by  a  linear  compensaUv.  The  MBFLC  approach 
is  designed  to  take  advantage  of  the  strengths  of  both  Unear  systems  theory  and  FVizzy  Set  theory. 
Linear  dements  are  to  contribute  a  degree  of  a  priori  synthesis  and  analysis  capabiUty  and  the 
abiUty  to  utiUze  available  but  inaccurate/incomplete  system  models,  frizzy  Lo^c  dements  are  to 
contribute  robustness  in  the  face  of  modeUng  (and  Unearization)  errors,  enabUng  the  appUcation 
of  the  resulting  compensator  to  uncertain  Unear  and/or  nonUnear  problems.  The  final  form  of 
the  controUer  was  devdoped  in  the  last  chapter,  and  tests  were  performed  to  assess  the  rdative 
contributions  of  the  Fuzzy  and  Unear  dements. 

The  degree  to  which  the  desired  performance  characteristics  were  obtained  in  simulation  is 
the  subject  of  this  conduding  chapter.  The  succ^  of  MBFLC  is  assessed  in  Ught  of  the  estabUshed 
research  objectives. 

The  chapter  is  divided  into  three  sections: 

1.  Controller  Summary 
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Figure  6.1  MBFLC  Block  Diagram 


2.  Thesis  Conclusions 

3.  Recommendations  for  Future  Research 

4.  Summary 

6.2  Controller  Summary 

The  objective  of  this  thesis  was  to  develop  a  foil-envelope  contrdler  for  a  nonlinear  plant 
using  both  linear  control  systems  theory  and  F\iKzy  Lofpc.  Based  on  the  linear  and  nonlinear 
design  considerations  discussed  in  previous  chapters,  the  MBFLC  of  Figure  6.1  was  obtained. 

The  MBFLC  is  essentially  a  Fuzzy  Controller  with  only  two  rules: 

1.  If  Error  is  not  Zero  AND  Error  is  not  Positive  THEN  U  =  Unnit) 

2.  If  Error  is  Zero  AND  Error  is  Positive  THEN  V  =  Viin{t)  +  (7(1) 

That  is  to  say,  as  long  as  the  system  is  far  from  the  reference  value  and  still  climbing  {Error 
is  negative),  then  the  control  input  to  the  plant  is  Uti„{t).  When  close  to  the  reference  value  and 
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dropping  back  towards  the  reference  value  after  the  overshoot,  the  input  to  the  plant  is  C^Na(0+C(t). 
Neither  (/«„  nor  C(t)  are  directly  derived  from  linear  contn^  theory,  though  both  are  based  on  linear 
concepts.  These  signals  are  discussed  in  more  detail  below. 

6.2.1  Formation  of  Uu^{t).  As  shown  in  the  Figure,  ffu.  is  formed  in  two  stages.  First, 
the  Error  signal  is  fed  through  a  "squashing  function"  implemented  in  Fuzzy  Logic.  This  Frizzy 
Limiter  allows  small  error  signals  to  pass  through  unaffected,  but  rapidly  attenuates  larger  errors. 
The  Limiter  is  necessary  to  acknowledge  the  fact  that  the  linear  compensator  is  not  valid  over  the 
entire  range  from  Yo  to  Yr,J.  The  exact  form  of  the  FHizzy  Limiter  is  determined  through  simulation. 

The  modified  Error  signal  is  then  input  to  a  linear  compensator.  The  compensator  is  based  on 
a  linearized  model  of  the  nonlinear  plant  and  is  designed  to  the  performance  specifications  desired 
for  the  overall  simulation.  Be  .use  the  Limiter  is  designed  to  only  constrict  the  Error  signal  and  not 
to  amplify  it,  the  compensator  used  should  correspond  to  the  most  "sluggish”  part  of  the  nonlinear 
plant  envelope.  In  other  words,  the  compensator  with  the  highest  DC  gain  should  be  selected  for 
the  MBFLC. 

It  was  shown  through  simulation  that,  by  adjusting  the  Error  Magnitude  penalties  associated 
with  the  Fuzzy  Limiter,  a  closed-loop  response  with  the  correct  overshoot  and  initial  n.,  settling 
time  can  be  obtained  from  the  (/{,-„  signal  alone.  With  no  other  compensation,  however,  large, 
lightly- damped  oscillations  will  be  induced  in  the  nonlinear  plant.  It  is  the  purpose  of  the  added 
signal  C(t)  to  damp  out  these  oscillations. 

6.2.2  Formation  of  C{t).  The  correct  form  for  the  damping  signal  C{t)  is  also  drawn, 
in  a  non-rigorous  way,  from  linear  analysis.  Simulations  were  performed  to  test  the  feasibility  of 
driving  a  linear  plant  using  a  bank  of  two  identical  linear  compensators  in  closed-loop  simulation. 
It  was  found  that  control  of  the  linear  plant  could  be  transferred  from  compensator  (7|(s)  to  ^^(s) 
at  any  point  along  the  state  trajectory  provided  that  the  input  of  (?i(s)  be  driven  to  zero  at  the 
time  of  the  transfer  and  that  the  output  of  (7i(s)  be  allowed  to  continue  to  influence  the  plant  for  all 
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time  forward.  The  residual  signal  produced  by  <7i(a)  successfully  accounted  for  the  non-equilibrium 
condition  of  the  plant  at  the  time  of  the  compensator  transfer. 

This  effect  of  two  compensators  perfectly  complementing  each  other  was  only  observed  for 
identical  linear  compensators  driving  a  LTI  plant.  When  there  was  a  mismatch  between  the  two 
compensators  or  the  r  mpensators  and  the  plant,  lightly  damped  oscillations  resulted.  Linear 
analysis,  then,  infers  that  there  may  exist  a  combination  of  the  signals,  similar  to  the  linear/no- 
mismatch  outputs  of  Cji(s)  and  Gi{s),  even  for  the  case  when  modeling  errors  and  compensator 
mismatches  exist  (as  would  be  the  case  in  a  nonlinear  simulation).  The  transfer  functions  for  these 
functions,  however,  cannot  be  predicted  by  linear  analysis. 

In  applying  this  two-signal  paradigm  to  the  control  of  the  nonlinear  plant,  Uii„{t)  drives 
the  signal  from  and  so  should  resemble  the  output  of  the  Gi(s)  compensator  for  the  linear  case. 
Simulations  show  that  this  is  indeed  the  case,  and  that  the  Limiter/Linear  Compensator  pair  induce 
a  linear-like  error  history  from  the  nonlinear  plant.  C(t)  then,  by  analogy,  should  also  resemble  the 
linear  output  time  history  of  of  Cr3(s). 

The  exact  form  and  magnitude  of  ^3(4)  is  unclear  when  errors  are  introduced  due  to  nonlin¬ 
earities,  so  the  contribution  of  the  second  linear  compensator  is  emulated  using  Fuzzy  Lo^c.  The 
exact  form  of  C{t)  can  then  be  determined  through  simulation.  C{t)  is  given  by: 

C{t)  =  X  *  rnin{nBi,zero{E),  (61) 

where  x  is  a  weighting  function  and  min{nEuZero{.E),liBi$PotixivX^))  **  It®  smaller  of  two  time- 
varying  Fuzzy  Membership  Functions.  This  function  mtn(»,  •)  is  one  realization  of  the  Fuzzy  AND 
operation,  hence  the  two  Fuzzy  rules  stated  above. 

By  specifying  the  weight  x  and  the  (x,<7’)  parameters  which  define  the  Fuzzy  Sets  EisZero 
and  EisPoaiiive,  C(t)  can  be  modified  to  damp  out  the  induced  oscillations.  The  precise  values 
necessary  must  be  determined  by  simulation,  and  are  functions  of  the  input  step  magnitude.  That 
is,  the  farther  the  reference  signal  is  from  the  quiescent  condition,  the  larger  the  dynamics  that  will 
be  induced  in  the  plant.  The  universe  of  discourse  for  EisPositive  must  be  rescaled  to  accommodate 
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larger  state  derivatives  and  the  overall  force  of  C’(t),  established  by  Xi  must  be  larger.  The  extent 
to  which  the  oscillations  are  damped  oat  is  highly  sensitive  to  changes  in  any  of  these  parameters. 

The  MBFLC,  then,  is  analogous  to  a  bank  of  two  linear  compensators,  weighted  at  their  Error 
inputs,  as  was  shown  in  Figure  3.12.  Notice,  however,  that  in  the  Figure  6.1  there  is  no  mechanism 
to  drive  the  Error  signal  into  the  Fuzzy  Limiter  to  zero  when  C(t)  activates,  as  was  the  case  for 
the  linear  bank.  No  mechanism  is  necessary  in  this  case  because  C{t)  only  becomes  active  when 
the  system  output  approaches  Yr-J  (E  -»  0).  Therefore,  C{t)  activates  as  the  signal  into  the  linear 
compensator  goes  to  zero  naturally.  A  properly  tuned  C(t)  will  keep  the  error  at  zero  and  allow 
the  linear  compensator  to  approach  steady-state. 

Full-envelope  performance  is  exhibited  over  a  prespeciiied  range  of  operation,  as  shown  in 
Figured  6.2.  The  MBFLC  clearly  forces  the  nonlinear  plant  to  exhibit  linear-like  behavior  and 
to  perform  within  the  specifications  established  by  the  linear  compensator.  Only  the  disturbance 
rejection  specification  is  violated  slightly  by  the  final  design.  This  eliminates  the  need  for  gain 
scheduling  within  the  envelope  of  operation  established  by  the  Fuzzy  Logic. 

6.S  Thesis  Conclusions 

The  expected  benefits  of  MBFLC  are  restated  here  from  Chapter  1: 

1.  Force  a  nonlinear  plant  to  to  exhibit  linear-like  performance  and  conform  to  linear  design 

specifications. 

2.  Incorporate  available  models  into  controller  structure. 

3.  Enhance  robustness  of  controller  in  the  face  of  unmodeled  uncertainties  (system  damage, 

noise). 

4.  Obtain  full  envelope  operation  without  the  need  for  gain  scheduling. 

5.  Introduce  some  a  priori  synthesis  and  analysis  capability  into  Fuzzy  Controller  Design. 

The  success  of  the  MBFLC  controller  at  meeting  these  objectives  is  examined  below. 

6.3.1  Force  a  nonlinear  plant  to  to  exhibit  linear-like  performance  and  obey  linear  design 
specifications.  The  simulation  responses  shown  in  Figure  6.2  demonstrate  that  the  MBFLC  is 
capable  of  inducing  a  linear-like  response  from  the  nonlinear  plant  such  that  it  falls  within  estab- 
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Figure  6.2  Closed-loop  response  of  Nonlinear  plant  for  step  inputs  of  Yr,)  =  0.6-|-0.035u_  i(t),  0.6-{- 
0.05i*_,(t),0.6-J-  0.075ti_i(t),  and  0.6  +  0.1u_,(<) 

lished  design  parameters.  Only  the  disturbance  rejection  specification  is  violated  by  the  MBFLC. 
In  this  respect,  the  MBFLC  is  a  success. 

When  compared  to  the  other  approaches  developed  in  Chapter  2,  the  MBFLC  controller  also 
shows  promise.  Figures  6.3  and  6.4  compare  the  responses  derived  from  applying  QFT  and  dynamic 
inversion  with  that  from  the  MBFLC.  These  Figures  show  that  the  MBFLC  delivers  performance 
approximately  equivalent  to  the  DI  compensator  for  the  nonlinear  plant  given  a  specified  region  of 
operation.  The  Fuzzy  Logic  Compensator  developed  in  Chapter  2  could  not  be  compared  with  the 
MBFLC  because  its  tuning  is  only  valid  for  one  input,  Yr,f  =  «-  i{i).  Any  other  input  or  operating 
point  would  required  redesign  of  the  FLC. 
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Figure  6.3  Closed-loop  responses  for  QFT-based  design  and  MBFLC  design 
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Figure  6.4  Closed-loop  responses  for  Dynamic  Inversion-based  design  and  MBFLC  design 
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6.3.2  Incorporate  available  models  into  controller  structure.  The  MBFLC  dearly  utilizes 
driving  data  derived  from  a  linear  compensator  rather  than  from  a  look-up  table,  as  for  the  FLC  of 
Chapter  2.  The  settling  term  C(t)  is  also  based  on  linear  simulations,  though  its  exact  value  must 
be  determined  through  simulation.  The  contribution  ultimately  made  by  the  linear  compensators, 
however,  is  small.  The  response  generated  by  the  F\izzy  Limiter/Linear  Compensator  does  reflect 
the  figures  of  merit  for  which  the  linear  compensator  is  designed,  but  this  is  largely  due  to  the 
tuning  of  the  Fuzzy  Sets. 

Several  experiments  from  Chapter  S  support  this  conclusion.  First,  it  was  demonstrated  that 
a  MBFLC  based  on  a  single  linear  compensator  was  actually  superior  in  some  ways  to  one  based  on 
a  bank  of  linear  compensators.  This  suggests  that  the  additional  linearizations  do  not  contribute 
to  the  control  of  the  nonlinear  plant,  even  in  the  neighborhoods  about  which  the  linearizations 
are  based.  Second,  it  was  shown  that  linear  robustness  does  not  translate  directly  into  nonlinear 
robustness.  Finally,  it  was  shown  that  large  errors  in  the  assumed  model  of  the  nonlienar  plant 
could  be  compensated  for  within  the  Fuzzy  Lo^c  with  very  little  d^adation  in  performance. 
These  results  indicate  that  the  contribution  of  the  linear  compensators  in  the  MBFLC  is  minimal. 

It  shotild  be  pointed  out,  however,  that  removal  of  all  linear  compensators  results  in  a  very 
undesirable  response.  The  the  Fuzzy  Logic  elements  can  produce  a  stable  response  without  the 
linear  compensator,  but  the  plant  output  is  highly  irregular  and  violates  the  established  performance 
boundaries. 


6.3.3  Enhance  robustness  of  controller  in  the  face  of  uncertainties.  The  MBFLC  archi¬ 
tecture  is  very  tolerant  of  modeling  error,  given  the  opportunity  to  retnne  the  Fuzzy  Sets.  Once 
retuned,  the  effect  of  mismodeling  is  exhibited  in  the  initial  overshoot  and  time  to  reach  steady- 
state  operations.  All  information  on  the  intended  design  characteristics  is  contained  in  the  linear 
compensator.  If  the  compensator  begins  to  significantly  differ  from  any  of  the  linearizations  of  the 
plant,  the  correct  characteristics  are  distorted. 
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If  no  retuning  is  allowed,  the  MBFLC  is  slightly  more  sentitive  to  plant  mismodding  than 
QFT  or  FLC.  However,  the  system  response  is  stable  throughout  the  range  of  variation  for  the 
parameter  a. 

6.S.4  Eliminate  the  need  for  gain  scheduling.  Once  tuned,  the  MBFLC  induces  linear-like 
behavior  throughout  the  specified  operating  envdope.  This  eliminates  the  need  for  gain  scheduling 
in  the  conventional  sense.  The  extent  to  which  the  envelope  can  be  expanded  given  a  single¬ 
compensator  MBFLC  is  unclear.  For  very  large  envelopes  of  operation,  the  bank  of  compensators 
may  prove  superior  to  the  single  compensator.  Farther  research  will  be  required  to  resolve  this 
issue. 


6.S.5  Introduce  some  a  priori  synthesis  and  analysis  capability.  The  linear  element  of  the 
MBFLC  can  be  designed  outside  of  the  simulation  environment.  In  this  respect,  some  degree  of 
a  priori  synthesis  capability  is  derived  through  the  MBFLC  approach.  The  final  design,  however, 
relies  heavily  on  toning  of  the  Fuzzy  Sets  to  establish  proper  weighting  and  granularity. 

Though  the  MBFLC  has  only  a  limited  a  priori  synthesis  characteristic,  it  still  has  a  well- 
defined  internal  structure  which  separates  it  from  other  FHizzy  Logic  based  approaches.  The  FLC 
from  Chapter  2  and  elsewhere  [10,  11,  12,  19]  is  based  only  on  empirical  tuning  data.  Even  when 
some  data  is  available,  as  from  the  DI  design,  the  look-op  table  based  on  this  data  may  wdl  be 
inadequate.  The  MBFLC  has  a  systematically  derived  internal  structure,  assuming  the  existence 
of  some  form  of  plant  model.  This  greatly  simplifies  controller  development  over  the  traditional 
FLC  approach.  FVirther,  the  MBFLC  approach  can  be  applied  to  systems  which  exhibit  high-order 
system  dynamics. 

6.4  Recommendations  for  Future  Research 

The  topics  touched  on  in  this  thesis  are  rich  for  future  research.  Recommended  areas  of  focus 
are  given  below: 
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6.4- 1  Analysis  of  MBFLC  error  sources.  Fkirther  research  into  the  actual  workings  of  the 
MBFLC  are  warranted.  The  claim  made  in  this  thesis  is  that  the  Fuzzy  Limiter/Linear  Compen¬ 
sator  induces  essentially  a  ’’linear”  response  in  the  nonlinear  plant.  This  argument  is  then  invoked 
to  justify  the  form  of  the  damping  signal  C(t)  and  to  explain  why  the  figures  of  merit  about  which 
the  linear  compensator  is  designed  are  ultimately  exhibited  by  the  nonlinear  plant.  The  validity  of 
these  claims  are  borne  out  in  simulation,  but  difficult  to  justify  theoretically. 

6.4.2  Self-Tuning  of  Fuzzy  Sets.  As  used  in  this  thesis,  FHizzy  Logic  is  simply  an  alternate 
way  to  express  a  nonlinear  look-up  table.  However,  the  true  utility  of  the  Fuzzy  Set  Theory 
is  reasoning  under  uncertainty,  and  learning  by  example  [33,  34].  A  F\izzy  learning  algorithm 
integrated  into  this  controller  could  greatly  reduce  tuning  requirements  and  enhance  robustness  in 
the  face  of  modeling  error. 

The  greatest  tuning  burden  is  to  determine  the  required  C{i)  to  damp  out  oscillations.  The 
forms  of  the  C(t)  responses  for  each  magnitude  step  input  are  currently  encapsulated  in  look-up 
tables.  These  tables  show  consistent  general  trends  but  are  difficult  to  estimate  precisely  prior  to 
simulation.  A  fuzzy  learning  algorithm  could  be  developed  to  perform  the  tuning  on-line  based  on 
the  perceived  performance  of  the  system.  Of  course,  this  would  require  some  sort  of  optimal  state 
trajectory  or  an  integrated  model  to  compare  agsunst  the  actual  system. 

If  learning  were  successful,  very  little  tuning  would  be  required  ad  hoc.  The  system  would  tune 
itself,  thereby  "getting  to  know”  its  host  system  through  actual  operations.  A  form  of  parameter 
identification  could  be  performed  by  optimizing  the  desired  performance  of  the  plant  over  time 
through  Fuzzy  or  neuro-fuzzy  unsupervised  learning.  Self  tuning  crmtiollers  is  an  area  of  current 
active  research  [24,  23,  35]. 

6.4. S  Development  of  equivalent  G2{s).  A  great  deal  of  the  complexity  of  the  MBFLC 
comes  from  the  emulation  of  the  second  controller  in  Fuzzy  Logic.  The  output  signal  C(f),  however, 
is  very  similar  to  that  produced  by  a  linear  controller.  Therefore,  C(f )  could  possibly  be  produced  by 
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an  equivalent  linear  network  developed  through  input/output  parameter  identification  techniques 

121]. 

This  would  reduce  the  role  of  Fu2zy  Logic  in  this  compensator  to  that  originally  envisioned  - 
providing  appropriate  weighting  to  a  bank  of  linear  filters.  The  bank  itsdf  would  consist  of  the  FViuy 
Limiter/Linear  Compensator  Pair  to  initiate  the  trajectory,  and  then  a  bank  of  linear  compensators 
each  corresponding  to  the  damping  signal  required  for  a  different  reference  input.  This  design  may 
not  perform  as  well  as  the  current  MBFLC,  doe  to  the  sensitivity  of  the  closed-loop  response  to 
small  changes  in  C(t). 

6.4.4  Compensation  of  plants  at  non-zero  initieU  conditions.  Banked  compensator  analysis 
demonstrated  the  difficulties  associated  with  developing  compensators  for  plants  which  are  away 
from  equilibrium  at  t,.  A  technique  should  be  explored,  perhaps  based  on  the  linear  banked  analysis 
of  Chapter  3,  to  develop  linear  compensators  for  the  control  of  non-quiescent  plants.  For  the  linear 
case,  this  could  be  pursued  by  examining  both  the  homogeneous  and  the  particular  solutions  for 
the  plants  in  question.  This  is  a  very  basic  problem  in  linear  control  and  should  be  explored. 

For  the  nonlinear  case,  the  results  6.4.3  could  be  utilized.  The  equivalent  linear  controllers 
developed  through  simulation  would,  in  effect,  be  designed  to  dissipate  state  energy  to  establish  a 
quiescent  condition.  This  type  of  ’’damping”  contitdler  could  have  applications  in  other  branches 
of  control  as  well. 

6.4.5  Polynomial  approximations  to  look-up  table  data.  If  adequate  equivalent  linear 
approximations  of  the  G,(s)  signals  cannot  be  developed,  the  look-up  tables  upon  which  the  MBFLC 
currently  relies  could  be  approximated  using  linear  or  piecewise  linear  equations.  This  could  lead 
to  a  more  linear  controller  which  provides  equivalent  contnd  as  the  MBFLC,  as  wdl  as  provide 
insights  into  on-line  retuning  opportunities. 

6.4.6  Model-FoUowing  Fuzzy  Logic  Controllers.  Model-following  hybrid  controDers  were 
shown  to  have  advantages  in  performance  over  any  of  the  designs  explored  in  this  research.  Unfor¬ 
tunately,  this  performance  is  acheived  at  the  cost  of  unrealistic  input  ”spikes”  into  the  nonlienar 
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plant.  There  are  many  des^  trade-offs  that  coold  be  made  to  adapt  a  model-frdlowing  approach 
for  actual  systems. 

The  hybrid  linear/Ftiszy  model  fcdlower,  for  example,  could  be  used  in  coigunction  with  other 
controllers.  Given  an  existing  controller,  the  hybrid  contndler  could  provide  inputs  only  when  the 
system  response  deviates  from  the  established  model.  In  this  fashion,  the  FLc  would  only  be  active 
when  absolutely  needed,  perhaps  even  through  auxiliary  contnd  surfaces. 

Nonlinear  Mathematical  Analysis.  By  approximating  all  of  the  Fuzzy  Logic  sets 
within  the  MBFLC  as  nonlinear  nth-order  functions,  nonlinear  stability  analysis  techniques  could  be 
applied  to  examine  the  stability  of  the  MBFLC.  This  technique  is  discussed  in  [36].  This  technique 
could  lend  mathematically  rigorous  insights  into  how  the  MBFLC  can  be  improved. 

6.4.8  Test  controller  in  realistic  environment.  This  thesis  largely  addresses  non-statistical 
uncertainty  due  to  nonlinearity.  The  problem  was  framed  in  Chapter  1,  however,  for  an  environment 
which  includes  parameter  uncertainty  and  noise.  These  elements  should  be  added  to  a  MBFLC 
simulation. 

The  current  MBFLC  design  does  not  demonstrate  an  exceptional  ability  to  cope  with  unknown 
parameters  without  retuning.  System  noise  essentially  adds  an  uncertainty  into  the  firing  of  the 
C(t)  signal,  which  has  varying  effects  depending  on  the  desired  reference  value  of  the  controller.  A 
desensitized  MBFLC  (by  self- tuning,  linear  approximation,  application  to  linear  plants,  or  through 
bang/bang  compensation)  should  be  considered  for  these  more  "real  world”  simulations. 

6.4.9  Implement  Fuzzy  Supervisor  Between  Uunit)  ond  C(f).  One  of  the  factors  con¬ 
tributing  to  tuning  difficulty  is  that  the  relative  contributions  of  Unmit)  and  C(t)  must  sum  to  the 
reference  value  as  steady  state.  An  additional  level  of  control  should  be  investigated  to  diminate 
the  need  for  this  precarious  balance. 

This  could  be  accomplished  by  reinserting  preweighting  terms,  implemented  through  either 
switches  or  a  Fuzzy  Supervisor.  In  this  case  the  linear  compensator  could  be  forced  to  reach 
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a  steady-state  value  aad  all  remaining  control  responsibility  would  be  1^  to  the  fksxy  Logic 
emulation. 

S.4.IO  Apply  MBFLC  Approach  to  Uncertain  Linear  Systems.  Though  the  MBFLC 
demonstrated  limited  robustness  in  the  face  of  modeling  error  for  nonlinear  plants,  robustness 
may  be  exhibited  when  considering  uncertain  linear  plants.  In  the  linear  case,  the  information 
contmned  in  the  linear  compensators  is  much  more  complete  than  in  the  nonlinear  scenario,  with 
fewer  "modes’*  of  failure. 

As  long  as  the  FHizzy  Logic  is  actually  producing  a  signal  to  emulate  a  linear  element,  the 
overall  design  will  exhibit  some  sensitivity  to  modeling  errors.  The  conventional  FLC  (from  Chapter 
2)  gets  its  robustness  by,  in  effect,  over-controlling  the  plant.  An  approach  similar  to  the  FHizzy 
bang/bang  controller  would  probably  be  more  effective  at  minimizing  the  sensitivity  of  the  overall 
system  to  undesired  dynamics.  There  are,  however,  other  practical  considerations  which  discourage 
the  implementation  of  bang/bang  controllers  on  physical  systems  [14]. 

6.5  Summary 

The  initial  hybrid  linear/Fuzzy  controller  concept  for  this  thesis  was  based  on  the  use  of 
multiple  linearizations  of  a  nonlinear  plant.  Linear  compensators  based  on  the  linearized  plants 
were  to  produce  plant  control  data,  assuming  the  plant  was,  in  fact,  in  the  re^on  where  the 
compensator  was  valid.  These  multiple  control  signals  would  then  be  blended  together  using  FVizzy 
Logic  to  produce  a  composite  control  signal  which  would  be  used  to  drive  the  nonlienar  plant.  The 
questions  addressed  in  this  thesis  were  tailored  to  show  conclusively  why  this  approach  would  or 
would  not  work. 

Unfortunately,  analysis  showed  that  a  Linear/Fuzzy  compensator  based  only  on  a  bank  of 
linear  compensators  was  not  feasible.  This  is  due  largely  to  the  small  re^ons  of  validitiy  for  the 
linearized  models  and  energy  considerations  within  the  plant/controller  system.  This  analysis  itsdf, 
however,  suggested  an  alternate  form  for  a  hybrid  Linear/Fbzzy  ^proach,  which  was  developed  in 
to  the  MBFLC  controller. 
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The  Model- Based  Fuzzy  Logic  controller  adds  a  new  dimension  in  capability  over  conventional 
Fuzzy  Logic  designs.  Current  FLCs  rely  on  heuristic  input/output  relationships  to  derive  appro¬ 
priate  control  commands.  The  MBFLC,  in  contrast,  utilizes  linear  models  and  linear  waveforms 
which  are  then  tuned  through  simulation  to  induce  the  desired  closed-loop  behavior.  The  controller 
design  presented  in  this  thesis  is  capable  of  driving  a  nonlinear  plant  to  any  point  within  a  prede¬ 
fined  envelope  of  operation  utilizing  only  Fuzzy  Logic  and  linearized  plant  models.  Its  performance 
is  comparable  to  a  controller  derived  using  deterministic  nonlinear  techniques.  The  MBFLC  also 
proved  to  be  approximately  as  robust  to  modeling  errors  as  the  other  control  approaches  considered. 

It  was  shown  that  limitations  in  the  current  MBFLC  can  be  overcome  by  model-following 
type  designs,  or  by  introducing  a  self-tuning  capability.  In  addition,  many  elements  of  the  MBFLC 
could  be  implemented  as  linear  networks,  simplifying  the  existing  architecture.  All  of  these  areas 
should  be  explored  in  follow-on  efforts. 

The  development  of  "intelligent  controllers”  is  contingent  on  the  successful  integration  of 
artificial  intelligence  into  conventional  control  methodologies.  This  thesis  demonstrates  that  the 
introduction  of  Fuzzy  sets  into  a  linear  controller  yields  new  capabilities  not  exhibited  by  linear  or 
Fuzzy  Logic-based  control  alone.  The  fusion  of  Fuzzy  Logic  techniques  and  conventional  control 
ideas  offers  many  possibilities  for  productive  future  research. 
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